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Abstract 



This paper concerns smooth transonic flows of Meyer type in finite de Laval nozzles, which are 
governed by an equation of mixed type with degeneracy and singularity at the sonic state. First we 
study the properties of sonic curves. For any C^ transonic flow of Meyer type, the set of exceptional 
points is shown to be a closed line segment (may be empty or only one point). Furthermore, it 
is proved that a flow with nonexceptional points is unstable for a C^ small perturbation in the 
^ ■ shape of the nozzle. Then we seek smooth transonic flows of Meyer type which satisfy physical 

^^ I boundary conditions and whose sonic points are exceptional. For such a flow, its sonic curve must 

be located at the throat of the nozzle and it is strongly singular in the sense that the sonic curve 

fvj I is a characteristic degenerate boundary in the subsonic-sonic region, while in the sonic-supersonic 

region all characteristics from sonic points coincide, which are the sonic curve and never approach 
the supersonic region. It is proved that there exists uniquely such a smooth transonic flow near 

^^ . the throat of the nozzle, whose acceleration is Lipschitz continuous, if the wall of the nozzle is 

sufficiently flat. 
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1 Introduction 

In this paper we study the compressible Euler system of steady isentropic and irrotational smooth 
transonic flows in a class of two-dimensional de Laval nozzles. Such problems arise naturally in physical 
experiments and engineering designs, and there exist large literatures on experiments and numerical 
simulations. However, a rigorous mathematical theory remains to be completed. 

A flow is called transonic if it contains both subsonic and supersonic regions, which are usually 
connected by shocks and sonic curves. Roughly speaking, transonic flows are governed by partial 
differential equations of mixed type; furthermore, shocks and sonic curves by which subsonic and 
supersonic regions are connected are free interfaces in general. Two kinds of transonic flow patterns 
have received much attention. One is transonic flows past a profile. Bers [3] showed that for two 
dimensional flows past an arbitrarily given profile, the whole flow field will be subsonic outside the 
profile if the freestream Mach number Mq within a certain interval < Mq < M < 1; furthermore, 
the maximum flow speed will tend to the sound speed as Mq -^ M. A natural question arises whether 
there also exists a smooth flow past the profile for M < Mq < 1. By constructing explicit solutions in 
the hodograph plane, one can get some smooth transonic fiows for some profiles and for some values 
of Mq ([!]). However, smooth transonic flows past a profile do not exist in general and are unstable 
even they exist according to Morawetz |161 \T7\ [TS] . There are few studies on transonic flows with 
shocks past a profile and almost no rigorous results are known. The other kind is transonic flows in a 
nozzle. Transonic flows in nozzles were investigated first by Taylor [19] and Meyer [15], where some 
special solutions were shown by using power series expansions. One usually considers, as we shall do 
in this paper, a de Laval nozzle which is symmetric and whose cross section decreases first and then 
increases. There are two types of smooth transonic flows in such a nozzle, named Taylor and Meyer 
types. In a transonic flow of Taylor type, there are two supersonic enclosures of the type encountered 
in a transonic flow past a profile, and smooth transonic flows of such type do not exist in general 
and are unstable with respect to small changes in the shape of the nozzle even they exist ([!]). In a 
transonic flow of Meyer type, the sonic curve extends from one wall of the nozzle to the other, which 
is believed to be located near the throat of the nozzle (where the cross section is smallest). Stable 
transonic flows of Meyer type should be subsonic upstream and supersonic downstream, while a flow 
with the reversing direction is unstable with respect to small changes in the shape of the nozzle. Most 
known smooth transonic flows in nozzles are solved as solutions to the governing equations and no 
boundary conditions are discussed ([!]). So, the nozzles cannot be given in advance. It is noticed 
that Kuz'min [11] solved the perturbation problems of accelerating smooth transonic flows with some 
structural assumptions in a class of nozzles by using the principle of contracting mappings. However, 
the physical meaning of the boundary conditions for these smooth transonic flows is not clear. There 
also exists a stable discontinuous transonic flow pattern in a nozzle called transonic shocks, which are 
supersonic upstream and turn to subsonic across shocks. We refer to [3 (U [71 [51 [12l [131 El 122] for 
the existence and stability of such transonic flows in flat and curved nozzles. Particularly, [12 1 113 1 [14] 
solved the transonic shock pattern described by Courant and Friedrichs: Given the appropriately large 
receiver pressure, if the upstream flow is still supersonic behind the throat of the nozzle, then at a 
certain place in the divergent part of the nozzle a shock front intervenes and the gas is compressed 
and slowed down to subsonic speed ([9]). For these transonic shocks, the shocks are free and the flows 
are away from the sonic state. 

The main goal of the present paper is to understand the behavior of transonic flows near the sonic 
state. Thus we investigate smooth transonic flows in a class of two-dimensional solid de Laval nozzles. 
As usual, it is assumed that the throat of the nozzle lies in the y-axis and the nozzle is symmetric 



with respect to the j;-axis. For convenience, we consider only the upper part of the nozzle due to its 
symmetry and the upper wall is given by 

^upw -y = fix), l^<x<l+, 

where /_ < < /+ and / G C^{[l^,l^]) satisfies 

' < 0, /_ < X < 0, 
fix) I =0, x = 0, 

>0, < X < /+. 

For the given inlet 

r,^:x = g{y), < y < /(/_), 

we seek a smooth transonic flow of Meyer type in the nozzle whose velocity vector is along the normal 
direction at the inlet. Since the flow is supersonic downstream, we choose the outlet 



out 



■.x = tiy), 0<y </(/+) 



(x,y) £n, 


(1.1) 


{x,y) eO, 


(1.2) 


(x,y) en, 


(1.3) 



as a free boundary where the velocity vector is along the normal of the outlet. Let f] be the domain 
bounded by Fupw, the x-axis, Fin and Fout- The transonic flow in fi satisfies the steady isentropic 
compressible Euler system: 

^(H + |(P-) = 0, 

— {P + pu^) + —{puv)=0, 

— {puv) + —{P + pv^) = 0, 
ox oy 

where {u,v), P and p represent the velocity, pressure and density of the flow, respectively. The flow 
is assumed to be isentropic so that P = P{p) is a smooth function. In particular, for a polytropic gas 
with the adiabatic exponent 7 > 1, 

P{p) = -p^ (1.4) 

7 

is the normalized pressure. Assume further that the flow is irrotational, i.e. 

— = — , (x,y)GO. (1.5) 

oy ox 

Then the density p is expressed in terms of the speed q according to the Bernoulli law ( [9] ) 

/ 2x / 7-1 o\V(7-l) ^- / 2 \l/2 ^ ^ 

p{q^) = (l- ^^g'j , q = Vn^T^, < g < q^ax = [-^) ■ (1-6) 

Summing up, the flow is governed by the system (jl.ip - (jl.6p . It is well known that the system (II. ip ^ 
(jl.6p can be transformed into the full potential equation 

div(p(|V(/.|2)V(^) = 0, {x,y)en, (1.7) 



where if is the velocity potential. 




Figure: the de Laval nozzle 

For a smooth transonic flow of Meyer type, the governing equation (jl.7p is elliptic in the subsonic 
region and hyperbolic in the supersonic region. It should be noted that such a flow is singular near 
the sonic curve in the sense that (jl.7p is degenerate elliptic in the subsonic-sonic region while non- 
strictly hyperbolic in the sonic-supersonic region. These factors cause the essential difficulties for 
the mathematical analysis. Another difficulty is that the sonic curve is usually believed to be free. 
So, understanding the behavior of sonic curves is not only important for physical applications and 
engineering designs but also one of the keys to study smooth transonic flows. Bers [21 H] studied the 
continuation of a flow across a sonic curve, where a subsonic-sonic flow was assumed to be given ahead. 
He found out that the existence and the structure of exceptional points play an important role in this 
continuation. The deflnition of exceptional points is referred to [4] (also § 2.2). For a C^ transonic flow, 
a sonic point is exceptional is equivalent to that the velocity vector is orthogonal to the sonic curve 
at this point. That is to say, exceptional points are characteristic degenerate in the subsonic-sonic 
region. If there is not any exceptional point, then, the flow can be continued in a unique way across 
the sonic curve as a supersonic flow without discontinuity. Whereas, if there is a unique exceptional 
point, the flow will be uniquely determined only in a strict subset of the same domain; furthermore, 
if the flow can be continued at all into this excluded region, this continuation will not be unique. 
Bers also mentioned in [1] that it would be interesting to know whether exceptional points are always 
isolated. Here, the supersonic flow is solved as a Cauchy problem from the sonic curve and there is 
no prescription of boundary conditions on the walls. For a smooth transonic flow of Meyer type in 
a nozzle, the supersonic flow satisfles not a Cauchy problem but an initial-boundary value problem. 
By some precise analysis on sonic curves, it is shown in the paper that for any C^ transonic flow, the 
set of exceptional points is a closed line segment (may be empty or only one point). Furthermore, 
exceptional points are strongly singular in the sonic-supersonic region in the sense that there are two 
different characteristics from each nonexceptional point in the nozzle, while all characteristics from 
interior exceptional points coincide and they never approach the supersonic region locally. Then, it 
is proved that a flow with nonexceptional points is unstable for a C^ small perturbation on the wall 
(even if the wall is still smooth). This instability is weak since it is unknown whether the flow is 
unstable if the wall is perturbed in C^ or other smooth spaces. So we seek a smooth transonic flow 
of Meyer type whose sonic points are exceptional. For such a flow, its sonic curve must be located at 
the throat of the nozzle and the potential on its sonic curve equals identically to a constant. Another 
motivation arises from our early paper [20], where the structural stability problem of a symmetric 



continuous subsonic-sonic flow in a convergent nozzle with straight solid walls was proved and the 
sonic curve of the continuous subsonic-sonic flow is shown to be a free boundary where the potential 
equals identically to a constant. In |20] the continuous subsonic-sonic flow is singular in the sense that 
while the speed is continuous yet the acceleration blows up at the sonic state. This singularity arises 
from the geometry of the nozzle and it seems to vanish if the sonic curve is located at the throat of a de 
Laval nozzle. It is noted that for a smooth transonic flow of Meyer type whose sonic curve is located 
at the throat of the nozzle, the potential at the sonic curve, which equals identically a constant, is free 
although the location of the sonic curve in the nozzle is known. 

For smooth transonic flows, the governing equation is elliptic in the subsonic region and hyperbolic 
in the supersonic region; furthermore, it is degenerate and singular at the sonic state. So, to study a 
smooth transonic flow problem, one must determine and control precisely the speed of the flow near 
the sonic state, which plays an essential role in the mathematical analysis. In p.7p . the speed of 
the flow is the absolute value of the gradient of a solution. So, generally speaking, it is very hard to 
estimate precisely the speed of the flow near the sonic state in the physical plane even if the location 
of the sonic curve is known. It is more convenient to study smooth transonic flows in the potential 
plane, where the speed of the flow is a solution to the governing equation. After confirming the sonic 
curve to be a free interface where the potential equals identically to a constant, we decompose the 
smooth transonic flow problem into a smooth subsonic-sonic flow problem with free boundary and a 
smooth sonic-supersonic flow problem with fixed boundary in the potential plane, which can be solved 
separately. For the smooth subsonic-sonic flow problem, we encounter the two main difficulties in the 
study on the continuous subsonic-sonic flow problem in [20] . One is that the problem is a characteristic 
degenerate free boundary problem and the degeneracy occurs just at the free boundary. The other is 
that at the inlet we should prescribe a Neumann boundary condition instead of a Robin one with which 
the problem may be ill-posed. However, besides these difficulties, there is an additional disadvantage 
in this study. Different from the problem in |20] . there are no background solutions, which play an 
important role to determine and control the rate of the flow from the subsonic region to the sonic 
state. Furthermore, the boundary conditions at the inlet and the wall are nonlinear, nonlocal and 
implicit. In the paper, we try to seek a subsonic-sonic flow which tends uniformly to the sonic state 
along different streamlines although there are no background solutions. To this end, we assume that 
the convergent part of the nozzle slopes so gently that 

/"(x) = o(x2), x^O^. (1.8) 

Indeed, (jl.Sp is necessary to guarantee that the change of the speed along the stream direction is 
infinitesimal of higher order than the one along the potential direction near the sonic state. By a 
fixed point argument with many very precise elliptic estimates, we obtain a subsonic-sonic flow which 
tends uniformly to the sonic state along different streamlines. For the smooth sonic-supersonic flow 
problem, the governing equation is a non-strictly hyperbolic equation with strong singularity at the 
sonic curve, where its two eigenvalues coincide and the eigenvector space reduces to a one-dimensional 
space. Furthermore, the singularity at the sonic curve is so strong that all characteristics from sonic 
points coincide, which are the sonic curve and never approach the supersonic region. And the boundary 
condition on the wall is nonlinear, nonlocal and implicit. As mentioned above, a crucial step to solve 
the problem is how to determine and control the rate of the flow from the sonic state to supersonic 
region. We seek a sonic-supersonic flow which moves uniformly away from the sonic state along 
different streamlines in the paper. For the sonic-supersonic flow in the divergent part of the de Laval 
nozzle, 

fix) = o{x^), X ^ 0+ (L9) 



is necessary to guarantee that the change of the speed along the stream direction is infinitesimal of 
higher order than the one along the potential direction near the sonic state. We solve this sonic- 
supersonic flow problem under the assumption (jl.9p by a fixed point argument and the method of 
characteristics. Due to the strong singularity, the computations of the flow near the sonic curve are 
quite complicated. Through some precise calculations and optimal estimates, we are able to get a 
desired sonic-supersonic flow. 

In the present paper, by solving a smooth subsonic-sonic flow problem with free boundary and a 
smooth sonic-supersonic flow problem with fixed boundary separately, we get a smooth transonic flow 
of Meyer type in the de Laval nozzle il, whose velocity vector is along the normal direction at the inlet 
and which satisfies the slip condition on the wall, under the assumptions (jl.Sp . (|1.9p and |/±| being 
sufficiently small. Different from the examples of smooth transonic flows of Meyer type by using power 
series expansions, where the exceptional point is isolated, the sonic curve of the smooth transonic flow 
in this paper is located at the throat and each sonic point is exceptional. Thus, this transonic flow 
pattern is strongly singular in the sense that the sonic curve is a characteristic degenerate boundary 
in the subsonic-sonic region, while in the sonic-supersonic region all characteristics from sonic points 
coincide, which are the sonic curve and never approach the supersonic region. It is surprising that 
there is a smooth transonic flow for this pattern with so strong singularity. Indeed, we get a smooth 
transonic flow in the sense that the acceleration is Lipschitz continuous. Furthermore, the transonic 
flow of this pattern is also shown to be unique. It is noted that (jl.Sp and (|1.9p are almost necessary for 
a C^ transonic flow whose sonic curve is located at the throat of the nozzle. Moreover, the assumption 
\l±\ being sufficiently small can be relaxed. In the convergent part and the divergent part of the nozzle, 
we get a smooth subsonic-sonic flow and a smooth sonic-supersonic flow, respectively. Using the same 
methods, one can get a smooth transonic flow of Meyer type in a long de Laval nozzle if it possesses 
the same properties near the throat ( (jl.Sp and (II. 9p . which are almost necessary) and its walls slope 
gently away from the throat. More generally, we believe that there also exist smooth transonic flows 
for general flnite de Laval nozzles satisfying (II. Sp and (11. 9p . In the upstream, by more complicated 
elliptic estimates than the ones in the present paper, we can get a subsonic-sonic flow for a very general 
convergent nozzle. In the downstream, the flow can be extended for a class of long divergent nozzles 
by solving an initial-boundary value problem of a strictly hyperbolic equation. These questions will 
be dealt with in our forthcoming studies. 

The paper is arranged as follows. In § 2 we analyze the structure and the location of the sonic 
curve for C^ transonic flows and show the instability of transonic flows with nonexceptional points, 
which motives us to seek a smooth transonic flow of Meyer type whose sonic points are exceptional. We 
formulate the problem of smooth transonic flows and state the main results (existence and uniqueness) 
of the paper in § 3. The smooth transonic flow problem can be decomposed into a smooth subsonic- 
sonic flow problem with free boundary and a smooth sonic-supersonic flow problem with fixed boundary 
in the potential plane, which are investigated in § 4 and § 5, respectively. 

2 Sonic curves and instability of transonic flows with nonexceptional 
points 

In this section, we analyze the structure and the location of sonic curves for C^ transonic flows. And 
we always assume that the sonic curve belongs to both the boundary of the subsonic region and the 
boundary of the supersonic region. 



2.1 Governing equations 

Let us first rewrite tlie system (jl.ip - (jl.6p in the physical plane as the Chaplygin equations in the 
potential plane. 

Since the conservation of momentum, ()1.2p and (11.30 . can be derived from (II. ip . (I1.4p ~ (ll.6p . the 
flow is just governed by the conservation of mass and the condition of irrotationality: 

d , . d , , du dv , . 

-J,u) + -(,.) = 0, ^-^ = 0, (2.1) 

where p is related to the velocity {u,v) by the Bernoulli law (jl.6p . The sound speed c is defined as 

c =P{p)= p^ =1 —q . 

At the sonic state, the speed is 

1/2 



^7 + 1' 

which is called the critical speed in the sense that the flow is subsonic {q < c) when < g < c*, sonic 
{q = c) when <? = c* and supersonic {q > c) when c* < g < qmax- 

Define a velocity potential (p and a stream function ■0, respectively, by 

dip dip dip dip 

ax ay dx dy 

which are 

d^ a ^'P ■ a 9^ ■ a ^^ a 

-—-=qcosU, -—- = qsint), —— = —pqsmO, —-= pqcosti 

dx dy dx dy 

in terms of polar coordinators in the velocity space, where 6, which is called a flow angle, is the angle 
of the velocity inclination to the x-axis. Direct calculations show that the system ()2.ip can be reduced 
to the Chaplygin equations: 

de ^ p{q^) + 2q'p'{q^) dq _^ 1^ _ J_^ ^ (2 2) 

dtp qp^ifp) dtp ^ qdip p{q^) dp> 

in the potential-stream coordinates {ip^ip). Note that 

5((/3, ip) dip dif) dip dip 2 

5(x, y) dx dy dy dx 

So the coordinates transformation between the two coordinate systems is valid at least in the absence 
of stagnation points. Eliminating 6 from (j2.2p yields the following second-order quasilinear equation 

d^Ajq) d^Bjq) _ 

~^^ + ~^^ ~ ^' ^^■'^' 

where 

„2\ I 0,2^/^,2 



A{q) 



/ 2^2^ ^ ^' ^^'^) = / ^^^^^^' < g < (?„ 



Here, B{-) is strictly increasing in {0,qmax), while A{-) is strictly increasing in (0, c*] and strictly 
decreasing in [c^,qmax)- The inverse function oi B{-) is denoted by B~^{-), while the inverse functions 
of A{-) lying in (0, c*] and [c^,,qmax) are denoted by AZ (•) and A'^ (•), respectively. 

It can be checked easily that both (I1.7P and (I2.3P are elliptic in the subsonic region ((7 < c*) and 
hyperbolic in the supersonic region (q > c*), while degenerate and singular at the sonic state (g = c^=). 
Therefore, the governing equations (11.7P and (|2.3p for transonic flows are both second-order mixed 
type quasilinear equations and are both degenerate and singular at sonic states. 

2.2 Exceptional points in the physical plane 

We begin with a description of sonic curves in [1]. Let 5 be a sonic curve of a C^ transonic flow. The 
positive direction on <S is defined by requiring that, if one moves along S in this direction, the subsonic 
region is located on the left. Then, Og, the derivative of 9 with respect to the arc length on S, satisfies 

es = -'^^^, (2.4) 

where 1/ is the unit normal of S pointing into the supersonic region and i? is the angle between the 
velocity vector and u. 




{u,v) 
subsonic I supersonic 



Figure: the orientation of the sonic curve 
It follows from ([23]) that 
Lemma 2.1 9 is nonincreasing along S. 

As in ^, points, where Og = 0, are called exceptional. According to (|2.4p . a point is exceptional if 

dq 
and only if at this point either 7— = or the velocity vector is orthogonal to the sonic curve. Indeed, 

ou 
dq 
at a point on S, if — — = 0, then the velocity vector is orthogonal to the sonic curve owing to the 

ou 
following lemma by Gilbarg and Shiffman |10j . 




Figure of Lemma [27 

Lemma 2.2 Assume that q = c^ at a point P on the circumference of a circle in whose interior 
q < c.^. If the velocity vector is not along the direction from the center of the circle to P, then at P 
the derivative of q along this direction is positive. 



Proof. The case when the circle is the interior of a subsonic-sonic region was considered in [10] . 
The authors derived a hnear elUptic equation which is degenerate at the sonic state from ()1.7p . A 
comparison principle holds for this degenerate equation, which can be proved by the same argument 
for uniformly elliptic equations. Then, the lemma follows from the same proof of the Hopf lemma and 
the auxiliary function has been constructed in [1^. Here, it is used that the velocity vector is not along 
the direction from the center of the circle to P, which shows that this direction is not characteristic 
for the degenerate equation. The proof is standard and thus omitted. D 

Therefore, one has 

Proposition 2.1 For a C^ transonic flow, a point at the sonic curve is exceptional if and only if the 
velocity vector is orthogonal to the sonic curve at this point. 

This proposition shows 

Remark 2.1 For a C^ transonic flow, exceptional points are characteristic degenerate in the subsonic- 
sonic region. 

Exceptional points are regarded to be singular in the following sense ([!]). Assume that a subsonic 
flow is defined in some domain whose boundary contains a smooth sonic curve S. If there is no 
exceptional point on S, then, the flow can be continued in a unique way across 5 as a supersonic flow 
without discontinuity. The flow will be determined in some neighborhood of S contained between the 
two Mach lines from the two endpoints of S. However, if there is a unique exceptional point on S, 
the flow will be uniquely determined only in a subset of the same neighborhood, where the points 
between the two characteristics from the exceptional point are excluded; furthermore, if the flow can 
be continued at all into this excluded region, this continuation will not be unique. In this section, 
we will show that the set consisting of exceptional points for any C^ transonic flow of Meyer type 
is a closed line segment (may be empty or only one point), while there is no exceptional point for 
any C^ transonic flow of Taylor type. Furthermore, the singularity of exceptional points in the sonic- 
supersonic region is so strong that there are two different characteristics from each nonexceptional 
point in the nozzle, while all characteristics from interior exceptional points coincide and they never 
approach the supersonic region locally. In the present paper, we get a smooth transonic flow of Meyer 

10 



type whose sonic points are all exceptional. This transonic flow pattern is strongly singular in the 
sense that the sonic curve is a characteristic degenerate boundary in the subsonic-sonic region, while 
in the sonic-supersonic region all characteristics from sonic points coincide, which are the sonic curve 
and never approach the supersonic region. It is surprising that there is a smooth transonic flow for 
this pattern with so strong singularity. Indeed, we get a smooth transonic flow in the sense that the 
acceleration is Lipschitz continuous. This transonic flow pattern also answers the problem mentioned 
in [3] that it would be interesting to know whether exceptional points are always isolated. 

2.3 Exceptional points in the potential plane 

It turns out to be more convenient to analyze exceptional points in the potential plane. Transforming 
Proposition 12.11 from the Cartesian coordinates to the potential-stream coordinates, one can get 

Proposition 2.2 For a C^ transonic flow in the potential plane, a sonic point is exceptional if and 

dq 
only if — — = at this point, 
dip 

Proof. On the sonic curve, it holds 

dq dq f dq . 2^ dq\ f dq 2^ 9q\ ,2x2^9 

Therefore, at a sonic point the velocity vector is orthogonal to the sonic curve in the physical plane if 

Qq 
and only if --— = at this point. Then, the proposition follows from Proposition 12. 1[ D 

dtp 




Figure of Lemma [X 

Similar to the proof of Lemma 12.21 it can be shown that 

Lemma 2.3 Consider a C"^ flow in the potential plane. Assume that q = c^ at a point P on the 

circumference of a circle in whose interior q < c^. If the direction from, the center of the circle to P 

dq 

is not parallel to the Lp-axis, then — — > at P, where fi is any direction which forms an acute angle 

dfj, 
with the direction from, the center of the circle to P. 

Proof. Without loss of generality, we assume that P = ((/JqjV'o) is the unique sonic point on the 
circumference of the circle. Otherwise, we can consider the flow in a smaller circle which intersects 
the original circle at P. For convenience, it is assumed that the circle is centered at the origin with 

11 



radius r. Since the direction from the center of the circle to P is not parallel to the ip-axis, "00 7^ 0. 
Set 



with positive constants e (< c*/2) and M to be determined, where 

Bi = {(V9, V^) : (^2 + ^2 ^ ^2|^ ^^ ^ |(^^^) . (^ _ ^^)2 ^ (^ _ ^^)2 ^ 1^2|^ 

Note that 



y <g(c^,V') <c*, {^,^)eBir\B2 



and 



d'^Ajq) d'^B{q) 

<A'{q{ip,^P))e{2M - 4MV^)e"*^^^'+'^'^ + B'{q{ip,iP))e{2M - 4MV^)e"*'^^^'+'^') 

< - Me{B'{q{^,^))iMi;^, - 2) - 2A'{q{^,^)))e-^'^'^"+^"^ 

< - Me({Mtp'^ - 2) min B'{s) - 2 max A'(s))e"*''('^'+'^'\ (v?,'^) G -Bi n B2. 

Therefore, there exists a sufficiently large M > such that g is a supersolution to (j2.3p in Si n i?2 for 
each e G (0, c*/2]. Choose e € (0, c*/2] so small that 

q{ip,^p) <q{ip,ip), {ip,^p) eBindB2. 

Then, it follows from a comparison principle that 

qi(p,ip) <qi(p,'ilj), {ip,^;) e Bir\B2, 

which leads to _ 



Here, the comparison principle is not the classical one since (j2.3p is degenerate at the sonic state and 
it can be proved similar to Proposition 3.1 in ^0]. D 



As a corollary of Proposition 12.21 and Lemma 12.31 it holds that 

Lemma 2.4 For a C^ transonic flow in the potential plane, if an interior point in the nozzle is 
exceptional, then the sonic curve near this point is a graph of function with respect to ip. 

Qq 
Proof. If the lemma is false, then — — = at this exceptional point. Note that Proposition 12.21 

Oif 

dq . . 

implies — — = at this exceptional point. So the derivative of q at this exceptional point along any 

oil) 

direction is zero. Since this exceptional point is an interior point in the nozzle, Lemma [2. 31 shows that 
the normal of the sonic curve at this exceptional point is parallel to the 93-axis, which yields that the 
sonic curve near this point is a graph of function with respect to '(/'■ D 

If an exceptional point is on the wall of the nozzle, one can get that 
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Lemma 2.5 Consider a C^ transonic flow in the potential plane with the sonic curve given by 

S■.^ = Sl{s),^l^ = S2{s),{S[{s))^ + {S!,{s))^>0, < s < 1, 5i, 52 G ^^([0, 1]). 
Assume that «S'2(0) = inf ^2 < sup 52 = 52(1) and the subsonic region is located on the left of S. 

(0,1) (0,1) 

(i) If S[{0) > and (5i(0), 52(0)) is exceptional, then S near (5i(0), 52(0)) is a graph of function 
with respect to ip. 

(ii) If S[{1) < and (5i(l), 52(1)) is exceptional, then 5 near (5i(l), 52(1)) is a graph of function 
with respect to ip. 

Proof. We will prove (ii) only since (i) can be dealt with similarly. If 5j(l) = 0, then 52(1) 7^ 
and the conclusion of the lemma is obvious. Therefore, if (ii) is false, then 5[(1) < and 5 near 
(5i(l), 52(1)) is a graph of function with respect to c/9, which is denoted by 

V' = /i(v?), 5i(l)<99<5i(l)+r (r>0) 

with 

/i'(5i(l)) = 0. (2.5) 

It follows from 52(1) = sup 52 that 

(0,1) 

h{^)<h{Si{l)), 5i(l)<99<5i(l) + T, 

which, together with (12. 5p . leads to 

-00 < /i"(5i(l)) < 0. 

If —00 < /i"(5i(l)) < 0, then there exists a circle whose interior is located at the subsonic region 

Qq 
and which intersects with 5 at (5i(l), 52(1)). Thus Lemma [2^ shows — — (5i(l), 52(1)) > 0, which 

dip 

contradicts (5i(l), 52(1)) is exceptional by Proposition 12.21 Therefore, 

h"[Si{l)) = -oo. (2.6) 

Due to (|2.5p and (j2.6p . there exists a number e S (0,r] such that 

/i'((^)<0, 5i(l)<^<5i(l)+e, 

which shows that 5 near (5i(l), 52(1)) is a graph of function with respect to V'- D 

Remark 2.2 As will he shown later (see Remark \2.5\) . neither S[{0) > in (i) nor 5[(1) < in (ii) 
can occur. 

Remark 2.3 In Lem.ma \2.5\ if either S[{0) < in (i) or S[{1) > in (ii), then the similar proof is 
invalid and it is unknown whether the conclusion of the lemma holds. At this moment these two cases 
cannot be excluded for transonic flows in a nozzle. 
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2.4 Behavior of sonic curves near exceptional points 

Consider a C^ transonic flow with the sonic curve 5 in the potential plane. To study the behavior of 
the sonic curve near exceptional points, we start with a geometric property of the sonic curve. 



v- 



X 



subsonic I supersonic 
Figure of Lemma 12.61 (i) 



^ 



supersonic subsonic 



"f 



Figure of Lemma 12.61 (ii) 



Lemma 2.6 Assume that Sq : ip = ho{ip) (■(/'i <^< V'2i ipi < 4^2) is a portion of S. 



then 



then 



(i) If the flow is subsonic in {{(fjip) : fJ-i < f < /io(V')i V'l ^ "ip ^ ^2} with a number ^li < min Hq, 

ho{ip) > min{/io(V'i), ho{ip2)}, ipi < ip < ip2- 

(ii) If the flow is subsonic in {((/9, -0) : hQ{Tp) < (/? < /i2) "01 ^ "0 ^ "02} with a number ^2 > rnax ho, 

bPi,'>p2] 



hoiip) < ma■K{ho{'^pl),ho{^p2)}, ipi < ip < ^^2 



Proof. We prove (i) only and the proof of (ii) is similar. Otherwise, there exists ■00 £ (V'i)'02) 
such that 



ho{ipo) = min /i < min{/io('0i), ^o(V'2)}- 

['/'l,V'2] 



Set 



q{ip,i;) = Aj{6{ip - /lo(V'o))), ^J'l <^ < /lo(V'o), -01 < V' < ^2, 
where (5 is a positive number so small that 

sup A(g(/ii,-)) < -(5(/io(V'o) -pti), sup {A{q{-,tl;i)),A{q{-,'4)2))} <-S{ho{tl;o)- ^ii). 

(■>Pi,it'2) {fiiMii'o)) 

Then, it follows from a comparison principle that 

qiV'^ip) < q{^,il^) = -4l^(5(y? - /lo(V'o))), f^i < ^ < /lo(V'o), V'l < "0 < "02, 

dq 
which leads to -^— (/io(V'o))'0o) = +00. Here, the comparison principle is not the classical one since 

D 



2.3p is degenerate at the sonic state and it can be proved similar to Proposition 3.1 in |20j . 
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v- 



subsonic |(/io(V'o),V'o) 



Figure of the proof of Lemma 



The following lemma and proposition describe the behavior of the sonic curve near an exceptional 
point. 



^ 



r(/io(V'2),^2) 

subsonic/: supersonic 

subsonic \i supersonic 

I(/io(V;i),^i) . 



9? 



Figure of Lemma 12.71 



Lemma 2.7 Assume that Sq : ip = /lo(V') (V'l < ^ ^ "02; V'l < "02) is a portion of S and the flow is 
subsonic on the left of Sq . 
dq 
(i) // — --(/io('0i)) V'l) = and h'Q^tjJi) = 0, then there exists a positive constant e (< "02 — V^i) such 

that 



ho{i)) < /lo(V'i), i^i<^ <4'i+e. 



dq 



(ii) If TT— (/io(02)5 0^2) = and h'r,(tjj2) = 0, then there exists a positive constant e (< 4'2 — V'l) 
oil) 

such that 

ho{ip) < /io(V'2), ip2-e <iIj < 02- 

Proof. We prove (ii) only and (i) can be proved similarly. Otherwise, Lemma 12.61 shows that 
there exists ipo € [ipi,ip2) such that 



ho{iP) > /io(V'2), K{iP) < 0, ipo<ip < -02. 
Furthermore, ()2.7p . Proposition 12.21 and Lemma 12.31 lead to 

dq 



^(/ioW,V')>o, Vo<V'<V'2 



(2.7) 



(2. 



and 



\ip & (V'Oi V'2) : T^(^o(V');^) > > is dense in any left neighberhood of 'i/'2- (2-9) 
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Owing to ()2.3p , one gets that 



B'{c,)^{hoW,^) = -A"{c,)(^^{ho{i^),i;)) -B"{c.)(^^{hoW,^)) > 0, Vo < V' < ^^2. 



Therefore, 



^V' 



9g 

9(^ 



(/l0(V'2),V'2)=0. 






Otherwise, (|2.10p yields 7772(^0(^2)5^2) > and hence 



d ( dq 



di/j \di{j 



ihoW,i^) 






d'^q 



(/l0(V'2),V'2)/ld(V'2) + ^(/io(V'2),V'2) > 0, 



Qq 
which contradicts Tr-(/io(V'2),''/'2) = and (12.8 
dip 



(2.10) 



(2.11) 



^ 



■\\(/lo(V'2),'02) 

subsonic ^ (^\ A supersonic 
ho - ^\ \ \/io 
"h 

Figure of the proof of Lemma 12.71 



Assume that ^ > 0, < r < ^2 — V'o and /i G C^(['(/'2 — t, 1^2]) satisfy 

/i'(V2) = 0, /lo(^) - ^ < /i(V) < /io(V'), /i"(V')>0, V2-T<^<V2, 

g(C(5,'0) < C*, /lo(V') -^f < 95 < /l(^), ^2 -T <'(/'< ^2, 



c^V' 



(y?, ■02 - t) > 0, /lo(V'2 - t) - ^i < 99 < /lo(V'2 - t). 



(2.12) 
(2.13) 

(2.14) 



Set G = {{Lp,ip) : ho{'tp) — ^ < ip < h{ip), 0^2 — t < V' < V'2} and 

q{ip, V^) = c* + (5(0 - V2)(^(0) - y'), (f/?, -0) e G, 

where 5, which depends only on fi, r, sup q{ho{) — n,-), inf -^(•,'02 — "J") and 

{i>2-T,ip2) {hoiip2-r)-fJ.,ho{'^2-r)) Oip 

ll^'llL°°((V'2--r.i/'2))5 is a positive number so small that 

q{ho{ip) - f^,ip) <q{ho{ip) - fi,ip), 02-T<'0<02 (2.15) 

and 

dq do 

— ((/?, '02 -r) > —(v:), 02 -r), ho{Tp2 - t) - fi < ip < h{ip2 - t) . (2.16) 
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Noting 

0(99,^) = 0, ^{^A) = 5(2/i'W + (V' - i^2)h"m < 0, {^,^P) G G 

due to (12.121) . one gets that g is a supersolution to ()2.3p in G. Since q is strictly subsonic in G from 
(j2.13p . the classical comparison principle holds. Then, it follows from (|2.15p and (j2.16p that 

q{^,i^)<c,+6{^-'^2){hii;)-^), /io(V')-M<'/5</i(^), ^2-T< V<V'2. (2.17) 

We now derive a contradiction. Owing to (|2.7p - (|2.9p . there exist fi>0, 0<T<ip2 — tpo, 
{/ifc}r=i C C^{[iIj2 - T, V2]) and K}~=i C (0,t) satisfying 

ho{i;)-fi<hk{i;)<hk+i{i')<hoW, ^2-T<V<V'2, A; = l,2,---, (2.18) 

/i',(^2) = 0, h'liiP) > 0, V2 - r < V < ^-2, fc = 1, 2, • • • , (2.19) 

q{ip,ip)<c^, hoiip) - fi<ip < ho{Tp), ip2-T <ip <ip2, (2.20) 

|^(c^,V'2-r)>0, /io(V'2-r)-/i<v9</io(V'2-r), (2.21) 

lim Tn = and lim hk{ip2 — Tn) = ^0(^2 — Tn) for n = 1, 2, • • • . (2.22) 

Note that (|2.18p - (j2.2ip imply that /i^ satisfies (|2.12p - (|2.14p for each positive integer k. It follows from 
(12171) that 

q{ip, ip) <c^ + 6{ip - il)2){hk{ip) - if), ho{ip) - fi < if < hk{tp), ^^2 - r < ip < ip2, k = l,2,--- , 
where (5 > is independent of k. This, together with (j2.22p and (j2.1ip . leads to 

Qq 

^-(^o(V'2 -Tn),'lp2 -Tn) > ^Tn, n= 1,2, ••• 
dip 



and thus 

d^q 



(/io(^2),V'2)<-<5. (2.23) 



difdip 
Due to (I2.10p and (12.23p . there exists a number e G (0, r) such that 

A.(^^{hoW,i^)) = ^{h,{^li,)h',{^l:) + 0(/lo(V'),^) > -^/i'o(V), V'2 - £ < V^ < ^2, 

which leads to 

||(/10(^),^) - ||(/io(^2),V'2) < -^(/lo(^) -/10(^2)), i^2-e<ij<iJ2. 
Qn 

This, together with Tjy(/io(^2),^2) = 0, dMI and ([2S]), yields that 

— (/io(^),V) =0, 11^2 - e < il) < 11)2, 
which contradicts (j2.9p . D 
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Let us remove the restrictions /ig('0i) = and /io(v^2) = in Lemma |2.7[ 






subsonic /supersonic 

subsonicVsup er sonic 

(/io(^i),^i) V^ 
Figure of Proposition [27 



Proposition 2.3 Assume that Sq : if = hQ{'ip) {^i <ip< 1^2, "01 < ^^2) is 0. portion of S and the flow 
is subsonic on the left of Sq . 
dq 

(i) // — — (/io('i/'i),'0i) = 0, then there exists a positive constant e {< Tp2 — V'l) such that 
dip 

/lo(V') < ^o(V'i), V'l < V' < "01 + £• 

dq 

(ii) // T7--(/io('02)) "02) = 0, then there exists a positive constant e (< ip2 — V'l) such that 
dip 

ho{ip) < /io(V'2), -02 - e < ^ < ^2- 

Proof. We prove (ii) only and the proof of (i) is similar. If the lemma is false, then Lemma 12.71 
implies that — 00 < /iq(^2) < and thus 

p-iho{ip2),i^2) = 0. (2.24) 

Since the flow is subsonic on the left of {ho{ip2),tp2), (|2.24p implies 

|4(^o(V'2),V'2)<0. (2.25) 

Choose a positive number r (< '(/'2 — V'l) such that 

h'oiiP) < 0, V2 - T < V < V'2. (2.26) 

Qq 
It follows from Lemma [2731 and (j2.26p that -jr—{ho{ip2 — T),i^2 — t) > 0. Let fi he a positive number 

so small that \{ip,ip) : ho^ip) — fJ. < f < hQ{ip),Tp2 — t < ip < ip2\ belongs to the subsonic region and 

o 1 ^ 

— (99,-02 -t)> 2 ^(^0(^2 -t),V'2 -t), /lo(V'2 -t) - fl<ip < /lo(V'2 -t). 

We will use a similar method as in Lemma [2.71 to complete the proof. Set 

^ ={h € C^{[iP2 - T, V2]) : hoi^P) -fi< h{iP) < hoiiP), h'iiP) < 
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and {11)2 - '4^)h"{'4') > 2h'{'ip) for each V2 - t" < V" < V'2}- 
For any < e < r, if /i € C'^{[ip2 — £, ^2)) satisfies 

{1P2 - ^)h"{ip) < 2/i'(^) < 0, V2 - e < V' < ^^2, 
then h must be unbounded. Hence, there exist {/ifcjfc^i C ^ and {r„}5^]^ C (0, r) such that 

hm T„ = and lim /ia:(^2 - t„) = /io(^2 - t„) for n = 1, 2, • ■ ■ . (2.27) 

n—i-oo fc->-oo 

For each k > 1, set 

Qfc(V5,V') = C* +5(V' - V'2)(/ifc(^) -V'), ^oW - f^<^ < hkilp), '4)2 - T < Ip < 1p2, 

where 5, which is independent of /c, is a positive number so small that 

g(/io(V') -Z^,-^) <9fc(^o(V') -Ai,V'), V2-'r<^<^2 (2.28) 

and 

— (99,V'2 -r) > — (99,?/'2 -r), ho{ip2 - t) - fi < ip < hk{4'2 - r). (2.29) 



Note that 
dip'^ 



(yj, V) = <^(2/lUV') + (V' - V'2)/il'(V')) < 0, /lo(V') - /U < (^ < /ifc(V), iP2-r<i)<ij2 



due to the definition of ^. Thus g^ is a supersolution to (j2.3p . Then one can get from the classical 
comparison principle, together with (|2.28|) and (|2.29|) . that 

q{(p,ip) < c* + 6{ip - tl)2){hk{'>P) -^), ho{ilj) - fi < (p < hkit/j), V'2 - ^ < 1/^ < 1/^2 , ^ = 1,2, ••• , 
which, together with ()2.27p . yields 

dq 

■7:-{ho{ll^2-Tn),ll)2-Tn)>dTn, n = 1, 2, • • • . (2.30) 

dip 
Then (|230|) . l^iM), (12:251) and -oo < /io(V'2) < imply that 



Additionally, ((OD yields 

B'ic,)^ihoiiP2),^2) = -A"(c.)(|^(/io(V'2),V'2))' - i?"(c,)(^(/io(V'2),V'2))' > 0. (2.32) 
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It follows from I^Mil, (12:32]) and -oo < /io(V'2) < that 



dq 



_ d'^q 



do 

(/lo('02),V'2)/lo('02) + ^(/io(V'2),^2) G (0,+00]. 



(2.33) 



Due to (I2.33P and ttt (ho (ip2),''p2) = 0, one can get that for sufficiently small e S (0, r) 

OTp 

dq 

— {ho{ij),ip)<0, ^2 - e < V' < -02, 



which contradicts (j2.26p and q is subsonic on the left of 5*0 . □ 

According to Lemmas 12. 4112.51 and Proposition [231 one can get the structure of exceptional points. 

Remark 2.4 Consider a C^ transonic flow in the potential plane. If an interior point in the nozzle 
is exceptional, then near the exceptional point, the sonic curve can he represented as (p = h{ip) which 
achieves a local maximum (if the left of this point is subsonic) or minimum (if the right of this point 
is subsonic) at this point. 

Remark 2.5 Under the assumptions of Lemma \2.5[ one further gets S[{0) = in (i) and S[{1) = 
in (ii). 



exceptional point 




subsonic ♦ supersonic 



y^ 



Figure of Remark [27 



^ 



subsonic 



(gl(l).g2(l)) 



subsonic 



supersonic 

exceptional point 
supersonic 




(Si (0), 52(0)) ^ 
Figure of Remark 12.51 



Remark 2.6 We never use supersonic regions in the discussion in § 2.2, 2.3, 2.4 except the orientation 
of sonic curves. So, the results in these three subsections (except Lemma \2.1\) also hold for any C^ 
subsonic- sonic flow or subsonic-sonic-subsonic flow whose sonic points is a curve. 

2.5 Structure of sonic curves 

Consider a nozzle in the potential plane of the form 

G = {{^,^) ■■ gi{ip) <^ < 52(1/^), -01 < -0 < -02}, 

where ipi < 'ip2 and 

9i{4')<92W, ipi < ip < '>p2- 

The structure of the sonic curve of a C^ transonic flow in G is as follows. 
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V'l 




subsonic 




subsonic 



supersonic 



Figure of Theorem 12.11 




Figure of Theorem [27 



Theorem 2.1 Let S be the sonic curve of a C^{G) transonic flow of Meyer type. Then, S is a 
disjoint union of three connected parts (may be empty): Sg, 5"+, S^, where Se is the set of exceptional 
points, while 5+ and S- denote the other two connected parts approaching the upper and lower walls, 
respectively. Furthernnore, 

(i) Se is a closed segment parallel to the ip-axis; 

(ii) 5+ and S- are two graphes of function with respect to ip, respectively. Along the positive 
direction of S, ip is strictly decreasing on 5+ while strictly increasing on 5_ . 

(iii) // Se is empty, then S = S+ or S = S-. 

Proof. Denote 

S:^ = Si{s),i; = S2{s),{S[{s)f + {S'^{s))^>0, < s < 1, Si, S2 e C\[0,1]). 

Without loss of generality, it is assumed that ^2(0) = ipi and ^2(1) = ^2, which is equivalent to that 
the subsonic region is located on the left of S. Note 

^(Siis), S2{s))S[{s) + ^iSiis), S2{s))S',is) = 0, < s < 1. 



It is clear that 



if S[ (s) = for some s G [0, 1], then ^(SUs), S2(s)) = 0. 

dip 



Lemma 12.41 and Proposition 12.31 show that 

dq 



if -jiSiis), 52(s)) = for some s G (0, 1), then S[{s) = 0. 



(2.34) 



(2.35) 



Owing to (|2.34|) and Proposition 12.21 5^(0) 7^ if (51(0), ^2(0)) is not exceptional. There are four 
cases to be considered. 

(i) The case that (5'i(0), 52(0)) is exceptional and ^{(O) > 0. Then, it follows from Propositions 
[221 O and Lemma [23] that S[{0) = 0. Set 

s* =sup{0 < s < 1 : S[{s) = 0}. 

Due to ()2.34p and Lemmas 12.61 12. 7^ one gets that 

S[{s) = 0, 0<s<s*. (2.36) 
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Indeed, if (|2.36|) is false, there exist si,S2 € [0, s*] such that si < S2, S'i{si) = 5'[(s2) = and 

S[{s)^0, si<s<S2. (2.37) 

It follows from (j2.34p and Lemma 12.71 that 

Si{s) < S'i(si), si<s<si+e, (2.38) 

Si{s) < Si{s2), S2-e<s<S2 (2.39) 

for some positive number e (< S2 — si). Additionally, Lemma 12.61 yields 

Si{s)>mm{Si{si),Si{s2)}, si < s < S2. (2.40) 

Due to ()2.38p - p.40p . Sj = either on [si, si +e] or on [s2 — e, S2], which contradicts p.37p . Owing to 
(lOGjl . (i23i]l and Proposition O (5i(s), 52(s)) is exceptional for each s e [0,s*]. If s* = 1, 5^ = 5 
is a closed segment parallel to the ^-axis. Otherwise, < s* < 1. Lemma 12.71 and the definition of s* 
imply 

S[{s) < 0, s* < s < 1. 

It follows from Proposition 12.21 and ()2.35p that {Si{s),S2{s)) is not exceptional for each s E (s*,l). 
Furthermore, Propositions 12.2 1 12.31 and Lemma 12.51 show that (5i(l), 52(1)) is not exceptional either. 
That is to say, Se = {{Si{s), S2{s)) : < s < s*} and S+ = {{Si{s),S2{s)) : s* < s < 1}. 

(ii) The case that (51(0), ^2(0)) is exceptional and 5*5^(0) < 0. Let us prove that S[{s) < for 
each s G [0, 1] by contradiction. Otherwise, there exists sq G (0, 1] such that S[{so) = and S[{s) < 
for each s £ [0, sq). Then, ()2.34p and Lemma [2771 show that Si{so — e) < Si{so) for sufficiently small 
e € (0,so), which contradicts that S[{s) < for each s € [0,so)- Therefore, 

S[{s) < 0, < s < 1. 

It follows from Proposition 12.21 and (|2.35p that {Si{s),S2{s)) is not exceptional for each s G (0,1). 
Furthermore, (S'i(l), ^2(1)) is not exceptional either due to Propositions 12.21 12.31 and Lemma [2.51 
That is to say, Se = {(5i(0), 52(0))} and S+ = {(5i(s), ^2(5)) : < s < 1}. 
(iii) The case that (51(0), £'2(0)) is not exceptional and 5^(0) > 0. Set 

s, = sup{0 < s < 1 : S[{s) > 0}. 

Then < s^, < 1. Let us prove 

S[{s) > 0, < s < s* (2.41) 

by contradiction. Otherwise, there exist sq G (0, s*) and eq G (0,s* — sq] such that S[{so) = and 
S[{s) > for each s G (so,so + eo)- Then, (|2.34p and Lemma [2.71 show that 5i(so + e) < 5i(so) 
for sufficiently small e G (0, £9)1 which contradicts that S'i{s) > for each s G (sojSq + eo)- Thus 
(5i(s), ^2(5)) is not exceptional for each s G (0,s*) due to (12.4ip . (12.35P and Proposition 12.21 If 
s* = 1, then either 5_ = 5, or 5_ = {(5i(s), ^2(5)) : < s < 1} and Se = {(5i(l), 52(1))}. 
Otherwise, < s* < 1 and 5((s*) = 0. Set 

s* = sup{s* < s < 1 : S[{s) = 0}. 
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Then, < s* < s* < 1. Similar to the discussion in (i), one can get that 

S;(s) = 0, s^<s<s*; 

furthermore, 5_ = {{Si{s), S2{s)) : < s < sj and Se = {(5i(l), 52(1)) : s* < s < 1} if s* = 1, while 
5_ = {{Si{s),S2{s)) : < s < sj, S^ = {(5i(l), 52(1)) : s, < s < s*} and 5+ = {(5i(s), 52(s)) : 
s* < s< 1} if s* < 1. 

(iv) The case that (5i(0), 52(0)) is not exceptional and 5^(0) < 0. As the proof of (ii), one can 
prove that S[{s) < and (5i(s), 52 (s)) is not exceptional for each s € [0, 1]. Then S-^- = S. D 

Theorem 2.2 For any C'^{G) transonic flow of Taylor type, each sonic point is not exceptional and 
if is strictly monotone along the sonic curve. 

Proof. Without loss of generality, we assume that 5 is the sonic curve intersecting the lower 
wall. Denote 

5:(^ = 5i(s), ^ = 52(s), (5Us))' + (5^(s))2>0, < s < 1, 5i,52 € C2([0, 1]) 

with 52(0) = 52(1) = ^2 and 5i(0) < 5i(l). Lemma EH implies that 5^(0) > and 5^(1) < 0. 

First we show that there is no exceptional point in the nozzle. Otherwise, it is assumed that 
(5i(so), 52(so)) is an exceptional point for some sq € (0,1). Owing to Remark 12.41 52(so) > and 
5i(so) achieves a local maximum if the left of (5i(so), 52(so)) is subsonic, while 52(so) < and 
5i(so) achieves a local minimum if the right of (5i(so)i 52(so)) is subsonic. Note that 52(0) = 52(1) 
and 5i(0) < 5i(l). If the left of (5i(so), 52(so)) is subsonic, then Lemma [2T71 vields that there exists 
si G (so, 1) such that 52(si) > 0, 5i(si) achieves a local strict minimum and the left of (5i(si), 52(si)) 
is subsonic. Similarly, if the right of (5i(so), 52(so)) is subsonic, then there exists S2 G (0,so) such 
that 52(^2) < 0, 5i(s2) achieves a local strict maximum and the right of (5i(s2), 52(52)) is subsonic. 
Each case contradicts Lemma [2.6[ Therefore, (5i(s), 52(s)) is not exceptional for each s € (0,1). 
Furthermore, Proposition 12.21 implies that 5i is a strictly increasing function on [0, 1]. 

subsonic subsonic 

, . N _ \subsonic subsonic / _ V. . 
subsonic 1^ ~ ^^ \ / ^ ~ ^^\Ji^^som.c 
^ supersonic\ /supersonic X..^ 

the left of {Si{so),S2{so)) is subsonic the right of (5i(so), 52(so)) is subsonic 

Let us show that (5i(0), 52(0)) is not exceptional by contradiction. Otherwise, 5j^(0) > and 
Lemma [2 . 5 1 vield that 5 near (5i(0), 52(0)) is a graph of function with respect to if). Then, Propositions 
12.21 12.31 implv that 5i(s) < 5i(0) for sufficiently small s > 0, which contradicts that 5i is strictly 
increasing on [0, 1]. Similarly, one can prove that (5i(l), 52(1)) is not exceptional either. D 

Theorems 12.11 and 12.21 can be described in the physical plane. Assume that the nozzle in the 
physical plane is of the following general form 

Q = {{x,y) ■■ fi{x) <y < f2{x)M <x < h], 
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where /i < I2 and 

fiix) < f2{x), h<X<l2. 

The counter part of Theorem 12.11 is the fohowing 

Theorem 2.3 For any C^ transonic flow of Meyer type in the nozzle Q, its sonic curve S is a disjoint 
union of three connected parts (may he empty): Se, 5+, 5_, where Se is the set of exceptional points, 
which is a closed line segment, while S^ and S- denote the other two connected parts approaching the 
upper and lower walls, respectively. Moreover, if Se is empty, then S = S+ or S = S-. 

Similarly, in the physical plan, Theorem 12.21 becomes 

Theorem 2.4 For any C^ transonic flow of Taylor type in the nozzle Q , each point on its sonic curve 
is not exceptional. 





subsonic 
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2.6 Characteristics from sonic points 

Characteristics play an important role in understanding the behavior of supersonic flows near sonic 
curves. So we first study characteristics from sonic points. 

As will been shown in § 3.3, in the supersonic region, it holds that 



dW dW 

— ((/?, V') + /3(9?, ^)^(¥', V') + ^{V, i^)W{^, ^P) = 0, 

dZ dZ 



(2.42) 
(2.43) 



where 

Ld{ip,1p) 



-^'(^^^'^^)^^^^^-^^±^.^(.,^)-l)^^7l-^.^(.,^) 



B'{q{^,i^)) 



-l/(7-l)-l/2 



:l±ig^(^, V^) m^q\,^, V) - 1) ^' (l - ^q\^, ^) 



2 



1/(7-1) da 



-(-*)-'(*-«)|(-«-1f^l^(-^)- 
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z(.,,)^-.',„.,,))|(..,)-^;M|(.,«^ 

Note that /3 is positive in the supersonic region and vanishes at the sonic curve. The positive and 
negative characteristics of the system (I2.42P , (I2.43P are governed by 

^'+(V') = /3KW>V') and <I>'„(V') = -/3(<I>'„W,V), 

respectively. 

The following a priori estimates hold along characteristics. 

Lemma 2.8 Assume that S+ (S_) is a positive (negative) characteristic. Then, 9{^±{tp),'ip) =F 
H{q{^±{ip),ip)) is invariant on T,±, where 



H{q) = r ( - ^{s)B'{s)) ^'^ds, c,<q<q„ 



Proof. On S±, it holds that 

r/ f)f) BO 

-^'(g($±(V),^))^($±(V),V') 

= ±^if(g($±(V),V)). 

So, 6'(<l>±('(/'),'0) =F-f^(9(*^±(V')>V')) is invariant on S±. D 

Consider a C^ transonic flow of Meyer type in Q whose sonic curve S intersects the upper wall at 
(xi, /i(xi)). Assume that the subsonic region is located on the left of 5, 5 = S+USe, 5+ is not empty 
and Se is a line segment with nonempty interior. In the potential plane, the lower wall corresponds 
to ^ = and the sonic curve is given by 

S:^ = Si{s),i^ = S2{s),{S[{s)f + {S'2{s)f>0, < s < 1, 81,82 e C\[0,1]) 

with 52(0) = and 82(1) = m. Set 

s* = sup{s : 8'i{s) =0,0 < s < 1}. 

Then, as proved in Theorem 12.11 < s* < 1, 

S[is) = 0, |^(5i(s),52(s)) = 0, 0<s<s*, (2.44) 

S[{s)<0, ^{8i{s),82{s))>0, s*<s<l. (2.45) 

Let us first study characteristics from (5i(s), 5*2(5)) for < s < s*. 
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Lemma 2.9 Assume thatO < s < s* . Then the positive and negative characteristics from, {Si{s), S2{s)) 
coincide and are given by 

$+(V') = ^-W = Si{s), 0<i;< S2{s*). (2.46) 

Proof. It is not hard to check that the functions given by (|2.46p are positive and negative char- 
acteristics from (5i(s), 82(8)), respectively. Let us verify the uniqueness of the positive characteristic. 
Assume that $ is a positive characteristic from (Si{s),S2{s)). Since $ is a nondecreasing function, it 
suffices to verify that $ coincide with the functions given by (J2.46P in a right neighborhood of 5*2(5). 
Owing to (123]) and dTIi]) . one gets that 

^"(c*)(^(5i(5), V))' + B'ic,)^iSi{sU) + B"(c,)(^(5i(s),V^))' = 0, 0<^< S2is*) 
and 

-{S^{sU)=0, - 
which imply 



^(5i(5),V) = 0, ^(S,{s),i;)=0, 0<V<52(5*), 



|^(Si(s),V)=0, 0<^<S2(s*). (2.47) 

Fix a small positive number 5 such that (S'i(s), Si{s) + 5]x [52(s), ^2(5) + 5] belongs to the supersonic 
region. Due to (|2.47p . there exists M > such that 

0<;5((^,V) <M((^-5i(s)), Siis) < ^ < Siis) + 5, S2is) < i^ < 82(3) + 5. (2.48) 

Choose a number e G (0, 5) such that l>(52(s) + e) < Si{s) + 6. Owing to (12.481) . I> satisfies 

0<l>'(^)<M($(V^)-5i(s)), S2{s) < i^ < 82(8) + e; ^S2{s)) = Si{s). (2.49) 

The unique solution to ()2.49p is 

^i;) = Si{s), S2{s)<iP<S2{s)+e. 

The uniqueness of the negative characteristic can be proved similarly. D 

We now turn to characteristics from (S'i(s), S2{s)) for s* < s < 1. 

Lemma 2.10 Assume that s* < s < 1. Then, there exist uniquely a positive and a negative charac- 
teristics from (5i(s), 52(5)), which contain no other sonic point except (5i(s), 52(s)). 

Proof. According to (|2.45p . the sonic curve from {Si{s*),S2{s*)) to (S'i(l), 52(1)) is a graph of 
function with respect to (p, which is denoted by 

iP = S{ip), Si{l) < ^ < Si{s*) 

with 

S'iv) = -^{^,S{^))[^{ip,S{^))y' e (-00, +00), 5i(l) <^< Si{s*). (2.50) 
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First we show that {Si{s), 82(3)) is the unique intersecting point between any characteristic from 
(S'i(s), <S'2(s)) and the sonic curve. Assume that $ is a characteristic from {Si{s),S2{s)), which is 
defined in a right neighborhood of S2{s). Owing to ()2.50p and $'(52(s)) = 0, the path of $ leaves for 
the supersonic region after ip = S2{s) locally. Therefore, $ is strictly monotone near ip = 5*2 (s) and 
its inverse function is denoted by ^, which solves 

if $ is a positive characteristic, while 

if $ is a negative characteristic. Since ^' = 00 at sonic points, it follows from (j2.50p that the path of 
^ never approaches {{Si{s),S2{s)) : s* < s < 1} after it leaves {Si{s), 82(3)). Let sq € [0, {s + s*)/2] 

satisfy 

{s + s* "1 
82(3) : < s < — - — >. 

Then, sq > s* since S £ C'^ and S^ is a segment parallel to the -i/^-axis. Owing to the monotonicity 
of ^, the path of \1' never approaches {{Si{s), 82(3)) : < s < sq} if <^ is a positive characteristic, 
while never approaches {(5i(s), 5*2(5)) : < s < s} if $ is a negative characteristic. Summing up, 
(jS'i(s), 5*2 (s)) is the unique intersecting point between the path of ^ and the sonic curve. 

Below we prove the local existence and the uniqueness of negative characteristics under the addi- 
tional assumption 8'{8i{s)) < 0, and the other cases can be discussed similarly. Since q{Si{s), 5*2(5)) = 
c^, there exists a positive number tq such that 

Go = {(v3,V) ■■8iis)-To<^<8i{s),8{ip) <ij<8iip) + To} 

belongs to the supersonic region and 

^ >1- min 8', ((^,V)eGo. (2.51) 



I3{ip,l/j) [Si{s)-To,Si{s)] 

For < e < tq/2, q{8i{s), 82(3) + e) > c*. Therefore, there exists uniquely a solution ^^ to the 
following problem 

K{^) = - .. I , ,, , ^<8i{s)- ^,{8i{s)) = 82{s) + e. (2.52) 

Let [ife, 5*1(5)] be the maximum interval of existence for the solution to the problem ()2.52p on Go- It 
follows from (j2.5ip that 

S{^) + {8i{s) - if) < 82{s) + e + (l - ^ min 5') (5i(s) -(/.)< *,((^), ^e<V<8i{s), 

(2.53) 
which shows that the path of ^^ does not intersect OGq n 8. Therefore, for < ei < £2 < to/2, 

V?£i < V'ea < 9^x0/2 and ^^.(v?) < ^£2(^3) for each 93^^/2 <(/?< 5*1(5). (2.54) 
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Let 

^{if) = lim ^e{^), ipro/2 <V<Si{s). 
e— >-0+ 

Then, it follows from (j2.52p - (|2.54p that ^ solves the following problem 

^V) = -^^^4py' "^-0/2 < V < Siis); ^(S.is)) = S2{s) (2.55) 

and satisfies 

S{ip) + {Sl{s)-^)<^{^)<^ro/2{'p), Vro/2 < ^ < Sl{s). (2.56) 

Note that (l235l) and (12361) imply 

^'{^) < 0, ^ro/2 <V< Si{s). 

Thus ^ is strictly decreasing and its inverse function is a negative characteristic from (5'i(s), 5*2(5)) 
due to ([235]) . 

Let us turn to the uniqueness. Assume that $ and ^ are two negative characteristics from 

dq 

(5'i(s), 5*2(5)), which are defined in a right neighborhood of 52(s). Since — — (5i(s), 52(s)) > from 
oyj 

(j2.45p . there exists a positive number ri (< tq) such that 

^|^(5i(.),52(s)) < |i((^,^) < ^^(5i(.),52(s)), (^,V) G Gi, (2.57) 

where 

Gi = {(99,^) : Si{s) -n<ip< 5i(s),5(^) < V < S{ip) + n}. 

Since $'(52(s)) = $'(52(s)) = 0, there exists a positive number (5 so small that 

{^{^P),^P),{^{^),^P)eGl, S2{s) < ip < S2is) + 6 (2.58) 

and 

dq 



< sup 

Gi 

< sup 

Gi 



9g 



(?(/? 



(5i(s) - '^m < l^iSiis), S2{s)){i; - S2{s)), S2{s) < ^ < 52(s) + 6, (2.59) 
(5i(s) - ^{^P)) < l^iSiis), 52(s))(^ - 52(s)), 52(s) < ^ < 52(s) + 5. (2.60) 



Set 

/i(V) = ($(V) - '^(V'))', ^2(5) < V < S2{s) + 5. 

Then, /i solves 

/i'(V')=C(V')('?(C(^),V')-c.)-^/'/i(V), 52(s)<V<52(s) + (^; /i(52(s)) = 0, (2.61) 

where ^ G C([52(s), 52(s) + 5]), while 

min{$(^), $(V')} < C(V') < max{$(^), $(V')}, 52(s) < ^ < 52(s) + 6. (2.62) 
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It follows from (12371) - ([2:60]) and (|232|) that 

(?(CW,V)-c*<sup 1^ (5i(.)-CW) + ||^(5i(s),52(s))(^-52(s)) 

^7l7('5i(^)' '52(s))(V' - S2(s)), S2(S) < ^ < S2is) + 5, 
4 aV' 



(2.63) 



g(C(V), V') - c* > - sup ^ (Siis) - cm + ]:Pr{Siis), S2{s)){i; - S2{s)) 
>^||(5i(5), 52(s))(V' - S2{S)), S2{S) <^< S2{S) + 5. 



(2.64) 



Substituting ([TM]) and (l23il) into (f23T]) yields 

|/i'('0)| <M(V'-52(s))"i/2;^(V'), S2{s) < ^ < S2{s) + 5- /i(52(s)) = (2.65) 

for some positive constant M. Since (j2.65p admits only a trivial solution, the uniqueness is proved. D 



For fixed s G (s*, 1), according to Lemma [2. 101 there exist uniquely a positive and a negative char- 
acteristics from (5i(s), 5*2(5)); furthermore, for each characteristic, the lower endpoint is {Si{s), S2{s)), 
the upper endpoint is located at the upper wall or the outlet of the nozzle. Without loss of generality, 
the nozzle is assumed to be so long that all upper endpoints of the characteristics from (S'i(s), 5*2(5)) 
for s G (5*,!) are located at the upper wall. Therefore, the positive and the negative characteristics 
from (5i(s), 52(5)) are governed by 

^ ,,, ;^<I>.,H_(V)=/3(^s,+ (V),V'), 52(5)<^<7n, 



and 



^.,+(^2(5)) = 5i(5), $,,+(m) = ^s,+ 
'^-$,,_(^) = -/3(*,,_(V'),V'), ^2(5) <iJ<m, 



$,,_(52(5)) = 5i(s),$,,_(m) = (^,,_ 

respectively, where s G (s*, !)■ 

For s* < si < S2 < 1, the uniqueness of the characteristic shows 



and 



Set 



fsu+ > V'S2,+ , fsu^ > ^ 



S2, 



'^su+W > '^S2,+W^ '^s^,-W > ^s2,-('0), ^2(52) <'il^<m. 



(/?_ = lim ifs-, ip+ = lim ips^+. 



(2.66) 
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Remark 2.7 It is not clear at this moment whether the positive and negative characteristics from 
(5i(s*), 5*2(5* )) are unique. However, 



Si{s*), 



< ^ < 52(s*), 



lim ^s +W, S2{s*) Kip <m 



and 



Siis*), 



< V < S2is*), 



lim ^s -{tp), S2{s*) <il) <m 



are the maximal positive and the minimal negative characteristics from (S'i(s*), 5*2(5*)), respectively. 

m 




subsonic 



(5lf5*),52(5*)) 

supersonic 



Figure: characteristics for the transonic flow in the potential plane 

Let X* be the point such that (x*,/i(x*)) is transformed to {(p-,m) in the coordinates transfor- 
mation from the physical plane to the potential plane. That is to say, 

X* = sup {x € [^1, ^2] '■ {x, fi{x)) is the endpoint of the negative characteric from a point on 5+}. 

(2.67) 

Proposition 2.4 For each x G [a;i,x*], /f(x) > 0. 

Proof. Fix s G (5*, 1). It follows from the definition of W that 

1 dq 



A'(.(^,*,,(^)))=a^(^'^--^(^)) 

dq 



(3{^,-i>s,+ {^))di; 



{(p, ^s, + {(p), Si{s) <(p< ips,+ . 



Since /3(5i(s), 52(s)) = and — — (5i(s), S2{s)) > 0, there exists a positive number 5 (< (ps,+ — 5i(s)) 
such that W{Si{s) + 6, ^,,+ (5i(s) + 6)) > 0. Then, ^Mh yields 

Wi(ps,+,m)>0. (2.68) 

Similarly, one can get 

Z{ips-,m)>0. (2.69) 

Due to to the arbitrariness of s G (s*, 1), it follows from (j2.68p and (j2.69p that 

W{(p,m)>0, Si{l) < if < (p+, W{(p+,m)>0 

and 

Z{ip,m) > 0, 5i(l) < if < (^_, Z(c^_,?n) > 0. 
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Thus, 



9-0^^' ' 2A'{q{^,m)) 



{W{^,m) + Z{ip,m))>Q, 5i(l) < (/? < 99^ 



Moreover, ()2.45p gives 



dq 

dip 



(5i(l),m)>0. 



(2.70) 



(2.71) 



Assume that (x,/i(x)) is transformed into {^{x),m) for li < x < I2 in the coordinates transformation 
from the physical plane to the potential plane. Then, 



^'ix)= ^^,"^/^f,^\, =g(rr,/i(x))Jl + (/((x))2>0, h<x<k 
coso[x, Jl[X)) V 



and 



n^)B'iqmx),m))^mx),n.) = ^^^ = TT(^' ^^ < ^ < '^ 
The conclusion of the proposition follows from (|2.70p - ()2.73p . 



(2.72) 

(2.73) 
D 





Figure of Proposition 12.41 
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Remark 2.8 Proposition \2.4\ still holds if Se is empty or a single point set. 



Similar to Lemmas 12.91 12.101 Proposition 12.41 and Remarks 12. 7| 12. 8| one can prove the following 
two general conclusions. One is on characteristics from sonic points, and the other is on the geometry 
of walls. 





Figure of Theorem [27 
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Theorem 2.5 Let S = 5+ U 5e U S- be the sonic curve of a C^ transonic flow of Meyer type in Q. 
(i) There exist uniquely a positive and a negative characteristics from each interior point of 5+ U 

(ii) All positive and negative characteristics from points of Se \ d{Sj^ U S-) are the same line Se- 
(iii) There are positive and negative characteristics from points of Se H 9(5+ U5_). Particularly, 
the maximal positive and the minimal negative characteristics are unique. 



Characteristics from boundary sonic points are as follows. 



Remark 2.9 Let S = 5+ U 5e U S^ he the sonic curve of a C^ transonic flow of Meyer type in Q. 
Assume that the subsonic region is located on the left. For the sonic point on the upper wall, 

(i) If it belongs to 5-|_, then there is not any characteristic from this point. 

(ii) If it belongs to Se \ dS- , then the positive and negative characteristics from this point are the 
same line Se- 

(iii) // it belongs to Se H dS^, then there are positive and negative characteristics from this point. 
Particularly, the maximal positive and the minimal negative characteristics are unique. 

(iv) // it belongs to S-, then there exist uniquely a positive and a negative characteristics from 
this point. 

For the sonic point on the lower wall, 

(i) If it belongs to S- , then there is not any characteristic from this point. 

(ii) If it belongs to Se \ dSj^ , then the positive and negative characteristics from this point are the 
same line Se- 

(iii) // it belongs to Se H 05+, then there are positive and negative characteristics from this point. 
Particularly, the maximal positive and the minimal negative characteristics are unique. 

(iv) // it belongs to 5+, then there exist uniquely a positive and a negative characteristics from 
this point. 

Theorem 2.6 Let S = 5+ U 5e U S- be the sonic curve of a C^ transonic flow of Meyer type in Q . 
Assume that the subsonic region is located on the left and S intersects the upper and lower walls at 
(xi,/i(xi)) and (x2, 72(3^2)); respectively. 

(i) If S+ y^ 0, then /('(x) > for each x G [xi,x*], where x* is given by ()2.67p . 

(ii) If S^ 7^ 0, then f2{x) < for each x G [x2,x*], where 

X* = sup |x G [/i, ^2] : {x, fiix)) is the endpoint of the positive characteric from a point on 5_}. 

(2.74) 

As a consequence of Theorem 12.41 and the proof of Proposition 12. 4^ we can obtain 

Theorem 2.7 Consider a C^ transonic flow of Taylor type in Q . Assume that the sonic curve in- 
tersects the upper wall at (xi^_, /i(xi^_)) and (xi^+, /i(xi^+)), while the lower wall at {x2-,f2{x2-)) 
and (x2,+ , /2(x2,+)) with Xk,~ < Xk,+ {k = 1,2). Then, 

fi{x) > 0, xi__ < X < xi^+ 

and 

/2'(x) < 0, X2,- < X < X2,+ . 
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Figure of Theorem 12.71 



Remark 2.10 For a smooth transonic flow of Taylor type, it has been mentioned in ^ that the bound- 
ary of a supersonic enclosure on a wall cannot contain a straight segment. Theorem \2. 7| strengthen it 
to be that the curvature of this boundary must be nonzero. 



2.7 Instability of transonic flows with nonexceptional points 

According to Theorems 12.11 and 12. 6| we are ready to show that C^ transonic flows with nonexceptional 
points are unstable with respect to small changes in the shape of the nozzle. 

/('(xi) = y\ y = fi{x) 




f^{x,) = o y = M-y 



Figure of Theorem [27 

Theorem 2.8 Consider a C^iG) transonic flow of Meyer type, whose subsonic region is located on the 
left. Assume that the sonic curve S = 5+U5eU5_ intersects the upper and lower walls at (xi, /i(xi)) 
and {x2, f2{x2)), respectively. If there is a nonexceptional point on S, then (p is unstable with respect 
to small changes in the shape of the nozzle in the following sense. 

(i) The case S+ ^ 0. Let x* be given by (j2.67p . If fi S C°°{[li,l2]) is a C^ small perturbation 
of /i and satisfies /('(xi) = with some xi < fi < x* , then there is not a C'^{Q) transonic flow of 
Meyer type with Q being the corresponding domain when the upper wall is replaced by /i, which is a 
C^ small perturbation of the background flow. 

(ii) The case S^ ^ 0. Let x* be given by (|2.74|) . If f2 € C°°{[li,l2]) is a C^ small perturbation 
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of f2 CLnd satisfies /2(^2) = with some X2 < X2 < 2;*, then there is not a C'^{G) transonic flow of 
Meyer type with Q being the corresponding domain when the lower wall is replaced by f2, which is a 
C^ small perturbation of the background flow. 

Proof. We prove (i) only and the proof of (ii) is similar. Consider a C'^{Q) transonic flow of Meyer 
type, whose sonic curve is denoted by 5 = 5+ U 5e U S- according to Theorem 12.11 If the transonic 
flow is close enough to the background flow in C^'-norm, then 5+ 7^ and its two endpoints are small 
perturbations of the two endpoints of 5+, respectively. Thus, xi and x* are small perturbations of xi 
and X* , respectively, and it follows from Theorem 12.61 that 

fi{x) > 0, xi < X < X*, 

where xi and x* are defined in a similar way as for xi and x*, respectively. This contradicts that 
/f (5i) = 0. D 

Remark 2.11 The instability in Theorem \2.8\ is weak. It is unknown whether the flow is unstable if 
the nozzle wall is perturbed in C^ or other smooth spaces. 

Remark 2.12 Theorem \2.8\ motives us to seek a smooth transonic flow of Meyer type whose sonic 
points are exceptional. Indeed, we do prove the existence of such a smooth transonic flow for the 
symmetric de Laval nozzle Q in the paper. 

Moreover, for C^ transonic flows of Meyer type in Q,, one can prove 

Proposition 2.5 Assume that /"(O) = 0. Consider a C^ transonic flow of Meyer type in Q, whose 
each streamline in the supersonic region is a graph of function with respect to x. If there is a nonex- 
ceptional point at the sonic curve, then not all streamlines are convex. 

Proof. It is assumed that the subsonic region is located on the left without loss of generality. 
We transform the flow from the physical plane to the potential plane. It is assumed that (0, /(O)) is 
transformed to (0, m) and the sonic curve is transformed to 

S:^ = Si{s), ^ = 52(s), {S[{s)f + {S'^{s)f > 0, < ., < 1 

with 52(0) = and 5'2(1) = m. Set 

s* = sup{s : S[{s) =0,0 < s < 1}. 

The symmetry of the nozzle and Proposition 12.21 show that (51(0), 5*2(0)) is an exceptional point and 
0(S'i(O), 5*2(0)) = 0. Similar to the discussion in §2.5, one can prove that < s* < 1 and 

5;(s) = 0, ^(5i(s),52(s)) = 0, e(5i(s),52(s)) = 0, < s < s*, 

S[{s)<^, ^{Si{s),S2{s))>0, 9{Si{s),S2{s))<0, s* < s < 1. 

Lemma 12.81 implies that 

e{ip-,m) + H{q{ip^,m)) = 0, e{^+,m) - H{q{ip+,m)) = 0, 
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where (^± are defined by (j2.66p and H is given in Lemma 12.81 Thanks to q{^±,m) > c*, one gets that 
0{f-, m) < < 9{ip-^,m) and thus (/?_ < < c^+. As shown in Proposition [231 one can get W{0, m) > 
and thus Z{0,m) = —W{0,m) < 0. Denote S_ to be the negative characteristic from (0,m). Due 
to Lemma [2.81 G~\~H{q) is invariant on S_, which equals identically to 9{Q, m) + H{q{0, m)) > 0. Note 
that 

e{S,{s),S2{s))+H{q{S,{s),S2(,s))) = 9{Si{s),S2{s)) < 0, < s < L 

Thus S_ never approach the sonic curve. Assume that S_ intersects the lower wall at (ipi,0). Then 
5i(0) < ifi. Since the flow is sonic at (5i(0), 0) and supersonic at ((/9i, 0), there exists ip2 G (5i(0), c^^i) 
such that W{ip2,0) < 0. Let 5]_|_ be the positive characteristic from ((/?2,0) and denote {if*,ip*) to 
be the intersecting point of T,^ and S_. It follows from (|2.42p and ()2.43p that W{ip^,'tlj^:) < and 
Z{ip^,,il)*) < 0, which imply 



dq 

dip 



{v^*,ip* 



"2A'(g((^„V'*)) 



{W{lp^,iIj^) + Z{lp^,iIj^)) <0. 



Assume that y = /*(x) is the streamline across (c/^*, ^*) and (x^,, /*(x*)) corresponds to (99^,, ^*) in the 

dq 



coordinates transformation. Then, as shown in (j2.72p and (j2.73p . /^'(x*) and -^(^3*, '(/'*) are of the 
same sign. Thus /^'(x*) < 0. D 




\ supersonic 



Figure of the proof of Proposition 12.5 



Formulation of the smooth transonic flow problem and main re- 
sults 



As mentioned in Remark I2.12[ we seek a smooth transonic flow of Meyer type, whose sonic points 
are exceptional, in the symmetric de Laval nozzle VL. In this section, let us formulate this smooth 
transonic flow problem both in the physical plane and in the potential plane, and state the existence 
and uniqueness theorems. 



3.1 Formulation of the transonic flow problem in the de Laval nozzle 

To seek a smooth transonic flow in Vt whose sonic points are exceptional, the first step is to determine 
the location of the sonic curve. It is not hard to verify from Lemma 12.11 and Proposition 12.11 that 
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Proposition 3.1 For a C'^{Q) transonic flow of Meyer type, the following statements are equivalent 
(i) The endpoint of the sonic curve on the wall is located at the throat point; 
(ii) The flow angle always equals to zero at the sonic curve; 
(iii) Every point at the sonic curve is exceptional; 

(iv) The velocity vector is along the normal direction at the sonic curve; 
(v) The potential is a constant at the sonic curve; 
(vi) The sonic curve is located at the throat of the nozzle. 

Therefore, the sonic curve of the sought smooth transonic flow should be located at the throat. 
That is to say, the flow should be subsonic in the convergent part, then sonic at the throat and then 
supersonic in the divergent part. For such a smooth transonic flow, the velocity vector is along the 
normal direction at the sonic curve and thus the mass flux is 

m=/(o)c:+'/(^-^). (3.1) 

Furthermore, the potential of the flow is a constant at the sonic curve and we normalize it to be zero in 
the paper. Then, the potential of the flow at the inlet, which is a constant, is free. So is the potential 
at the outlet. 

Thus, the problem of a transonic flow in Jl, whose velocity vector is along the normal direction at 
the inlet and the outlet, which satisfies the slip condition on the wall and whose sonic curve is located 
at the throat, can be formulated as 

dw{p{\Vip\^)Vip) = 0, 

'f{g{y),y) = ^11, 

^(x,0) = 0, 
dy 

|^(x,/(x))-/'(x)||(x,/(x)) = 0, 

^{t{y),y) = Cont, 

\Vip{0,y)\=c.,ip{0,y)=0, 
|Vv3(x,y)| < c*, 
|V(^(x,y)| > c^, 

where Cin,Cout ^ ^ are free, the outlet Tout : x = t{y) (0 < y < f{l+)) is free, Q- and 0+ are the 
convergent part and the divergent part of the nozzle, respectively, i.e. 

r2_ = {{x,y) G : X < O}, U^ = {{x,y) gQ : x > O}. 

Since the sonic curve of the smooth transonic flow is located at the throat, we can decompose the 
transonic flow problem (J3.2p - (j3.9p into a subsonic-sonic flow problem and a sonic-supersonic flow 
problem as follows 

div(p(|V(^|2)V(^) = 0, (x,y) G 0_, (3.10) 

v{g{y),y)=Cin, 0<y </_(/_), (3.11) 

^{x,0)=0, g{0)<x<0, (3.12) 

dy 
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{x,y) en, 


(3.2) 


0<y<f{l-), 


(3.3) 


g{0) <x< t(0). 


(3.4) 


I- < x < Z+, 


(3.5) 


< y < /(/+), 


(3.6) 


o<y< /(o), 


(3.7) 


{x,y) gO_, 


(3.8) 


{x,y) £ n+, 


(3.9) 



|^(x,/„(x)) - f'_{x)^ix,f^{x)) = 0, Z_ < X < 0, (3.13) 

\Vip{0,y)\=c„ip{0,y)=0, 
|V(/?(2;,y)| < c^, 



and 



0<y</_(0), 


(3.14) 


(x,y) gO_ 


(3.15) 


(x,y) en+, 


(3.16) 


0<y< /+(0), 


(3.17) 


< X < i(0), 


(3.18) 


< X < /+, 


(3.19) 


< y < /+(/+), 


(3.20) 


{x,y) G 0+, 


(3.21) 



div(p(|V(/.|2)V(/.) = 0, 
\V^{0,y)\=c„ip{0,y)=0, 

§^(x,0)=0, 
|^(x,/+(x))-/;(x)||(x,/+(x)) = 0, 

'p{t{y),y) = Cout, 

|Vc/?(x,y)| > c*, 

where /_ and /+ are the walls of the convergent part and divergent part of the nozzle, respectively, 
i.e. 

/_(x) = /(x), /_<x<0 and /+(x) = /(x), < x < /+. 

Remark 3.1 For the subsonic-sonic flow problem (13.10p - (l3.15p . there are two boundary conditions 
at the sonic curve (see (j3.14p ). Indeed, the potential at the sonic curve is a free constant, which 
is normalized to be zero and thus the potential at the inlet is free (see (13.1ip . where Cin is a free 
parameter). 

It should be noted that to study a subsonic-sonic flow problem and a sonic-supersonic flow prob- 
lem, one must determine and control precisely the speed of the flow near the sonic state, which plays 
a very essential role in the mathematical analysis. In ()3.2p . the speed of the flow is the absolute value 
of the gradient of a solution, which is very hard to estimate precisely. It turns out to be more con- 
venient to solve the subsonic-sonic flow problem ()3.10p - (13.15p and the sonic-supersonic flow problem 
()3.16p - f3.2ip in the potential plane, where the speed of the flow is a solution to ()2.3p . We will thus 
formulate these two problems in the potential plane in next two subsections. 

3.2 Formulation of the subsonic-sonic flow problem in the potential plane 

We will show that there exists a subsonic-sonic flow in the divergent part of the nozzle il_. In this 
subsection, we formulate this subsonic-sonic flow problem in the potential plane. 
Assume that /_ G C^{[1^,0]) (-1 < /_ < 0) satisfies 

6i^^{-x)^- < f'l{x) < 62,^{-x)^- , L<x<0 (3.22) 

with positive constants A_, 5i- and (52,- such that A_ > 2 and (5i,„ < 62-- Choose the inlet Tin as a 
small perturbation of the arc 

ro:x = go{y), 0<y </_(/_), 

where 



,, ; Ul~) [^2 2 n^ ^fn^ T, /-(/-)V(/i(^-))^ + l 

9o[y) = 1-- TTTTT ~ V ^0 - y ' < y < /-(/-), Rq - 
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More precisely, 

Ti,,:x = g{y), < y < f^{l^) 

with g € C3([0, /_(/_)]) satisfying 

<7'(0) = 0, 5(/-(L)) = /_, 5'(/_(Z„)) = -/:(/„) (3.23) 

and 

(3.24) 



where < e(s) < 1 for — 1 < s < and Hm e(s) = 0. 

Since the velocity vector of the flow is along the normal direction at the inlet and the wall of the 
nozzle is solid, the flow angle at the inlet and at the upper wall can be expressed as 

0in(y) = -arctang'(y), < y < /_(/_) 

and 

0_(x) = arctan/i(a;), /_ < x < 0, 

respectively. Let the speed of the flow at the inlet and at the upper wall be denoted by 

qigiy),y) = ^Uy), < y < f_{l^) 

and 

q[x, f{x)) = ^_(x), /_ < X < 0, 

respectively. The incoming mass flux is given by 

™,. = /'-"■' ^^<^^MfM.,. ,3.5) 

Jo COSein(y) 

Since the potential at the sonic curve is normalized to be zero, the potential at the inlet is given by 

C- = f~ -^#kdx. (3.26) 

Jo cosB_(x) 

At the inlet, the stream function is 

V^(5(y),2/) = ^in(y), 0<y </_(/_) 
satisfying *in(0) = 0, *in (/-(/-)) = min and 

cos Bin (y) 



Thus, ^in is expressed as 



^in(y) = r ^i£(£M^iMd,, < y < /_(/_). (3.27) 

Jo COs6in(s) 
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Denote by Y[^ the inverse function of ^in, i.e. 

YinW = KnW, < VJ < lUin- 

At the upper wah, the potential function is 

(p{x,f-{x)) = ^^{x), /_ < X < 
satisfying <!>-(/-) = C-, *^-(0) = and 



(3.28) 



$_ x) 



Thus, $_ is expressed as 



<I>_(x) 



cosG_(x 
^_(s) 



> 0, l_ <x <0. 



-ds, l^ < X < 0. 



/Q cos0_(s) 
Denote by X_ the inverse function of $_, i.e. 

Therefore, the subsonic-sonic flow problem in the potential plane is formulated as follows 
d^A{q) d'^B{q) 



(3.29) 



(3.30) 



dA[q) 

dip 
dq 
dip 
dB{q) 



+ 



9-02 



(C-.V^) 



(99,0) = 0, 



0, 



e;„(y)cosei„(y) 



('/',"T.in) 



^Uy)p{^l{y)) 



G'_(a;)cose_(x) 



y=yUi') 



dil) J2-{x) x=X-{tf)' 

g(O,'0) = c*, 

^in(?/)=9(0,V) , ^ ,,, 

"^-(x) =g((/?,min) 

73=<l?_(a:) 

Remark 3.2 A natural boundary condition at the inlet seems to be 



(<^,V^)G(C-,0) 


X (0,min), 


(3.31) 


V^ G (0,min), 




(3.32) 


¥'G(C-,0), 




(3.33) 


'/'G(C-,0), 




(3.34) 


V' G (0,min), 




(3.35) 


ye [0, /„(/„)], 




(3.36) 


x€ [/_,0]. 




(3.37) 



dip 



/(C-,V)Kg'(C-,V')) 



(3.38) 



instead of (|3.32p . However, since — ( — ;— ttt ) < m g € (0,c*) anc? — 0(„(yin(V')) cos 0in(11n('0)) ^ 

dq^qpiq^)^ 
l/(2i?o) /rom (j3.24p . ii seems difficult to obtain the well-posedness of the problem (|3.3ip . (j3.38p . 

(|3.33p - (j3.35p . Therefore, we prescribe the boundary condition at the inlet by (|3.32p . 
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The problem (J3.3ip - (j3.37p is a boundary value problem for a second-order quasilinear degenerate 
equation in a rectangle with two free parameters. Furthermore, the degeneracy is characteristic and 
the boundary conditions (|3.32p and ()3.34p are nonlinear, nonlocal and implicit. As far as we know, 
there is no known theory to show the existence and uniqueness of solutions even if the problem is 
uniformly elliptic. Thus we will use a fixed point argument to prove the existence of solutions to the 
problem (I3.3ip - (l3.37p as follows: For given ^\^ and S-, we solve the problem (|3.3ip - (13.35p . then 
define a mapping by (I3.36P and (I3.37p . To this end, one must determine and control precisely the 
rate of the solution tending to c* in order to solve the problem (j3.3ip - (|3.37p . However, there is no a 
background solution to this problem to suggest what the rate is. And we will carry out some precise 
elliptic estimates in this paper to determine and control this rate. 

Since (I3.31|) is degenerate, it is more convenient to introduce weak solutions to the problem (|3.3ip - 
()3.35p . After determining the rate of the solution tending to c*, one can establish the regularity of 
weak solutions. 



Definition 3.1 A function q € L°°(((^_,0) x (0, rjiin)) is said to he a weak solution to the problem 
(lOTD - dSlSSD . if 

< inf q < sup q < c^ 

(C-,0)x{0,mi„) (C_,0)x(0,min) 



and 



o2 A o2 c 



dip' 

" 9;„(j/)cos9in(^/) 
^in(y)p(^fjy)) 



y=Yinii>) 






9'_(x)cos9_(rc) 



^^(x) 



x=X- (if) 



C{(p,min)d(p = 



for any ^ € C^([C-,0] x [0, min]) with 



|<"») 



(C-,0) 






(C-,0) 



and 



|<c- 



(O.min) 



e(o,' 



(O.min) 



0. 



Though the condition (j3.22p seems quite strict, yet, it is almost necessary for the existence of C^ 
solutions as shown in the following remark. 



Remark 3.3 If the problem ()3.3ip - ()3.35p admits a C^([C-,0] x [0,r7iin]) solution, then 

f'^{x) = 0{x^), L < X < 0. 
Proof It follows from HSMh and (l3:35|) that 

2 



^"(e.,(|(„,,4)%B'(„)0...-.>.o" 



(0.*) + B"(c.)(^(0,«) 



and 



0(°-« = »' 



d'^q 

dip'^ 



(0, V') = 0, < ^ < mi, 



which imply 



|(0.^) = 0, 



< ip < mi, 



(3.39) 



0, < V < "li: 
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and thus 

fin 

-{ip,7P)=0{ip), q{ip,7P)=c,-0{ip'), (99,^) G (C-,0) X (0,?nin). 



dip 
Hence 

(v?,V')G(C-,0) x(0,min), 
which, together with ([33T]) . ([333]) and ([33^1) . yields 

- -[^,'4')d^ = 0{ip^), 99G(C-,0). (3.40) 



Then, ([H:HU]) follows from ([2301), e„(0) = 0, ^_(0) = c, > 0, X_(0) = and X'_{{)) = 1/c, > 0. D 

3.3 Formulation of the sonic-supersonic flow problem in the potential plane 

A sonic-supersonic flow is expected in the divergent part of the nozzle ^2+. In this subsection, we 
formulate this sonic-supersonic flow problem in the potential plane. 
Assume that /+ e C^{[Q,l+\) (0 < /+ < 1) satisfies 

5i,+x^+ </+(x) <(52,+x^+, /"'(x)>0, 0<x</+ (3.41) 

with positive constants A+, 5i^+ and ^2,+ such that A-|- > 2 and 5i^+ < (^2,+ . 

Remark 3.4 The second condition in (j3.4ip can he relaxed by 

f"'{x) > -ujx^+~^, < X < /+, 

where uj is a positive constant depending only on 7, ni, A+, (5i^+ and (^2,+ (see Remark \5.S\ for the 
detail). 



Set 






Q{^,^)=A{q{ip,i:)). 


Then, (]2.3]) takes the form 


Q^^ + {K'{Q)Q^)^ = 0, 


where 


K{s) = B{A^\s)), s<0 


and 


K'is)^{l-'^^\Y{'-'~2'^f''~"^' 


1 <0' 

q=A-\s) 


K"is)H,^l),il-'^,\Y{'-'~2'^1 


3/(7-l) + l 


1 <0' 

q=A-\s) 



(3.42) 



s <0, 



s < 0. 



41 



Since lini K (s) = — oo, (j3.42p is a singular nonlinear wave equation in the sonic-supersonic region. 
We now transform (j3.42|) into a first order system to use the method of characteristics. Set 

U{ip, i^) = Q^{(p, ip), V{ip, ip) = Q^{ip, tp). 
Then, (j3.42p is transformed into the following system 





U^ - {b{Q)V)^ = 0, 








V^-U^ = 0, 






i.e. 










(U\ / -b{Q)\/u\ /O -p{Q)Q^^ 






where 










b{Q) = -K'{Q) > 0, piQ) = -K"{Q) > 0, 


g<o. 




Set 









(3.43) 



Multiplying (j3.43p by R~^{Q) on the left side, one gets 



which is equivalent to 



W\ / Vm) Ww^\ / -p{Q)Q^ \ U\ dR~\Q) W 



( VW) \ afi-i(Q) /^^ 

^V -76(0)7 ^V^ Uy 



where 



\z) V V \-u- Vmv) [ -Q^ - VmQ^, 

Therefore, (|3.43p becomes 

W^ + b'^\Q)W^ = ^b-\Q)p{Q){Q^ - b^/HQ)Q^){W + Z), 
Z^ - b'/\Q)Z^ = ^b-\Q)p{Q){Q^ + b'/\Q)Q^){W + Z). 
The flow angle on the upper wall is 

0+(x) = arctan/^(x), < x < /_|_. 
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Let the speed of the flow on the upper wall be given by 

q{x,U{x)) = J2+{x), 0<x</+. 

Then, the potential function on the upper wall is 

(^(x,/+(x)) = ^+(x), 0<x</+ 



satisfying <l*+(0) = and 
Thus, $+ is expressed as 



cos6+(2;) 



*-«=/„ i;^)'"- "S-Si.. (3.44) 



Denote by X^ the inverse function of <!>+, i.e. 

X+iif) = <^+\ip), < (/. < C+, (3.45) 

where 

C+ = $+(/+). (3.46) 

Therefore, the sonic-supersonic flow problem in the potential plane can be formulated as 



Qw - iKQ)Q^h = 0' 


((^,V^)G(0,C+)x(0,m), 


(3.47) 


Q(0,V')=0, 


V'G (0,m), 


(3.48) 


Q^{o,^p) = o, 


V'G (0,m), 


(3.49) 


Q^{^,0) =0, 


9'e(o,C+), 


(3.50) 


eV(x+(v9))cose+(x+((/,)) 

^'^^'^'"'^ 6(Q(v.,m)))^;HQ(^,H) ' 


V^G(0,C+), 


(3.51) 


^+(x) = A;i(Q($+(x),m)) 


xe [o,Z+] 


(3.52) 



or equivalently as 

W^ + b^/\Q)W^ = h-\Q)p{Q){Q^ - b'/^{Q)Q^){W + Z), (yp, V) € (0, C+) x (0, m), (3.53) 

Z^ - h''\Q)Z^ = \h~\Q)p{Q){Q^ + h^'\Q)Q^){W + Z), {^, ij) e (0, C+) x (0, m), (3.54) 

M^(0,V') = 0, '(/'G(0,m), (3.55) 

Z(O,'0) = O, 2pe{0,m), (3.56) 

W{ip,0) + Z{ip,0) = 0, <^G(0,C+), (3.57) 

wr ^ ^ zr ^ 2eV(X+(y.))cos9+(X+((^)) 

W^((/p,m)+Z((/j,m) = — — 1 r-, v?e(0,C+), (3.58) 

6V^(Q((p,m))A+ {Q{ip,m)) 

Q^{ip,i^) = l{W{ip,iP)-Z{ip,iP)), ((/p,V)G(0,C+)x(0,m), (3.59) 
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Q(0,V)=0, 



ip € (0, m), 

xe[o,l+]. 



(3.60) 
(3.61) 



The problems (|3.47p - (j3.5ip and (j3.53p - ()3.60p should be regarded as initial boundary value prob- 
lems, where v^-direction plays the role of the time. Since lim b(Q) = +00, (I3.47p . (|3.53p and (13.54p 

are singular at sonic points. Furthermore, as shown in Theorem 12. 5| all positive and negative char- 
acteristics from any sonic point are the same line {0} x [0, m]. That is to say, the singularity at the 
sonic curve is so strong in the sense that there is not any characteristic moving from a sonic point to 
the supersonic region. The source terms in (]3.53p and (j3.54p are also singular at sonic points because 

lim b~ {Q)piQ) = +00. Moreover, the boundary conditions (|3.5ip and (|3.58p are nonlinear, nonlocal 
Q->0- 

and implicit. We will use a fixed point argument to prove the existence of solutions to the problem 
(J3.53p - (j3.6ip as follows: For given Q and =S_, we first solve the problem (J3.53p - (j3.58p . then define a 
new Q by solving the problem p.59p . (I3.60p and a new <S_ by (I3.61|) . To get a fixed point, one must 
choose a space with a precise behavior near the sonic curve for Q and then show that the new Q also 
belongs to this space by some elaborate calculations and optimal estimates. 

Since (|3.47p . (j3.53p and (j3.54p are singular, we first seek weak solutions, which are defined as 
follows. By more complicated and precise estimates, one can get smooth solutions. 



Definition 3.2 A function Q is said to be a weak solution to the problem (j3.47p ^ p.51|) if Q ^ 
C°'i([0,C+] X [0,m]) satisfies dOS]) and 



{Q^{ip, i^)i^{^, -0) - KQi'f^ '0))Q^(¥', ^)C^(¥', '>P))dipdTp 



<;+eV(x+((/,))cose+(x+((^)) 



A~\Qiip,m)) 



^{ip,m)dip = 



for any ^ G C^{[0,(+] x [0,m]) with ^(C+, ■ 



(0,m) 



0. 



Definition 3.3 A triad of functions {W,Z,Q) is said to be a weak solution to the problem ()3.53p - 
msmh ifW,Z e L°°((0,C+) x (0,to)) and Q G C°'i([0,C+] x [0,m]) satisfy 
(i) For any ^,77 G C^([0,C+] x [0,m]) with 



(0,m) 

it holds that 



viC^ 



(0,m) 



0, e(-,o) 



(o,C+) 



-^(•,0) 



(o,C+) 



^(•,m) 



(o,C+) 



-ri{-,m] 



(o,C+) 



(W^^ + b'/\Q)WC^ + -b-^'\Q)p{Q)Q^Wi 
+ \h'\Q)p{Q){Q^ - b''^{Q)Q^){W + Z)iy^d^ 
+ j J [Zv^-b'/\Q)Zr,^--b-^/\Q)p{Q)Q^Zr^ 

+ h-\Q)p{Q){Q^ + b^/\Q)Q^){W + Z)riyipd^ 

^+ 0V(X+(yp))cos9+(X+((^)) 
A-\Q{^,m)) 



S,{(p,m)dip = 0; 
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(ii) (j3.59p holds in the sense of distribution; 
(iii) Q satisfies (|3.60p . 

As the discussion for the subsonic-sonic flow problem in Remark 13.31 the condition (j3.4ip on the 
wall of the nozzle is almost necessary for the existence of C^ solutions. 

Remark 3.5 If the problem (I3.47p - ()3.51|) admits a C^([0,C+] x [0,m,]) solution, then 

f'l{x) = 0{x^), 0<x<l+. 

3.4 Main results of the transonic flow problem in the de Laval nozzle 

We now state the main results on the transonic flow problem. Since the sonic curve of the transonic 
flow studied in the paper is located at the throat of the nozzle, the results of the transonic flow can 
be stated for the subsonic-sonic part and the sonic-supersonic part separately. 

For the subsonic-sonic flow in the convergent part of the nozzle, the first result concerns the 
existence and regularity estimates of such a solution. 

Theorem 3.1 Assume that /_ G C^{[1_,0]) satisfies (KWi and g G C^{[0, f-{l-)]) satisfies ^JM 
and (|3.24p . There exists a positive constant Jq,- depending only on 7, m, A_, 6i^~, ^2,- ^iT'd e(-)) 
such that for any — (5o,_ < ^- < 0, the problem (J3.3ip - (j3.37p admits a weak solution q € C°°((C-,0) x 
(0,min)) n C^{[C-,0) X [0,min]) n C([C_,0] X [0, min]) with mjn = m satisfying 



dq 



iip,^P) < Cl,„(-C-)^"^/'(-y^)^/^ {ip,^P) G (C~,0) X (0,m), 



(3.62) 



A{qi^',,p'))-A{qi^",r)) <C^,_i\^'-^"\'/' + \^l^'-r\), (</^', V'), (V'", V'") G (C-,0) x (0,m), 

(3.63) 
c.-C3,-(-v^)^-/2+i <g((/p,V) <c*-C2,-(-v^)^-/2+i, {^,^l^) G (C_,0) x (0,m), (3.64) 

where Ci-{i = 1,2,3)(C2,- < C3-) are positive constants depending only on 7, m, A_, 61- and 
(52,-. Moreover, 

(i) // /_ satisfies 



\f'l'ix)\ < (53,_(-x)^-/^+i/2^ ^_ < X < 



(3.65) 



with a positive constant 6^- additionally, then any weak solution satisfies q G C^([C-,0] x [0, ?7T,in]) 
and 



dq 

dip 



{^,i^) <C4,-(-(^) 



A_/4+l/2 



dq 



{ip,i;) < C4,-(-(^)^-/2+i, {ip,^p) G (C-,0) X (0,m), (3.66) 



where the positive constant C^^- depends only on 7, m, A„, (5i,_, 82 ~ and 5^^. 
(ii) ///_ G C"^([/_,0]) satisfies (IHTHHIl and 



\f':^\x)\<5. 



4,-, 



/_ < X < 



(3.67) 



with a positive constant 84,- additionally, then any weak solution satisfies q G C^'^{{(^-,0] x [0,m]) 
and 



d'^q 
dif"^ 



iv,i^) 



<c. 



d'^q 
dfdip 



i^,^) <c,A-^)'-^'^'^', 



d\ 

dtjj'^ 



(ip,^) <C5,-(-¥.)^-/2+i, 
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((^,V')G(C-/2,0)x(0,m) (3.68) 

with a positive constant C^- depending only on 7, m, A_, 5i,-, ^2-, ^2,- o.'^'d (54,-. 

(iii) If g ^ C^'"([0, /_(/_)]) and /_ G C^'"([/_,0)) for a number a G (0, 1) additionally, then any 
weak solution satisfies q G C^'"([C-,0) x [0,m]). 

The next result yields the uniqueness. 

Theorem 3.2 There exists at most one solution ip G C^'^{i^-) to the problem p.l0p - (|3.15p . 

For the sonic-supersonic flow in the divergent part of the nozzle, we start with the existence and 
regularity estimates. 

Theorem 3.3 Assume that f^ G C^([0, /+]) satisfies (I3.4ip . There exists a positive constant 60,+ > 
depending only on 7, m, A+, (5i,+ and 62,+ , such that for any < /+ < (5o,+, the problem (j3.47p - ()3.52p 
admits at least one weak solution Q G C*^'^([0, C+] x [0,m]) satisfying 

-C2,+(^^++2 < Q((^, ^) < -Ci,+/++2, (c^, V;) G (0, C+) X (0, m), (3.69) 

(3.70) 

wit/i positive constants Ci,+ and 6*2,+ (Ci,+ < 6*2,+) depending only on 7, m, A+, 5i,4. and (52,+ . 
Furthermore, assume that /+ satisfies 

I/+ (a;)| < <53,+x^+-\ 0<x</+ (3.71) 

wit/i a positive constant (53,+ additionally. Then there exist positive constants 5o,+ '^'^'^ ^"3,+, 6oi/i 
depending only on 7, m, A+, (5i,+ , (52,+ and (53,+, such that if < l^ < (5o,+, ^/len the problem 
(f3?T7l) -- ([332]) admite a^ /easi a solution Q G C^'^([0,C+] x [0,m]), w/iic/i satisfies ([3^ . (fSTTOl) and 

|Qw(¥'>V')l < C3,+(/.^+, |Q^V'(</''V')I < C3,+v.^^+/^+l/^ IQ^^(¥',V')I < C3,+(/p=^^+/2+\ 

((/p,V)G(0,C+)x(0,m). (3.72) 

The following theorem shows the uniqueness of weak solutions. 

Theorem 3.4 Assume that /+ G C*{[0,1+]) satisfies ^JT\f and (f3TT|) . T/ien the problem ^J7\f - 
()3.52p admits at most one weak solution Q G C'^'"'^([0, (;^+] x [0,m]) satisfying ()3.69p and (I3.70|) mi/i 
any given positive constants Ci^-\- {i = 1, 2). 

Let us connect the subsonic-sonic flow and the sonic-supersonic flow to get a global smooth 
transonic flow. 

Theorem 3.5 Assume that /„ G C"^([/_,0]) satisfies iK22\i . (IHIBSll and JKm . g G C3'<^([0, /_(/_)]) 
satisfies ([3:23]) and (fSlMD . i(;/ii/e /+ G C^([0,/+]) satisfies ([OT]) and (fSTTI) . // -Z_ and /+ are 
sufficiently small, the problem (j3.3ip - (|3.37p admits a unique solution g_ G C"'^'^([(^_, 0] x [0,m]) sat- 
isfying (j3.62p - (13.64p . (13.66P and (j3.68p . w/ii/e t/ie problem (J3.47p - (j3.52p admits a unique solution 
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(5+ G C^'^{[0,C+] X [0,?7i]) satisfying (|3.69p . (j3.70p and (|3.72|) . Connect q^ and Q^ in the following 
way 

f^-(9',V'), (99,^^) G [C-,0] X [0,m], 

'^^'^' ^"\a;1(Q+(^,V')), ((^,V')G[0,C+]x[0,m]. 

Then, q £ C^'^([C_,C+] x [0,?7i]) is a solution to (I2.3p . Furthermore, — — = — — = ai ifte sonic curve 

dip dip 

Q^q 
and ——77 is continuous at sonic points. Therefore, q satisfies (12. 3p at the sonic curve, i.e. 
oip'^ 






lim ( „'!,'^'' + -^ ^^^7^^^^ 1 =0, < V' < m. 



In the physical plane, Theorem 13.51 can be stated as follows. 

Theorem 3.6 Assume that /„ G C"^([/_,0]) satisfies <^M), (l3:65D and <^M), 9 € C3'°([0, /_(/_)]) 
satisfies IK2^ and ([321]), i^/i«/e /+ G C^([0,;+]) satisfies pTiTD anc? (ITtTD . // -L and /+ are suffi- 
ciently small, the transonic flow problem (I3.2p -( f3r9]l admits a unique classical solution ip G C^'^(Q), 
which satisfies (|3.62p - ()3.64p . ()3.66p and ()3.68p /or g = IV99I in Q_, w/iiZe satisfies (j3.69p . (|3.70p and 
(f3T2]l /or Q = A(|V(^|) in 0+. 




X 

Figure: the smooth transonic flow in the de Laval nozzle 

Remark 3.6 In Theorem \3.6\ we get a C^'^ transonic flow. However, this transonic flow pattern 
is strongly singular in the sense that the sonic curve is a characteristic degenerate boundary in the 
subsonic-sonic region, while in the sonic- supersonic region all characteristics from sonic points coin- 
cide, which are the sonic curve and never approach the supersonic region. 

Remark 3.7 In the de Laval nozzle 0,, there is also a subsonic-sonic-subsonic flow, whose sonic curve 
is located at the throat of the nozzle. 

Remark 3.8 For the nozzle Q, if there is a flat part between the convergent part and the divergent 
part, there is also a transonic flow, which is subsonic in the convergent part, sonic in the flat part, 
while supersonic in the divergent part. 

Remark 3.9 The geometry of nozzles is important for smooth transonic flows. As mentioned in 
Remarks O and EH 

f"{x) = 0{x^), l^<x<l+ 
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is necessary for a C^ transonic flow whose sonic curve is located at the throat of the nozzle. Moreover, 
in \20^ . we get a continuous subsonic- sonic flow in a convergent nozzle with straight solid walls. 
However, this subsonic-sonic flow is singular in the sense that while the speed is continuous yet the 
acceleration blows up at the sonic state, and there is no way to extend it to be a transonic flow or a 
subsonic-sonic-subsonic flow. 

4 Subsonic-sonic flows in the convergent part of the nozzle 

In this section, we will establish the well-posedness of the subsonic-sonic flow problem in the conver- 
gent part of the nozzle. For convenience, we abbreviate A_ and 5^- by A and 5^ for k = 1,2,3,4, 
respectively, and use Cj, c, /Xj, kj, Mj to denote generic positive constants. Furthermore, a parenthesis 
after a generic constant means that this constant depends only on the variables in the parentheses. 
Owing to (l3:22]) . /_(/_) and Rq satisfy 

^^_l„2/(7-l) ^ _^_(^_^ ^ ^^_l_2/(T,-l) ^ ^^^ _^^ ^ ^^^^_^ ^ Q 
and 



1+2/(7-1) . 1+2/(7-1) 



(A + l)(m + 52ci+'/(^~'))y^fTT ^0 (A + l)m 

with m given by (j3.ip . 

In order to solve the problem ()3.3ip - (j3.35p by a fixed point argument, one needs to specify £2\a 
and J2- in advance as follows: =Sin G C'^'^([0, /_(/_)]) satisfies 

max{|,c, -C2(-/-)^/2+i| < ^.^(y) < max{|,c. -Ci(-/„)V2+i}, |^(^^(y)| < i, 

< y < /-(/-) (4.1) 

and ^_ € C([/_,0]) nC^/2Q/_^;_/2]) satisfies 

C < ^-(X) < C, L < X < 0, [^-]l/2;(«_,Z_/2) < Co, (4.2) 

where Cj (i = 0, 1,2) and c(Ci < C2, c < c*) will be determined. Then, rriin, C-i ^in ^'^d <!>'„ satisfy 
mi = ^7nc;'/(^"'V(c^/4) < min < ma = c./5(c2)(mc;'"'/(^~') + 52)^<5i + 1, 



c*a/5| + 1/_ <C- <cl-, 



\c.p{cl/A) < ^^{y) < c,p{cl)^5l + l, < y < /_(/_), 



c<^'_{x)<c^j5l + l, /_<x<0. 
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Moreover, one can verify that 



^^(_(^_)A+i < Qhi(^)cosQin(y) 



^Uy)p{^l{y)) 



y=yu^) 



<IJ2{-C-r^\ ^G(0,min), 



'@{^iy)coseUy) 






y=yinW 



X ^ 0'_(x) COS 0_(x 



^_(x) 



< ^^2(^2 + l)(e(C-) + e)(-C-)'^/^ ^ G (0,mi,), 



x-=X_((p) 



(4.3) 
(4.4) 
(4.5) 



where ii\ and /i2 (/^i < /^2) depend only on 7, m, A, 5i, 82 and c. 
4.1 A comparison principle 



We will state a comparison principle for weak solutions to the problem (J3.3ip - (j3.35p . More generally, 
instead of (j3.35p . one imposes the following boundary condition 



g(O,-0) =c, V e (0,min), 



(4.6) 



where < c < c*. Weak solutions, supersolutions and subsolutions to the problem (|3.3ip - (|3.34p . (j4.6p 
are defined in the following sense. 

Definition 4.1 A function q G L°°((^_,0) x (0, ?7iin)) is said to be a weak supersolution (suhsolution) 
to the problem (l3:3T|) - (l33iD . (jMD, if 

< inf q < sup q < c^ 

{C-,0)x{0,mi„) {C_,0)x(0,min) 



and 



t-niin 

C- ^0 






Mliv. V'))^(<^, ^) + B{q{^, ^)) 5wi^('^' ^) ) ^'^^^ 



- e;^(y)cosGin(y) 



+ 



Jo d<f Jo ^in(y)p(^i';(y)) 

e'„(x)cose_(x) 



y=YU^P) 



e(C-,V')# 



C- 



^_(x) 



x=X_ ((p) 



^{(p,min)d(p < (>)0 



/or any nonnegative function S, G C^([C-,0] x [0,min]) with 



i'-») 



(C-,0) 9-0 



= and -^(C 
(C-,0) d^p' 



(O.min) 



e(o,' 



(O.min) 



0. 



Furthermore, q € L°°((<^_,0) x (0, min)) is said to be a weak solution to the problem ()3.3ip - ()3.34p . 
if q is both a weak supersolution and a weak subsolution. 

Let q € L°°((C-, 0) x (0, min)) be a weak solution to the problem (fSrsB - dS^iD . (jM]). If q satisfies 



sup q < Cf, 

{C-,0)x{0,mi„) 



(4.7) 



theng G C°°((C-,0) x (0, rriin)) nC^([C-, 0] x [0,min)) nC([C-,0] x [0,?nin]) due to the classical elliptic 
theory. 

It follows from the proof of Proposition 3.1 in |2U] that the following comparison principle holds. 
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Lemma 4.1 Assume that q+ € L°°((C-,0) x (0,77iin)) and q- G L°°((C_,0) x (0,77iin)) are a weak 
supersolution and a weak subsolution to the problem (I3.31j) - (l3.34p . (j4.6p . respectively. Then 

q+{^,ip) > q-{(p,ip), {(p,^^) G (C~,0) X (0,min). 

4.2 Solutions to elliptic boundary problems 

We now study solutions to the problem (j3.3ip - (|3.34p . (j4.6p with < c < c*. Due to Lemma HTTl such 
a solution is unique. First we show the existence for sufficiently small c. 

Lemma 4.2 There exists ki G (0, 1] depending only on 7, m, X, 61 and 82 such that for any —ki < 
C_ < 0, the problem (J3.3ip - (l3.34p . (j4.6p admits a unique solution with (j4.7p and 

sup g(C-,-) < c. - (e(C-) +e)(-C~)^/'+' (4.8) 

(0,min) 

ifO<c< c,/2. 
Proof. Set 

Choose Ki G (0, 1] such that 



en = mm < ^ , mi ^,, , > . 

Il2m| c,/2<.<5c./6S'(s)J 



2 c* c* j^ ^ 2eomi 



Then, 



Kf < — , Ki < — , < < ^^^ inf B'(s). 



2 



is a supersolution to the problem (j3.3ip - (13.34p . (14. 6p if < c < c=f/2. In fact, 

1 5 

-c* < q+{v,'il^) < -c*, (99, V') G [C-,0] X [0,min], (4.9) 

2 b 

^^(^'^) + ^^(^'^)<^'('^+(^'^))^(^'^)+^V(v',^))^(^,^) 

= - 2A'{q+{ip, ^)) + 2eoB'{q+{ip, ^)) 
<0, (99,^)G(C-,0)x(0,7nin), 

^^^^+^-(C-,^) = -2(1 + C-)^'(g+(C-,V')) < 0, V e (o,mi„), 



5^(9+) 



((/p,0)=0, v^G(C-,0), 



dijj 
{(p,min) = 2eominB'{q+{ip,min)) > fj.2{-^)^, ^ e (C-,0). 
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Additionally, it is not hard to verify that 

q-{ip,ip) = Az\A{c) + fi2^), {ip,ip) E [C-,0] X [0, m;, 



is a subsolution to the problem (J3.3ip - (j3.34p . (j4.6p . Then, one can get the existence and the uniqueness 
of the solution to the problem (I3.3ip - ()3.34p by a standard argument in the classical elliptic theory. 
Finally, (gS]) follows from (gj]) and (e(C-) + C-)(-C-)^/^+^ < 2ki < c,/12. D 



Below, we investigate properties of the solution to the problem (j3.3ip - (|3.34p . ()4.6p with < c < c*. 
As ithe proof of Proposition 3.2 in [20], we introduce two functions. The first one is 



Clearly, E G C°°((-oo,0)) satisfies 



E{s) = A{B~\s)), s<0. 



^'(.)Hi-^.^)(i-^.^^^^"^"'" 



E"is) = -i^ + l)q^l 



7-1 2\-3/(7-l)-2 



q=B-^s) 

<o, 



>o, 



E"'{s)=-i^ + l)q\i + 3q^) 1 



q=B-Hs) 
7-1 2\-'^/(t-^)-3 



g=B-i(s) 



<0, 



s <0, 
s < 0, 
s < 0. 



The other function is 

G{s) = E'{E"\s)), s<0, 
which satisfies G G C°°((-oo,0)) and 
E"{E~Hs)) E"{t) 



G'{s) 



G"{s) 



E>{E-^{s)) E'{t) 
E"\t)E'{t) - {E"{t)f 



t=E-^{s) 



<o. 



'dW-i 



(E'it)) 



i=£-i(s) 



< 



{E"it)) 



'm^3 



(E'it)) 



t=_B-i(s) 



<o, 



s<0. 



s<0. 



Lemma 4.3 Assume that c*/3 < c < c^: and q, satisfying (j4.7p and ()4.8p . is i/ie solution to the 
problem (fSrsT]) - (l33i] l . gSD- ^/len 



g((^,^) > Al^(A(c,/3) + /i293), ('/^,'0) e [C-,0] X [0,miJ, 



9g 



^(<^,V) <A.3(C2 + l)(-C-)'-'/'(-<^)'/^ (<^,V^)G[C-,0]x[0,m, 



(4.10) 
(4.11) 



and 



^(g((^', VO) - ^(9(<^", O) < M4(C2 + 1)(|<^' - v^'f /2 ^ j^/ _ ^//|)^ 

(99',V''),(9'",Oe[C-,0]x[0,min], 
where fii = fJ-ii^j, ?n. A, 5i,52,c) for i = 3,4. 



(4.12) 
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Proof. First, (|4.10p follows from Lemma |4. II and the fact that 

q-{ip,ilj) =AZ^{A{c,/3)+fi2^), i^,i;) G [C-,0] x [0,mi^] 
is a subsolution to the problem (I3.3ip - (l3.34p . (I4.6p . To prove M.llh . one can adapt the proof of 



{^,ij), ((/J,V)€ [C-,0] X [0,mi: 



Proposition 3.2 in [20] as follows. Set 

dB{q) 



Then z solves the following problem 



dip 



d z d z dz dz 



(c/^,V')G(C-,0) x(0,min). 



dz 



(C-,^) + e(^)z(C-,^) 



i?'(<?(C-,v))/e;jy)cosei„(2/) 



dip 
z(v;,0) = 0, 

z((/?, mill) = 
z(0,V)=0, 



A>{q{c.,ip))\^Uy)p{^M) 



y=Yinii>) 



e'„(x)cose_(a;) 



^_(x) 



x=X_((/3) 



^ G (0,min^ 
¥'€(C-,0), 

^e(C-,o), 

V' G (0, min' 



(4.13) 
(4.14) 
(4.15) 
(4.16) 
(4.17) 



where for {(p,Tp) G (C-)0) x (0, min). 



ji{ip,iP) = G{Aiq{ip,i;))), hi^,^) = 2G'iA{qi^,^))) 



dA{q{ip,iP)) 
dip 



M9. W = -G'(A(,(^, *)))^^%^ , M9, « = a"(A(M^. m' "•"""• *' ^ ' 



and 



eW 



dijj 



E"iqiC_,i;)) e[^iy)cosQUy) 



y=YinW 



dip 
<0, ^G(0,min). 



It follows from M.lOh . ()4.8p and ()4.4p that there exists a positive number fl-^ depending only on 7, m, 
A, (5i, 62 and c such that 

B' {q{C.,^))fe[^{y) cos QUy)^ 



A'iqiC^,iP))\jS,4y)p{J2l{y)) 



y=YUi>) 



<^A3(C2 + l)(-C-)^-\ ^G(0,min). 



Note 



^J,i^,^l;){-^)-^/^-j,{^,^|;){-p,)y^ >V-Ji(¥',V')i4(</:',V')(-<^)-'/' 

> - ^J2(V5,V')(-'/')"'/', (V',V') G (C-,0) X (0,min). 
One can verify that 

z±((^, V) = ±fl3{C2 + l)(-C-)^~'/'(-'^)'/', (<^,V') € [C-,0] X [0,min] 



are super and sub solutions to the problem (J4.13p - (j4.17p . Then, (|4.1ip follows from the classical 
comparison principle. Finally, one can get (|4.12p from (j4.10p and (j4.1ip by a standard process as in 
the proof of Proposition 3.2 in [20]. D 
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4.3 The weak solution to the degenerate elliptic boundary problem 

We are now ready to solve the problem (j3.3ip - (|3.35p . First we determine c by 

c = Az\A{cj3)-^i2) G(0,c*/3). (4.18) 

Proposition 4.1 There exists K2 € (0,ki] depending only on 7, m, A, 61, 62, e(-) and C2 such that 
for any —K2 < C~ <0, the problem ()3.3ip - (l3.35p admits a weak solution q £ C°°((C-,0) X (0, min)) n 
CH[C^,0) X [0,min))nC7([C_,0) X [0,min]) satisfying 



dq 

dip' 



i^,i;) < ^3(C2 + l)(-C-)^-^/'(-c^)^/^ ((/',V') G [C-,0] X [0,m, 



(4.19) 



A{q{^', iP')) - A{qi^",r))\ < IIa{C2 + IW - V"\''^ + W - <|), 

((/.',V''),(9'",Oe[C-,0]x[0,mi,], 



(4.20) 



where /.is and ^g (^5 ^ ^e) depend only on ^, m, A, (5i and 52. 

Proof. Set 

^ = |cj3 < c < c* : the problem ([3:31]) - ([33i|) . (jOj) admits a solution g^ with ([32]) and ^^\. 

According to Lemmas 14. Il and l4.2[ S" (^ 0) is an interval. We claim that S = (c^,/3, c*) if —K2 < C- < Oi 
where K2 S (0, ^i] will be determined. 

Assume that c € <#". It follows from Lemma 14.31 that qc satisfies (|4.1ip and (|4.12p . Integrating 
(j3.3ip and (|3.32p over (0, rriin) with respect to ip and using (|3.33p and (j3.34p . one gets that 



A{qc{Lp,il)))d'>lj 



d'^^ Jo 

i^— I A{q^{ip,iP))diP 



e'_{x)cose-{x) 



^_(x) 



x=X_(</p) 



v^iC-,0), 



(4.22) 



v=C- 



-'"G;jy)cosein(y) 



^in(y)p(^!L(y)) 



J/ = >1nW) 



#■ 



Direction calculations yield 
/■o e'__(x)cose_(x 

and 



^_(x) 



X=X_(l/3) 



(ic/5 



° 9^_(x)cose-(x) ^, ^ 



e'_(x)dx = -e_(-/_ 



-'■^e;j2/)cosei„(2/) 



^Uy)pi^Uy)) 



/^-('-)ey^)cose^ 

f-ii-) 

e;„(y)dy = ein(/_(L)). 
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Therefore, 






95=0 



0. 



(4.23) 



It follows from ([422]), (1423]) and gSJ that 



A{qc{ip,'il)))dip = mi^A{c) 



/-o a' 



¥> •''/' 



e'_(x)cose_(x) 



=g_(x) 



x=X_(s) 



(is(i(^, (/? G [C-)0]) 



which, together with ()4.5p . leads to 

in 

A(gc(v^,V'))# < "T-in^lc) 



^1 



(A + l)(A + 2) 



-V,) 



A+2 



</'G[C-,0]. 



Therefore, there exists a Lipschitz continuous curve L : ip = uj{ip) (C- < (/^ < 0) in [C-,0] x (0,min) 
satisfying 



A(gc('/',w(9j))) < 



/^i 



(_^)A+^, (^e[C_,o]. 



min(A + l)(A + 2) 
Thus, there exists Mq = Mo(7, m. A, 51,62) such that 

which, together with (14.110 . leads to 

qc{C~,i^) < c, - Mo{-C^)^/^+^ + ^13{C2 + l)m2{-C-)\ ^ e [0,m: 
Decompose (C-)0) x (0, mjn) into Vli and ^^2 by the curve L, where 

Oi = {((^,^) G (C-,0) X (0,min) :^<w((^)}, !^2 = {(¥', V') e(C-,0) x (0, mi,,) ■.tP>u{^)]. 
Then, qc solves the following two problems 



(4.24) 



(4.25) 



+ 



0, 



ge(C-, V') < c* - Afo(-C-)^/'+' + //3(C2 + l)m2(-C-)\ 
dqc 



dij) 



((^,0)=0, 



ge(</^,c^(v?)) < c* - Mo(-^)^/2+\ 
<?c(0,^) = c, 



and 



d^A{q,) ^ d^B{q, 



0, 



(C_, V-) < c* - Mo(-C-)^/'+^ + /U3(C2 + l)m2(-C-)\ 

(99,^(99)) <c,-Mo(-(^)^/2+i^ 



^€(0,^(C-)), 

^G(C-,o), 
'/'e(C-,o), 

^G (0,^(0)) 

{ip,i{j) € i72, 

V' € (a;(C-),min) 

v'e(C-,o), 
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dBjqc) , . eL(x)cos9_(x) 

qc{0,ip) = c, V S (a;(0),min). 

Set _ 

and 

with U5 = — -. Note that qi and ^2 satisfy 

2(1 +1712) 

^0^{v,i^) + ^^i^,i^)<o, (v.,^)GOfc, k = l,2. 
Choose a number K2 G (0, ki] depending only on 7, m, A, 5i, 82, e(-) and C2, such that 

/^3(C2 + l)"^2^2 ^ ^^' sup e(-'S) + ^^2 ^ — T' /"2'^2 ^ /^5 inf -B'(s)- (4.26) 

2 0<S<K2 4 C<S<C. 

Owing to (j4.26|) and (j4.5p . one can get from the classical comparison principle that 

qc{^,ip) <qk{v,tp), i^p,ip)e0.k, A; = 1,2. 

Therefore, 

qdv, i^)<c,- /i5(-V?)^/'+\ (¥^, V') G [C-, 0] X [0, m,^]. (4.27) 

Owing to (I4.12|) and (14. 6j) . one also gets that 

A((7c(¥^,V')) <^(c)+/i4(C2 + l)(-(^)^/^ (</.,^) e[C-,0]x[0,mi„]. (4.28) 

Moreover, it follows from (|4.25|) and (|4.26p that 

Q,(C-,V')<c*-2(e(C_) + C-)(-C-)^/'+\ 0<V<mi„. (4.29) 

Set Co = sup(o'. We show that cq = c* by contradiction. Otherwise, cq < c*. Then, for any 
c*/3 < c < Co, c € (^ and gc satisfies (|4.27p - ()4.29p . where qc is the solution to the problem p.3ip - 
(1331, (USD. Set 

9co(V',^) = lim_gc(v',^), iv,'^) G [C-,0] X [0,min]. 

Then, it follows from ^^?Zh - KM\ that 

sup geo < c„ sup gco(C-,-) < c* - 2(e(C-) +Ci)(-C-)^/'+\ (4.30) 

(C-,0)x(0,min) (O.min) 

and it is not hard to show that ^cq is the solution to the problem (j3.3ip - ()3.34p . ()4.6p with c = cq. 
Hence cq € (o'. It follows from (|4.30p and the stability theory of uniformly elliptic problems that the 
problem (I3.3ip - (l3.34p . (|4.6p admits a unique solution with (14. 7j) and (|4.8p if cq < c < co + e with a 
sufficiently small positive number e. This contradicts that cq = suptf < c*. Therefore, to = (c*/3, c*). 
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For any c*/3 < c < c*, let Qc be the solution to the problem (|3.31|) - (|3.34p . (j4.6p with (j4.7p and 
(USD. Then, g^ satisfies (fiJT]) . (l4l2D . (11:271) and g29]). Set 

5(93, V') = lim_gc(v3,'0), (V'^V') e [C-,0] x [0,min]. 

One can show that q is the weak solution to the problem (J3.3ip - (j3.35p satisfying (j4.19p . (j4.20p and the 
second inequality in K2T\i . Moreover, q G C°°((C-,0) x (0,min)) n C^^-, 0) x [0,min)) n C([C-, 0] x 
[0,min]) thanks to the Schauder theory on uniformly elliptic equations. 

It remains to verify the first inequality in (j4.2ip . Similar to the proof of ()4.24p . there exists a 
Lipschitz continuous curve L : ip = u>{(p) (C- < V < 0) in [C-,0] x (0,min) satisfying 

g(V9,cD(v?)) > c. - Mi{-^)^/^+\ V? G [C_,0] 

with Ml = Mi(7, 771, A, 51,^2). Decompose (C-jO) x (0, min) into f^i and ^2 by the curve L, where 

ni = {{^,ij) G (C-,0) X (0,min) :^<a)(99)}, J^2 = {(^,V') e(C-,0) x (0, mi,,) : V > ^^M}. 

Prescribe a Dirichlet condition for q on L. Then the problem ()3.3ip - (l3.35p can be regarded as two 
problems on Qi and f22, respectively. Set 

q^{^,^l;) = c, - M2(-v.)^/2+i(i ^ ^2^ _ ^2)^ ^^^ ^) ^ ^^^ 
g^lv', V') = c* - M2(-(/?)^/2+^(l + m2„ - -02 ^ 2minV), (v', -0) S ^2. 

Owing to A > 2, (j4.3p and (|4.5p . there exist sufficiently small k,2 G (0,^2] and sufficiently large M2 
such that q is a subsolution to the problem of q on Q^ for A; = 1,2, where M2 depends on 7, m, A, 
5i and 62, while K2 also on e(-) and C2. Then, the first inequality in (j4.2ip follows from a comparison 
principle. Here, the comparison principle is not the classical one since p.3ip is degenerate at the sonic 
state and it can be proved in a similar way as for Proposition 3.1 in [20]. D 

4.4 Existence of subsonic-sonic flows 

We now choose 



Ci = ;U5C^/2+S C2 = /U6(c.^<5i + l)^^'^\ (4.31) 

Assume that /_ G C^{[l-,0]) satisfies ([3:22]) with -60 <l- <0, where 

c.y^fTi I ''V3(C2 + i)c.p(c2)y5fTT^ 'V2a.6^ r ^ ' 

which depends only on 7, m. A, 61, 62 and e(-). Assume that g G C^([0, /_(/_)]) satisfies ()3.23p and 
(13:21]) . Choose 

Co = /i4(C2 + l)(^c,^(^2 + lj sup|-— ^:c<s<c*-/X5(^-j |. (4.33) 

Set 

^ = {(^in,^-) G c°'i([o, /_(/_)]) X c([/_,o]) nci/2([z_,/_/2]) : 
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^in satisfies dH]) with (jOTD . while ^_ satisfies (02]) with KTE\\ and (I05D | 
with the norm 

\\{^in,^-)\\.y = max|||^in||Lcx>((oj_(i_))), ||^_||loo((;__o))|, i^in,^-) £ ^• 

For given (^in,^_) € ^, it follows from Lemma l4.ll and Proposition 14.11 that the problem (|3.31|) - 
([3:35]) admits a unique weak solution q € C°°((C-,0) x (0,min)) n C^([C-,0) x [0,min)) n C([C-,0] x 
[0,min]) satisfying the estimates (|4.19p - (|4.2ip . Set 

4n(y) = g(C-,*in(y)), 0<y </_(/_) 

and 

J_(a:) = g($_(a:),min), L < x < 0. 

It follows from KW - Km and (li3T|) - (l03|) that Jin € C°'i([0, /_(/_)]) satisfies dH}, and J_ G 
C([L,0]) n C ' {[l-,l-/2]) satisfies (14. 2p . Therefore, we can define a mapping J from ^ to itself as 
follows 

J((^in, ^_)) = (Jin, J_), (^in, ^-) £ ^. (4.34) 

Theorem 4.1 ^sstime ^/lai /_ G C3([/_,0]) satisfies IK22\i . and g G C3([0, /_(/_)]) satisfies ([5:23]) 
and (I323D. //-(5o < ^- < 0, then, the problem (f3:§T]) - (fSrsTj) arfmiis a wea^ solution q G C°°((C-,0) x 
(0,min)) nC^([C-,0) X [0,min]) nC([C-,0] X [0, rriin]) t«ii/i the following estimates 

||((^,V')| <M3(C2 + l)(-C-)^"'/'(-</^)'/^ (¥',V') G (C-,0) X (0,min), 

yi(g(^', V'')) - A{qi^",r))\ < i^a{C2 + iw - v"\'^' + w - r\), 

((^',V'),(9'",Oe(C-,0)x(0,min) 



and 



M6 



(_^)A/2+l < ^(^^^) < c, _^^(_^)A/2+l^ (^^^) ^ (^_^Q) ^ (0,min), 



where 6o and C2 are defined in (j4.32p and (j4.3ip . respectively, while fii{i = 3,4,5,6) are given in 
Lemma \4-3\ and Proposition \4-l\ Furthermore, 

min = m. (4.35) 

Proof. As mentioned above, the mapping J defined by (|4.34p is from ^ to itself. It follows from 
(|4.19p - (|4.2ip that J is compact. Therefore, the first part of the theorem follows from Proposition 14.11 
and the Schauder fixed point theorem provided that J is also continuous. 

Now let us show that J is continuous. Assume that {(cSj^" , =S_ )}J^q C =y satisfies 

lim ll^f^") - ^fn^||L-((oj^a_))) = 0, lim ll^i") - ^i°^|Uoo((;_,o)) = 0. (4.36) 
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Let mS^\ C-\ Mn\ ^!r\ ^-^ and xi"^ be defined by (l3:25]l - (f330D with ^i^ = ^f^^ and ^_ = ^i"^ 
for n = 0, 1, 2, • • • . Then, (fOHj) yields 

lim m'^^ = m^\ hm ci"^ = d°\ (4-37) 

Jim 11^;^") - ^L")||o,i;(o,/_a_)) = 0, Y^^ -^ y^f uniformly in {0,m^), (4.38) 

lim II^L"^ - $L°^ llo !■(/- 0) = 0, X^"^ — ^ Xi°^ uniformly in (d°\ 0). (4.39) 

It follows from the definition of J that 

where g(") € C7~((ci"\o) x (CmJ;;^)) n C7i([ci"\o) x [0, ?n["^)) n C7([d"\ 0] x [0,?n["^]) is the unique 
weak solution to the following problem 



aM(g(")) a2s(gW) 



+ .% ^ = 0' (V'-V^) e (d"\0) X (0,mJ::)), (4.40) 



(^H,^)_ ^B.(.)co^Ba.)^^ ^^^ ^.(0,^L")), (4.41) 



y=yir'(V') 



9A(^ („) ^ 9;„(j/)cos9in(j/) 

^(^,o) = o, ¥^e(d"\o), (4.42) 

aB(g(^ („) e^_(x)cos9„(x) („) 

gW(0,V')=c„ VG(0,m;;;)) (4.44) 

for 77, = 0, 1, 2, • • • . Set 

g(-)(^,7/;)=gW(d"V,rr7L")V^), ((^,7^) G (0, 1) x (0, 1), 77 = 0, 1, 2, ■ • • . 

It remains to verify that 

Jim (||gW(l, •) - gW(l, OIIl-CCo,!)) + lld"H-, 1) - 9(°)(-, 1)IIloo((o,i))) = 0, (4.45) 

which, together with (|i:37|) - (|iT39|) . imphes 

hm ||J((^J,"),^i"))) - J((^;°),^(°^))||,^ = 0. 

n— >-oo 

We prove (j4.45p by contradiction. Otherwise, there exist a positive number eq and a subsequence of 
{(?^"^}^xi denoted by itself for convenience, such that for each 77 = 1, 2, • • • , 

||gW(l, •) - ^W(l, ■)IIl°<>((o,i)) + lk"^"^(-, 1) - g^^^C-, 1)IIloc((o,i)) > £0. (4.46) 

Since g'""^ satisfies (j4.20p and (|4.2ip with r77,in = rrv^ and C_ = C- for each ri = 1, 2, • • • , there exists 
a subsequence of {(r'^'^^x-, denoted by itself again for convenience, such that 

lim llg^'^) - rilL-((o,i)x(o,i)) = 0, (4.47) 
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where 

q*{^,i^)=q*iC^'^0,m^i;), (0,^;) G (0,1) x (0,1) 

and q* G L°°((ci°\o) x (0,m[°^)) satisfies (fOOD and (fOT]) with nun = m|°^ and C- = C-^- Letting 
n -^ oo in (|4.40p - (|4.44p and using (|4.37p - (|4.39p and (|4.47p . one can get that g* solves the fohowing 
problem 

^.^^^0, (,,,,. (d»>0,x(„,™|.?, 



^(^,0)=0, ^G(d°\0), 



dB{q*) (o),_ 9^_(x)cose-(x) 



.(0) 



It follows from Lemma |4. II that 

q*{ip,i;)=q('\ip,i;), ((/p,V) G (d°\ 0) x (0,mj°)), 

which contradicts (j4.46p and (|4.47p . Hence (|4.45p holds. Finally, it follows from (|3.14p that the sonic 
curve of the flow lies at the throat of the nozzle and the velocity vector is along the normal direction 
at the sonic curve in the physical plane. Thus (I4.35P holds. D 

4.5 Uniqueness of subsonic-sonic flows 

We show the uniqueness of the subsonic-sonic flow in the physical plane. Note that the boundary 
condition ()3.14p means that the velocity vector of the flow is along the normal direction at the sonic 
curve. Thus, for the problem (|3.10p - (|3.15p . the solution satisfies both a Dirichlet and a Neumann 
boundary conditions at the sonic curve although there is a free parameter at the inlet. Hence the 
uniqueness theorem follows easily. 

Theorem 4.2 There exists at most one solution ip G C^'^(J7_) to the problem (I3.10p - (13.15p . 

Proof. Let (pi,ip2 G C^'^{Q-) be two solutions to the problem (I3.10|) - ()3.15p . Assume that 

fk{g{y),y) = Qn,/t, < y < /_(/_), A; = 1,2. 

For /c = 1,2, multiplying the equation with respect to ipk by {ipi — ^2) — (Cin,i — Cin,2) and then 
integrating over ri_ by parts, one can get that 

p(|Vv3fc(2;,y)P)Vv9fc(x,2/) • V(<y9i(x,y) - ip2{x,y))dxdy + (Cin,l - Cin,2)p(c^)c*/-(0) = 0. 

Hence 

{p{\Vipi{x,y)\^)\/ipi{x,y) - p{\\/ip2{x,y)\^)\/ip2{x,y)) ■ (Vipi{x,y) -Vip2{x,y))dxdy = 0, 
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Q_ 



which, together with (j3.14|) and (j3.15|) . leads to 



D 



4.6 Regularity of the subsonic-sonic flow 

The subsonic-sonic flow in Theorem 14.11 is only continuous near the sonic state. We investigate its 
regularity in this subsection. 



Theorem 4.3 Let q be the weak solution to the problem (I3.3ip - (j3.37p under the assumptions of The- 
orem l4n\ 

(i) If f- satisfies ()3.65p additionally, then q G C"'^([C-,0] x [0,m]) satisfies 

|i(v., V)| < M7(-9^)"/'+'/', ||^(</'>V')| < /Ur(-(^)"/2+S i^,ij) G (C-,0) X (0,m) (4.48) 

with fi7 = fi7{'y,m,X,6i,52,53). 

(ii) ///_ G C"'([/_,0]) satisfies (^3M\\ and (^Wi\\ additionally, then q G C^'H(C-,0] x [0,m]) 
satisfies 



d'^q 



dip'' 



ivii') 



< /^8, 



d'^q 



d(fd\(j 



(ip,^) <;.8(-v,)V4+l/2^ ^i^,^) <^8(-V^) 



dip'^ 
(</.,^)G(C-/2,0)x(0,m) 



\A/2+l 



(4.49) 



with fis = H8{l,'m,X, 61,62, 63,64). 

(iii) If g e C3'"([0, /_(/_)]) and /_ G C3'"([/_,0)) /or a number a G (0,1) additionally, then 
gGC2'"([C_,0)x[0,m]). 



Proof. Thanks to the Schauder theory on uniformly elliptic equations, one can get that ()4.48p 
holds in (C-,C-/2]x(0,m) in the case (i), g G C^'^((C-,0)x [0,m]) in the case (ii) and g G C2'°([C_,0)x 
[0,m]) in the case (iii). So it suffices to verify (jOSj) and (j09|) in (C-/2,0) x (0,m). 

For any < e < -C-/2 < 1/2, q G C2([-2e, -e/2] x (0,m)) n C^([-2e, -e/2] x [0,m]) satisfies 



!;("<-«£) -^(^<-«S)-- 



9^- 



((^,0) = 0, 



— {tf,m) = e{ip). 



(<^,V)G(-2e,-e/2) x(0,m), 

V^G(-2e,-e/2), 

¥^G(-2e,-e/2), 



where 



and 



a{ip, V') = ^'(g((/^, V')), Hip, V-) = S'(g((/J, V')), (</', V') e [-2e, -e/2] x [0, m] 



e{ip) 



e'_{x_{ip))cose-{x_{ip)) 

B'{q{ip,m))q{ip,m) 



ipe[-2e,-£/2]. 
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Shifting and rescaling the potential as 



ip = e 



-A/4-1/2 



i^ + e), ^e[-2e,-e/2] 



and setting 

q{^,i,) = 5(-e + eV4+i/2^^^) _ ^^^ (^^^) ^ j_^_V4+i/2^^-A/4+i/2/2] ^ fo,m], 

one can check that q G c2([-e-V4+i/2^^-A/4+i/2/2] x (0, ?n)) nC7^([-e-V4+i/2^g-A/4+i/2/2] x [0,m]) 
solves 



dq\ d 



d(p 

— (c^,m) = e((^), 



d(pJ dip 



dq 



dip 



^ e~'/'-M^,iP)^)+^ h{^,iP)^)=0, (^,V') e (_e-V4+i/2^, -A/4+1/2/2) x (0,m), 



^e(_e-A/4+i/.^£-A/4+i/./2), 



-A/4+1/2 g- A/4+1/2 
-A/4+1/2 -A/4+1/2 



^e(_£-^/4+V^^£-A/4+i/./2), 



where a and /i are defined on [— e •^/^+-^/^,e ^/^+'^/'^/2] x [0, m] by 

a((^, iP) = a{-e + ^^4+1/2^^ ^)^ /^(^^ ^) ^ ;,(_^ + eV4+i/2^^ ^)^ 

and 

e(^) = e(_e + eV4+l/2^)^ ^ ^ [_^-A/4+l/2^^-A/4+l/2/2]. 

It follows from the Holder continuity estimates for uniformly elliptic equations that there exists a 
number /3 G (0, 1) such that 

[9]/3;(-l/4,l/4)x(0,m) < ^fi ||g||ioo((_i/2,l/2)x(0,m)) < ^2e ' "^ , (4.50) 

where Mi = Mi{-f, m, A, (5i, ^2) for i = 1, 2. Then, ()i30]) and ([3:65]) yield 

^~ ~ ll^ll/3;{-l/4,l/4)x{0,m) + ll^ll/3;{-l/4,l/4)x{0,m) < ^Si l|e||/3;(-l/4,l/4) < ^-^3^ 

with M3 = M3(7,?Ti, A,5i, (52, (^s)- So one gets from the Schauder estimates on uniformly elliptic 
equations that 

ll^lll,/3;(-l/8,l/8)x(0,m) < ^4 ( ||g||L°°({-l/4,l/4)x(0,m)) + l|e||/3;(-l/4,l/4)) < M^E ' + , (4-51) 

where Mj = Mj(7, m, A, (5i, ^2, 5^) for i = 4, 5. It follows from ()4.5ip that 



with Mq = M6(7,r7i, A, 5i, (52,(53), which lead to (j4.48p owing to the arbitrariness of e G (0, — C_/2). 
Furthermore, we get from (|4.5ip . (|3.65p and (j3.67p that 



-A/2-1 II ~ I 



,/3;(-l/8,l/8)x(0,m) + ll^lll,;9;(-l/8,l/8)x(0,m) < ^7> P||l,/3;{-l/8,l/8) 



1/8) < MyeV^+i 



with My = M7(7, 771, A, 61,82,63,54). Using the Schauder estimates again yields 

ll^ll2,/3;{-l/16,l/16)x(0,m) < ^^8 ( ||9||L°°({-l/8,l/8)x(0,m)) + l|e||l,/3;(-l/8,l/8)) < MgE ' + , 

where Mj = Mj(7, m, A, 5i, (52, (53, (54) for i = 8,9. Therefore, 



d^q 



dif"^ 



i-e,^) 



< Mio, 



d'^q 



dipdip 



(-e,^) <Afioe^/^+i/2, |!|(-e,v.) <Afioe^/2+i, ^^€[0,^] 



5-0 



with Mio = Mio(7, "^i -^j 6i , 62,63,64), which yield ()4.49p owing to the arbitrariness of e G (0, —C,-/2). 
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5 Sonic-supersonic flows in the divergent part of the nozzle 

In this section, we will establish the well-posedness of sonic-supersonic flows in the divergent part of 
the nozzle. For convenience, A4. and 6k,+ will be abbreviated by A and 6k for fc = 1, 2, 3, 4, respectively. 
We will solve the problem (J3.47p ~ (l3.52p by a fixed point argument. Let ^^ be given in advance 
such that cS_|_ € C([0, /+]) satisfies 



c,<^+ix)<c= ^* + V^/(^ — ]1^ 0<x<l+. (5.1) 



Due to (j5.ip and the first formula in (j3.4ip . it holds that 



c* < ^+(x) < cJ6^ + 1, < X < 



■+) 



^' -(/.^ < GV(X+M) < ^^\ 0<(/.<C+. 



cA (52 + 1) A/2+1 ^^ - +^ ^^^^^ - c^ 

5.1 An iteration scheme 

The existence of solutions to the problem (J3.47p - (j3.52p will be proved by using the Schauder fixed 
point theorem. For given =S_|_ and Q in some suitable set, we solve the following problem 



Q^^ - {b{Q)Q^h = 0' 


((/J,V) e(o,C+) x(0,m). 


(5.2) 


Q{0,iP) = 0, 


tpe (0,m), 


(5.3) 


Q^{0,^l^)=0, 


ipe (0,m), 


(5.4) 


Q4^,o) = o, 


¥'e(o,c+), 


(5.5) 


^ ^ ^ GV(X+(^))cose+(X+((^)) 

0.1. in. m.) — ^ ■ — T- . 


10 <= (or. V 


C5.m 



^'''^' b{Q{ip,m)))Al\Q{ip,m)) ' 

If the problem (|5.2p - (|5.6p admits a unique solution Q, then we can define a mapping as 

J((^+,Q)) = (A;i(g($+(.),m)),Q(.,.)). 

A fixed point of the mapping J is a desired solution to the problem (|3.47p - (j3.52p . Note that the 
problem (|5.2p - (|5.6p is equivalent to 

W^ + b^^\Q)W^ = ^b~\Q)p{Q)WiW + Z), 

Z^ - b"''{Q)Z^ = -\h-\Q)p{Q)Z{W + Z), 

T^(0,V) = 0, 

Z(0,V') = 0, 

W^(99,0) + Z(v?,0)=0, 

W{(p, m) + Z{(p, m) = h{if), 

Q{0,iP)=0, 
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((^,V')G(0,C+)x(0,7n), 


(5.7) 


(^,V')G(0,C+)x(0,7n), 


(5.8) 


V'S (0,m), 




(5.9) 


V^e (0,m), 




(5.10) 


</'e(o,c+), 




(5.11) 


9'e(o,c+), 




(5.12) 


(^,^)G(0,C+)x(0,7n), 


(5.13) 


V^e (0,m), 




(5.14) 



where 

W{ip,tl;) = Q^{^,tl;)-b^/\Q{ip,^))Q^{ip,tP), (<^, V) G (0,C+) x (0,m), (5.15) 

Zi^,i,) = -Q^{ip,ij)-b'/\Q{ip,ij))Q^i^,i>), ((^,V')G(0,C+)x(0,m), (5.16) 

7, , 2ev(x+((^))cose+(x+(y)) 

Remark 5.1 For the system (j5.7p . (j5.8|) and (j5.13p . i/ie compatibility condition is of the form 

Q4^, V') = -^rV2(Q(^, ^)) (iy(^, ^) + z((^, V')) , ((^, V^) E (0, C+) X (0, m). (5.18) 

It turns out to be more convenient to solve the problem (J5.7p - (j5.14p instead of the problem (j5.2p ~ 
()5.6p . This win be done by solving first the problem ()5.7p - (l5.12p and then the problem (I5.13|) . ()5.14p . 
The problem (J5.7p - (j5.12p can be solved by using the contraction mapping theorem. For given (w, z) 
in some suitable set, we solve the following problem 



W^ + b^/\Q)W^ = lb-\Q)p{Q)W{W + 



(v9,V')G(0,C+)> 


< (0,m), 


(5.19) 


(v9,^)G(0,C+)> 


< (0,m), 


(5.20) 


^G (0,m), 




(5.21) 


^G (0,m), 




(5.22) 


V5G(0,C+), 




(5.23) 


^G(0,C+) 




(5.24) 



4 
Z^ - b^/\Q)Z^ = -^b-\Q)piQ)Z{w + Z), 

W{0,i;) = 0, 

Z{0,i;)=0, 

W{ip,0) + Z{ip,0) = 0, 

W{ip,m) + Z{(p,m) = h{(p), 

and obtain a solution {W,Z) = ^[{w,z)). Then, the unique fixed point of the contraction mapping 
^ is the unique solution to the problem (15.7p - (15.12p . 

In the following two subsections, we will solve the problem ()5.7p - (l5.14|) . and consequently, give a 
solution to the problem (l3:Tn) - (l332D . 

5.2 A hyperbolic system with singularity 

Fix =S+ G C([0,/+]) satisfying (|5.ip . Assume that Q G C°'^([0, C+] x [0, m]) is given and satisfies 

-^l2^^+^ < Q{^, V') < -^19^^+' {^, ^) G (0, C+) X (0, m), (5.25) 

-/32<^^+'<Q^(¥^,V')<-/3i/+\ |Q^(¥^,V)l</33¥^'^/'+\ ((^,^)G(0,C+)x(0,m), (5.26) 

where < /.ii < l-i'2, ^ < Pi ^ h-, Ps > 0. Here and what follows, fn, j3i, r, k and tj will denote 
generic positive constants. Owing to A > 2, (|3.4ip . (|5.25p and (|5.26p . there exists r G (0, 1] such that 
if < C+ < T", then 

^^^5A/4+l/2</^(^)<^^^5A/4+l/2^ /l'(^) > 0, (/.G(0,C+), (5.27) 

-/i5/32/+' < ^^(</', V') < 0, < Z((^, V^) < /i5/32<^^+\ (<^, ^) G (0, C+) X (0, m), (5.28) 

-//5/32(^"' < ^fe~HQ(v5, V'))P(Q(¥', V'))^(<^, ^) < 0, ((^, V) G (0, C+) X (0, m), (5.29) 

< \b-\Q{^, ^))piQ{v>, ^))Z{^, ^) < IJ^5l32^-\ i^, ^) G (0, C+) X (0, m), (5.30) 
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where the constants fii {i = 3, 4, 5) (/X3 < ^4) depend on 7, m, A, 5i and 82-, while r also on /3i and j3^. 
The problem (l5.7P " (l5.14p will be solved under the assumptions (I5.25p ~ (l5.30p in this subsection. 

Remark 5.2 The constants fii, /i2, /?i, P2 o.'i^d (3^, which will he determined in the proof of the 
existence theorem (Theorem \5.1\) . depend only on 7, m, \, 61 and 62- Therefore, the second condition 
in (|3.4ip . which is only used to guarantee the second estimate in (j5.27p . can be relaxed as noted in 
Remark\3^\ 



As mentioned in the end of the last subsection, we first consider the problem (I5.19p - (l5.24p . Note 
that the problem (j5.19p - (|5.24p can be decomposed into two problems: one is of homogeneous source 
terms and the other is of homogeneous boundary condtions, which are solved by the following three 
lemmas. For {w,z) G M^, we always use K'u;,^)! to denote max{|tt;|, |2;|} in this section. 

Lemma 5.1 The problem 

W^ + b''\Q)W^ = Q, (^,^)€(0,C+)x(0,m), (5.31) 

Z^-b^'^{Q)Z^ = Q, ((^.V') G(0,C+)x(0,m), (5.32) 

Ty(O,'0) = O, ^^£(0,^), (5.33) 

Z(0,^)=0, Ve(0,m), (5.34) 

Ty(^,0) + Z(99,0) = 0, V^G(0,C+), (5.35) 

Ty(^,m)+Z(99,m) = /i(v9), v9G(0,C+) (5.36) 



admits a unique weak solution (W,Z) G L°°((0, C+) x (0,r7i)) x L°°((0,C+) x (0, ?n)) satisfying 

lim 

¥3^0+ 



lim, |||(Ty,Z)(99,.)||L^((o^^))=0. (5.37) 



Furthermore, {W, Z) satisfies 

-^7(/j^+i < W{^, V) < -l^&^^^\ ^6^/^^+' < Z{^, ^) < ^i7^^+\ ((^, V) G (0, C+) X (0, m) 

and 

< W{ip, V) + Z{if, ^) < /i4¥^'^/"+^/^ (v?, ^) G (0, C+) X (0, m), 

where fiQ and ^j {fiQ < fj-j) depend only on 7, m, X, fii, fi2, Ms and /U4. 

Proof. First consider the uniqueness. Let {Wi,Zi) and (^^21-^2) be two weak solutions to the 
problem (|0T]l ~ (f06]l satisfying (fOT]) . Set 

(Ty,Z)(99,^) = (1^1 -1^2,^1 -^2)((^,V'), (V',^) e (0,C+) X (0,m). 

Then, (14/^, Z) satisfies 

^hm |l|(Ty,Z)(e,.)lL.((o,^))=0 

and is a weak solution to the following problem 

U^ + b^/^{Q)U^ = 0, (yp,^) € (e,C+) x (0,m), 
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V^ - b^'\Q)V^ = 0, {ip,^) € (e,C+) X (0,m), 

C/(e,V) = T^(e,V'), VG(0,m), 

y(e,^) = Z(e,^), V'e(0,m), 

[/(v?,0) + y(v9,0)=0, y^eCe.Cn-), 

[/(V9,m) + 1/(99, m)=0, <^€(e,C+) 
for any < e < ("+• Hence 

|(VF,Z)((^, V)| < |||(W^,^)(e, OlIL^ao,.^)), < e < y. < C+, V' G (0,m). 

Letting e ^ O"*" yields 

{W,Z){^,^)=Q, {^,ij) G (0,C+) X (0,m), 

i.e. 

{Wi, Zi){ip, V) = {W2, Z2)(V9, V), {V, ^) G (0, C+) X (0, m). 

We now turn to the existence. For any < e < C+, it follows from the classical theory for strictly 
hyperbolic systems that the problem 

W'^ + b'l\Q)WI, = 0, ((^.V') G (e,C+) X (0,m), 

Z^ - 6i/2(Q)^. = 0, {^, i;) G (e, C+) X (0, m), 

admits a unique weak solution {W'^,Z'^) G C^'^{[e,C+] x [0,m]) x C°'^([e,C+] x [0,m]). Moreover, it is 
clear from (|5.27p that 

W%^,i>) < 0, W^{ip,^P) < 0, W;{ip,iP) > 0, (y.,^) G [e,C+] x [0,m], 

Z^(^,7/;) > 0, Z^((^,^) > 0, Z;((^,V') > 0, {ip,iP) G [e,C+] X [0,m]. 

For any < ei < £2 < C+) it is not hard to show that 

W{ip,i;) < W'^ip,i;), Z''{^,^P) > Z'^{^,i^), (^,V) G [e2,C+] x [0,m]. 

Set 

(VF, Z)(^, V') = lim (Ty^ Z^)(^, V'), (</^, ^) G (0, C+) x (0, m). 



(5.38) 
(5.39) 
(5.40) 
(5.41) 
(5.42) 
(5.43) 



(5.44) 
(5.45) 



{ipk+i,tpk+i) 




the casg5that k is even 



(¥?fc+i,V'fc+ 




the case that k is odd 



Fix ((/9o,V'o) € (£)C+) ^ (0, m). We estimate W{^o,ipo) by the method of characteristics. Let 
ip = ^i{f) be the positive characteristic across {ipojipo), which approaches either the initial boundary 
{e} X [0, m] or the lower boundary (e, C+) x {0} at a point (ipijipi). If ipi > e, then there exists a negative 
characteristic ip = ^2(9^) from {(pi,ipi), which approaches either the initial boundary {e} x [0, m] or the 
upper boundary (e, C+) x {m} at a point {(p2,'4'2)- If 9^2 > &, then there exists a positive characteristic 
ip = ^3(1/3) from {(p2,ip2), which approaches either the initial boundary {e} x [0,m] or the lower 
boundary (e, C+) x {0} at a point (933, 7/^3). Since the system (|5.38p . ()5.39p is strictly hyperbolic, there 
exists a nonnegative integer k such that 



ti 



m > ^1 > ■ ■ ■ > fk > 'Pk+i = e 
0, I < j < k and j is odd, 



I jn, 1 < J < /c and j is even, 

f ^jC^) = ^^^^(<3(95, ^jiv))), fj<v< fj^i, 



< ipk+i < m, 

1 < J < ^ + 1 and j is odd, 

1 ^ J < ^ + 1 and J is even. 



Then, 



2<2j<k 



For 1 < j < A; + 1, it holds that 



V3-1 



b'/'{Q{s,^,{s)))ds = \^l;j.i-i;j 



Vj 



' € (0,m), if j = l, 
= 771, if 2 < j < /c, 

€ (0,m], if j = A; + 1. 



Thus, it follows from this and ()5.25p that 



¥'2 



A/4-1/2 



</'2j 



V2(j + 1) 



66 



/" 1^2(7 — 1) 

(^-//^-i/2(^20-i) - ^2j) > / s-'-'^-^'^ds > iljm, 2<2j<k + l, 

where /ig and jl-j are positive constants depending only on 7, A, //i and ^2- Hence 

2<2j<fc ^^"^ 2<2i<fe-l ^'^"^ -^^fe 

if A; > 2, 

^-^ ■' LLvm ^-^ J \J I J „ ^ / 

2<2j<k ^' 2<2j<k '^ ^ 

if fc > 0. 

These estimates lead to 

with fiQ and ^7 (/.fg < yt^?) depending only on 7, m, A, /Ui, fi2, /^s and /U4. Therefore, 

-^7(^^+1 < i^(^, ^) < -fi6ip^+\ (if, ^) G (0, C+) X (0, m). (5.46) 

Similar, one can get that 

//6¥'^+' < Z{^, V) < ^7<^^+\ (<^, ^) G (0, C+) X (0, m). (5.47) 

Moreover, it follows from ()5.44p and ()5.45p that 

< W'{ip, V) + Z'{^, ip) < h{ip), [if, ip) G [e, C+] X [0, m], 

which implies 

0<t^(<^,^) + Z(^,V)</U4<^'^/'+'/', (¥^,V)G(0,C+)x(0,m). (5.48) 

It is not hard to show that {W, Z) G L°°((0,C+) x (0,m)) x L°°((0,C+) x (0,m)) with (fSllGD - ffSlIHD is 
a weak solution to the problem ()5.3ip - (l5.36p satisfying (I5.37p . D 

We now turn to the problem with nonhomogeneous source terms and homogeneous boundary 
conditions, and consider first the strictly hyperbolic case and then the singular case. 

Lemma 5.2 Assume that < e < C+, Wq,Zq G L°°((0,m)), and w,z G L°°{{eX+) x (0,m)) satisfy 

Then, the problem 

W^ + b^/\Q)W^ = ^b-\Q)p{Q)W{W + z), (<^,V)G(e,C+)x(0,m), (5.49) 

Z^-b^/\Q)Z^ = -^b-\Q)p{Q)Z{w + Z), ((/p,V)G(e,C+)x(0,?n), (5.50) 
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W{e,i^) = Wo{ip), ^^€(0,™), (5.51) 

Z(e,V) = ^oW, VG(0,m), (5.52) 

Ty(<^,0) + Z(<^,0) = 0, 9JG(£,C+), (5.53) 

VI/((/j,m) + Z((^,m) =0, (^G(e,C+) (5.54) 

admits a unique solution {W,Z) € L°°((e,C+) x (0, m,)) x L°°((e,C+) x (0,m)). Moreover, (W, Z) 
satisfies 

\{W,Z){^,^)\ < -^_^^^^_^Kc^5A/4+i/2 ^ |||(|^,,^,)|||^^^^^^^^^^, (^,^) e (e,C+) X (0,m). 

Proof. According to the classical theory for strictly hyperbolic systems, the problem (j5.49p - 
(|5.54p admits a unique weak solution {W,Z) G L~((e,C+) x (0, m)) x L°°((e,C+) x (0, m)). 
We first estimate VF along a positive characteristic. Assume that 

S+ : ^ = ^+(93), i^<99<(^ (e<'^<'^<C+) 

is a positive characteristic of (|5.49p . i.e. 

[ < ^+(9?) < ?Tl, ip<ip<lf. 

On S+, VF satisfies 

^^(<^, $^(^)) _ h-\Qi^, ^+i^)))piQ{^, ^+{^)))Wi^, ^+{^))Wi^, *+((^)) 

=lb^\Qiv, ^+M))p(Q(</', ^+(v)))^(v', ^+M)2(</', ^+M), <^ < v^ < ^, 

i.e. 

1 /-^ 



VF((^, ^+(y.))exp{ --j b~\Q{s, ^+{s)))p{Q{s, ^+{s)))W{s, ^+{s))ds])' 
-z(v9,^+((^))(exp{ - i rri(Q(s,^+(s)))p(Q(s,^+(s)))W^(5,^+(5))d5})', 0<^<^. 



lip 
This leads to 

1 f^ 



1^(99, ^+iip))exp[ - i r fe-i(Q(s, ^+is)))piQ{s, ^+is)))W{s, M/+(s))ds})' 



K(^-"- ' 1/2 (exp{ - i I'' b-\Q{s, ^+{s)))p{Q{s, ^+{s)))W{s, ^+{s))ds})', ^ < ^ < ^. 



W{^, ^+(v?))exp{ - - / b~\Q{s, ^+{s)))p{Q{s, ^+{s)))W{s, ^+{s))ds] - W{0, ^+m 



Then, for any (f < (p < ip, it holds that 

1 r^ 

lifi 

<K 1^ t5A/4+l/2 (^^^p| _ 1 I* ^-1(Q(^^ ^i>_,is)))p{Q{s, ^+is)))Wis, ^+is))ds})'dt 
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(5.55) 



= K(^5A/4+l/2g^p| _ 1 |''6-l(Q(.,*+(s))MQ(s,M/+(s)))Ty(s,vI/+(s))ds} -K(^5A/4+l/2 

which yields 

(^^5A/4+i/2 _ |^(^^ xi,^(^))|)exp{ - \ J^ b-\Q{s, ^+{s)))piQis, ^+is)))W{s, ^+{s))ds} 

(5.51 
Similarly, we can estimate Z along a negative characteristic. Assume that 

S_ : -0 = ^-(v?), If < (f < (f {e < < ip < C+) 
is a negative characteristic of (|5.50p . i.e. 

r vlr'_(^) = -b'/\Q{^, M/_((^))), ^ < V^ < (^, 
[ < ^_((/7) < m, (^ < (/? < (^. 

On S_, Z satisfies 

(Z(^, M/„(99))exp{i I'' b-\Q{s, ^.{s)))p{Q{s, ^.is)))Z{s, ^.{s))ds})' 
= - w{if, M/_((p)) (expjl I" b'^Qis, ^^is)))p{Q{s, ^_(s)))Z(s, ^_(s))ds})', yi < ^ < ^, 
which yields 

(^Z{^, M/„((^))exp{i I" b~\Q{s, ^.{s)))p{Q{s, ^.{s)))Z{s, ^_(s))ds})' 
<K^'^/^+'/^(^exp{'^l''b-\Qis,^.{s)))piQis,^.{s)))Z{s,^.{s))ds]y, < if < ip. 
Therefore, for any (p < f < (p, it holds that 

iKip^^/''+'/'-\Zip,^i>.ipmexp{^J%~HQis,^.is)))p{Qis,^.{s)))Zi^^^ 

+ {lx + ^)^ l%'^/'~'/'exp{^ l\-HQ{s,^^{s)))piQis,^^is)))Zis,^^is))ds]dL 

(5. 



(5.56) 
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Next we estimate W by the method of characteristics. As shown in the proof of Lemma |5. 11 for 
any ((/3o,V'o) € (e,C+) ^ (0, m), there exists a nonnegative integer k such that 



V'i 



ipo > (pi > ■ ■ ■ > fk > fk+i = e 
0, I < j < k and j is odd, 



m, ^ < j <k and j is even, 



< ipk+i < in, 



^'jif) = b^^^Qif, ^jiv))), fj<V< fj^i, 



1 < j ^ ^ + 1 and j is odd, 
1 f^ i f^ ^ + 1 and j is even. 



I ^j ifj ) = "^i ' ^i (V'i-i ) = ^j-i . 

One can see the figures in the proof of Lemma |5. 11 Define a function r{s) (e < s < ip) as follows 
/ 1 



r{s) = < 



-b-\Qis,^,{s)))piQis,^,is)))Wis,^,is)), 



LPj < s < (pj~i, I < j < k + 1 and j is odd. 



-b~\Qis, ^,is)))piQis, ^,is)))Zis, ^,(s)). 



fj<s< (pj-1, 1 < J < A; + 1 and j is even. 
Then, for 1 < j < A; + 1, it follows from ([535]) and ([536]) that 

i^pfJi^'^' - |VF((/.,-i, V'.-ODexpj r~' r{s)ds} 
>AC(/.f /'+'/' - |T^(v^,-, Vi)| + (^A + ^)k r~' t5^/^-^/2gxp{ / r(s)ds}(it, if j is odd, 



(ft^y?^-/"^^ - |2'((/?j_i,?/'j-i)|)exp| / r(s)ds| 

which are equivalent to 

>(.^f ^-^^/^ - |M^(v^„ V',)! - (^A + i).£^ t^V4-i/2,,p| _ j^^ r{s)ds}dt) 

• exp-{ — / r{s)ds\, ifjisodd, (5.57) 



J IS even. 



>dt 



ifij 



Kip^j^l'"^ - \Z{ipj_i,tpj_i)\ 



5, 1 

4A + -.K 



fj-i 



t5V4-l/2g^p| 



1/5J-1 



■(s)(is} 



dt 
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• exp< — / r{s)ds[, if j is even. (5.58) 

If k is odd, one can get from ([5371) . ([538]) . ([53T]) - (l53i]l that 

>{^^IT'^' - \W{^,^^,i^,^^)\ - (^A + ^)-/"'" t^^/^-^/^exp{ - I""' r{s)ds}dt) 
■ exp< — / r{s)ds > 

=(Ke5V4+i/2 _ \Wo{iPk+i)\)ew{ - P r{s)ds] 

>^^5A/4+l/2^^p| _ p ^(^)^^| _ |^q(^,^^)|. (5.59) 

Similarly, if k is even, one can get that 

n^f"^''' - \W{^,M)\ - (^A + i)^^" t^V4-i/2^^p| _ ^"° ^(^)^^l^, 

>^^5A/4+l/2g^p| _ y ^(^)^^| _ |Zo(V'fc+l)|. (5.60) 

Then, it follows from (l539|) and (IS^Ol) that 

|W((^o, V'o)! <^^o'/'+'/' + IlKW^o, ^o)l|Lo.((o,^)) - ^£^'/'+^/'exp{ - j r{s)ds] 

<K(/?0 + II 1(14/0,-^0 j III ioo((om)) -'^V^o ^ 

_(^.^).^-,v.-v.(^)-,< 

^K-,n^^/^+^/^ -U lllCT/f/ 7 "llll ^,^-M5/32^5A/4+l/2+At5/32 

<K(^o + II 1(14/0,^0 j I ||ioo((o,m)) "'^'^0 ^ 

^A + 2 , ^5A/4+l/2 ,^-/^5/32^5A/4+l/2+M5/32^ 

"5A + 2 + 4/^5/32 ^^° -V'o ^ ) 

^ 4/X5/32 5A/4+1/2 , |||/... .7 s,||| 

- 5A + 2 + 4/.5/32 "'"" + llK^o,^o)||L.((o,„)). 

The estimate of Z can be proved similarly. D 

Lemma 5.3 Assume that {w,z) € L°°((0,C+) x (0, ?7i)) satisfies 

\{w,z){ip,tl^)\ < V^/4+'/', (¥',V') e (0,C+) X (0,m). (5.61) 
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(¥P,^)€(0,C+)x(0,m), 


(5.62) 


(^,V')€(0,C+)x(0,m), 


(5.63) 


V'€ (0,m), 


(5.64) 


V'G (0,?n), 


(5.65) 


<^e(o,C+), 


(5.66) 


v'e(o,C+) 


(5.67) 



Then, the problem 

W^ + b'^\Q)W^ = ^b-\Q)p{Q)W{W + z), 
Z^ - b^/\Q)Z^ = -^b-\Q)p{Q)Z{w + Z), 

ty(o,v) = o, 
z(o,v) = o, 

iy((^,0) + Z((^,0) = 0, 
W(c^,m) + Z{(p,m) = 0, 

admits a unique solution (W, Z) G L°°((0,C+) x (0, jn)) x L°°((0,C+) x (0,m)) satisfying 

^limJ||(T^,Z)(v.,-)l|L.o((o,™))=0. (5.68) 

Moreover, (W, Z) satisfies 

Proof. We start with the uniqueness. Let {Wi,Zi) and {W2,Z2) be two weak solutions to the 
problem ([KIMjl - dHIHT]) satisfying dSSH]). Set 

(W^,Z)((^,V') = (1^1 -1^2,^1 -^2)(9',V'), (v^,^) e (0,C+) X (0,m). 

Then, {W, Z) satisfies 

^limJ||(Ty,Z)(^,)||L.,o,„,,=0 

and is a weak solution to the following problem 

U^ + b'/\Q)U^ = \b~\Q)p{Q)WU, (¥P,V) G (e,C+) x (0,m), 

V^^ -^'/'(Q)^V' = -|?^-H0)P(Q)^^, ('/',V') G (e,C+) X (0,m), 

C/(0,V) = T^(e,V'), Ve(0,m), 

y(0,^) = Z(e,V), ^G(0,m), 

C/(v;,0) + y(<^,0) = 0, </'G(e,C+), 

f/(v9, m) + 1/(99, ?n) = 0, 99G(e,C+) 

for any < e < C+- It follows from Lemma 15.21 that 

|(VF,Z)((/p,V)| < |||(W^,Z)(e,.)lL^((o,^)), < e < v? < C+, V' G (0,m). 

Letting e ^ Q'^ yields 

{W,Z){^,'^) = (0,0), {^,,1:) G (0,C+) X (0,m), 

i.e. 

(Wi, Zi)((^, V) = (VF2, Z2)((/p, V), (y^, ^) G (0, C+) X (0, m). 
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We now turn to the existence. For any < e < C+7 it follows from Lemma 15.21 that the following 
problem 

W^ + b'/\Q)W^ = ^b-\Q)p{Q)WiW' + z), (y.,^) G {e,C+) x (0,m), 

Zj-6i/2(Q)z; = -^b~\Q)piQ)Z{w + Z'), {ip,iP) G (e,C+) x (0,m), 

W%e,ij)=0, ^G(0,m), 

Z%e,ij)=0, ^|J€{0,m), 

W^ip,m) + Z%y,,m) = 0, ^ ^ {e,C+) 

admits a unique solution {W^, Z^) G L°°((e, C+) x (0, m)) xL°°((e,C4.) x (0,m)); furthermore, {W'^,Z'^) 
satisfies 

l(P7^z^)(^,v)l < ^xttS^;:^''^''^''''^'' ^^'"^^ "^ ^''^+^ "" ^°'"'^- ^^-^^^ 

Define 

{W, Z'){^, ^) = (0, 0), {^, ^) G (0, e) X (0, m). 

Then, there exists a decreasing sequence {en\^=i C (0, C+), which tends to as n — )■ cx), and W,Zq 
L°°((0,C+) X (0,m)) such that 

{W'",Z^'') -^ {W,Z) weakly * in L°°((0,C+) x (0,m)) x L°°((0,C+) x (0,m)) as n ^ cx). 



Moreover, it follows from (15. 691) that 



It is not hard to show that {W, Z) with (|5.70p is a weak solution to the problem (j5.62p - ()5.67p satisfying 
(15:68]) . D 

Remark 5.3 In Lemma lS.cH if (j5.6ip is replaced by 

\w{ip,ip)\<K^p°', \z{ip,^p)\<K^p"', (c/?,^) G (0,C+) X (0,m) 

with a positive constant a. then the solution (W, Z) to the problem (J5.62p - (j5.67p satisfies 

\{W,Z){ip,^P)\ < f ^% Kip^, ((/.,V) G (0,C+) X (0,m). 

To solve the problem (|5.19p - (|5.24p . we introduce a norm 

\\{w,z)y=snp{ip-'^/^-'/^\{w,z){ip,n : {^,^P) G (0,C+) x (0,m)} 
and define 

<^={{w,z) gL°°((0,C+) x(0,m)) xL°^((0,C+)x(0,m)): ||(u; - P^*, z - Z*)|l.^ < +oo}, 
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where {W*,Z*) is the unique weak solution to the problem (j5.3ip - ()5.36|) satisfying ()5.37|) . It follows 
from Lemma l5. II that 



-/U7(^^+' < W*iip, ^) < -^69^^+', ^69^^+' < Z*{ip, i^) < ^Lr^^+\ (<^, ^) G (0, C+) X (0, m) 

(5.71) 

and 

< W*{^, V) + Z*{^, ^) < fi4^^^/^^'/\ ((^, V) € (0, C+) X (0, m). (5.72) 

Then we have 

Proposition 5.1 For given {w,z) G ^, the problem (j5.19p - (|5.24p admits a unique weak solution 
{W,Z) G ^. Furthermore, if \\{w -W\z- Z*)\\gg < ^8/32 additionally, then \\{W -W*, Z - Z*)y < 
IJ.8/32, where fis = 4:fi4fi5/{5X + 2). 

Proof. Since the uniqueness follows from a similar proof in Lemma 15.31 it suffices to consider 
the existence. Consider the following problem 

W^ + h^'\Q)n = \^~^ {Q)P{Q)W{W + {W* + Z*) + {z- Z*)) , {^, V) G (0, C+) X (0, m), 

(5.73) 

JS:-^ - h^'\Q)^^ = -h}-\Q)p{Q)Z{{w - W*) + {W* + Z*) + ^), (<^, V) G (0,C+) X (0,m), 

(5.74) 
^(0,'i/^)=0, V€(0,m), (5.75) 

^(0,^)=0, V'e(0,m), (5.76) 

^(y., 0) + 2f{ip, 0) = 0, ¥. G (0, Ch-), (5.77) 

W{ip,m) + ^{ip,m)=Q, (/7G(0,C+). (5.78) 

Thanks to Lemma[5?3l the problem (j5.73p - (|5.78p admits a unique weak solution {W ., 3f) G L°°((0, C+) x 
(0,m)) X L°°((0,C+) X (0,m)) satisfying 

K^' ^)(V'> V')! < ,x^tf! . (/^4 + IKu; - M^*, ^ - Z*)||,^)^5V4+i/2^ 
5A + 2 + 4^5/?2 

(y.,V')G(0,C+)x(0,m). (5.79) 

Therefore, (:F + W*, ^ + Z*) G ^ and it is not hard to verify that (^^ + V^*, ^ + Z*) is a weak 
solution to the problem (J5.19p - (j5.24p . Finally, the rest of the proposition follows from (j5.79p directly. 

D 

We are now ready to solve the problem (J5.7p - (j5.14p . 

Proposition 5.2 For given =S+ G C([0,/+]) with ()5.ip and given Q G C°'^([0, C+] x [0, m]) with 
(f5r25]) - (f5l30]) . the problem (fSTT]) - ([5TTi]) admits a unique weak solution {W,Z,Q) with {W, Z) G ^. 
Furthermore, the following estimates hold: 

-firV^-^' - /i8/32^'^/'+'/' < QAV, ^) < -^W^^^ + /X8/32y^'"/'+'/^ (V^, ^) G (0, C+) X (0, m), 

(5.80) 
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|Q^((/p,V)| </i9(/52 + l)'/53V2+i^ (^,V')e(0,C+)x(0,m) (5.81) 

and 

-^2^+2 - m^9(/32 + l)<^=^^/2+l < g(^^^) < -ai(^^+2 ^ ^^g(^^ ^ l)(p3V2+l^ 

(<^,V)G(0,C+)x(0,m), (5.82) 

where fig depends only on j, X, fj,i, ^2, ^^4 and fis, while ai and 02 (o"i < (T2) depend only on 7, m, \, 
5i and 82 ■ 

Proof. Let {w,z) a ^ he given. It follows from Proposition 15.11 that the problem (j5.19P " (|5.24p 
admits a unique weak solution (W, Z) € i^. Define the mapping 

^■.^ — >^, {w,z)^^{W,Z). 

For given (wi, zi), {w2, Z2) € SS, it follows from the definition of =^ that 

^{{wi,zi)){^,^) - ^{{w2,Z2)){^,ij) = {U,V){^,^), (99,^) G (0,C+) X (0,m), 

where {U,V) G L°°((0,C+) x (0, ?n)) x L°°((0, C+) x (0,m)) solves the following problem 

U^ + h^'\Q)U^ = h)-\Q)p{Q)W{U + {zi-Z2)), (95, V') G (0,C+) X (0,m), 

y^ - 6i/2(g)y^ = -'^-h-\Q)p{Q)Z{{wi - W2) + V), {^, 4,) G (0, C+) x (0, m), 

C/(0,^)=0, Ve(o,m), 

y(o,V) = o, Ve(o,m), 

[/(<^,o) + y(v;,o) = o, <^e(o,c+), 

[/(V9, m) + y((/p, m) = 0, V9G(0,C+) 



satisfying 



^limJ||(f/,y)(v.,.)l|L.((o,.))=0. 



Moreover, Lemma 15.31 leads to that {U,V) satisfies 

\{U,V){^,n < ^^—^^^^-^IK^i -«;2,Zi -z2)|LV^/^+i/2, (^,V^) G (0,C+) X (0,m). 



Therefore, 



\\^{{wi,Zi)) - ^{{w2,Z2))\\:^:g < ^)^j^2 + la B ^^^^^ -W2,Zl- Z2)l^, 



which shows that J7 is a contraction mapping on ^. Thus, the problem (|5.7|) " (|5.12|) admits a unique 
weak solution (W,Z) G ^; furthermore, one gets from Proposition 15.11 that 

\\{W-W\Z-Z*)\y<f,s/32. (5.83) 

Then, we can get the weak solution Q G C'^'^([0,C+] x [0, ?n]) to the problem (j5.13p . (j5.14p . And this 
{W, Z, Q) is the unique weak solution to the problem (|5.7P " (|5.14p with (W, Z) G ^. It follows from 
dSSSD, (I5TT]1 and ([5^ that 

-fi7^^+' - ^.8/329^'^/"+'/' < Q^{^,i^) < -/i69'^+' +/^8/32'/5'^/'+'/^ {^,i^) G (0,C+) X (0,m), 
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-^/i4V'''/'+^/' -/i8/32^''/'+^/' < by'{Q)Q4^,^l^) < /.8/32^^^/^+l/^ (c^,V) G (0,C+) x (0,m), 

which yield (|5.80p and (|5.8ip immediately. 

Finally, we prove (|5.82p . Owing to (|5.7p , (jS.Sp , (|5.13p and (|5.18p , it holds in the sense of distri- 
bution that 

= (6(Q(V7, V'))Q^((/p, V))^, (V5, V') G (0, C+) X (0, m). 

It follows from this formula, (jS.lip and (|5.12p that 

<f Tnr /w/ eV(X+M)cos9+(X+M) 

d^^ Jq A^\Q{ip,m)) 

which, together with (|5.80p and (|5.14p . yields 

"nr /w/ r /•* eV(x+(5))cose+(x+(5)) ^ ^^ 

Q{ip,ip)dip = - / ,_!,-, TT dsdt, 0<99<C+- 

io io ^+ (Q(s,m)) 

Then, there exist ai and (T2 (o"i < (T2) depending only on 7, m. A, Ji and 82 such that 

1 f™ 

-^2(^^+2 < - / Q((/J, il^)d'4: < -aiip^+^, 0<ip<C+, 
"^ Jo 

which, together with (fOT]) . leads to ^Mil- □ 

5.3 Existence of sonic-supersonic flows 

Based on the well-posedness of the linearized problem (j5.7p - (|5.14p . which is equivalent to (I5.2p - (l5.6p . 
we are going to prove the existence of weak solutions to the nonlinear problem (J3.47p - (j3.52p . 

Assume that /+ G C'^([0, Z+]) satisfies p.4ip with < /+ < 5o, where 60 is a positive constant to 
be determined below. Let J2+ € C([0, /+]) and ^ € C°'^([0, 1] x [0, m]) be given such that (fET]) holds 
and 

-2a2(c^<5i + l/+)^^V+' < ^(0,V') < -^cTi(cJ+)^+2(/>^+2, (c/,,V;) G (0,1) x (0,m), (5.84) 
-/32(c^<5i + l/+)^^%^+^ < ^^(</>,V) < -/3i(cJ+)^+2</,^+i, (,^,^) G (0,1) X (0,m), (5.85) 

\^^{(t>,i^)\ < /33(c^<^i + l^+)'^^'^V'^/'+S (0,^^) G (0,1) X (0,m), (5.86) 

where ai and cr2 are defined in Proposition 15.21 which depend only on 7, m, A, 61 and ^2) while 
/3j (i = 1, 2, 3) (/3i < /32) will be determined below. Define 



Q{ip,i>) = ^(^.V'), {^,4^) G (0,C+) X (0,m), 
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where C+ is given in I^M\\ . Then, Q G C°'i([0,C+] x [0,m]) satisfies 



-2a2(Hvl±l)'+%A+2 <Q(^^^) < -U(^^Y~'%'^', {^,i.) e (0,C-i-) X (0,m), 
_^ / cV^fTT xA+2^,^, < Q^{ip,^P) < -/3i( ^_^ )^^%^+\ ((p,^) e (0,C+) X (0,m), 

|g^(v^,V')l < U^^^^±iy^^'^'^^^/^+\ (c^,V^) e (o,c+) X (o,m). 

It follows from these estimates and (13.4ip that there exists ri G (0, 1] such that if < C+ ^ '^ii then 
^^^5A/4+l/2<;^(^)<^^^5A/4+l/2^ /i'M > 0, ^G(0,C+), 

-/^5/32/+' < M^lv?, V') < 0, < Zi^, ^) < /i5/32V5^+\ (V', 4>) e (0, C+) X (0, m), 

-//5/32¥'-' < ^fe"'(g(v5, ^))p{Q{^, i^))W{ip, ^) < 0, {if, V) e (0, C+) X (0, m), 

< \b~\Qi^, ^))p{Q{ip, ^))Z{^, ij) < ii^h^~\ {^, ij) G (0, C+) X (0, m), 

where h, W and Z are defined by (|5.17p . (jS.lSp and (|5.16p . respectively, /ij (i = 3,4,5) (/^s < ^4) 
depend on 7, m, A, 5i and (52, while ri also on /3i and /Ss. According to Proposition 15.21 the problem 
(f5Tl) ~ (|5.14p admits a unique weak solution {W, Z, Q) with (W, Z) G =^ if < C+ < n; furthermore, Q 
satisfies ()5.80p ~ ()5.82p . That is to say, the problem (|5.2p - (j5.6p admits a unique weak solution Q with 
{W, Z) € ^ if < C+ < Ti, and Q satisfies 

-^.7/+! - ^.8/329^'^/'+'/' < Q^((^,^) < -/i6/+' +/i8/32'/5'^/'+'/', (V^,^) G (0,C+) X (0,m), 

(5.87) 

|Q^(</.,^)|</i9(/32 + l)(p'^/'+\ (<^,^)G(0,C+)x(0,m) (5.88) 

and 

(v?,^)G(0,C+)x(0,m), (5.89) 

where ^i {i = 6, 7, 8, 9) (/ig < /^t) depend only on 7, tti, A, 81 and 52- Choose 

/3l = -/i6, /32 = 2;U7, /33 = /i9(/32 + 1) (5.90) 

and 

. // /i6 \ 4/(^-2) . ^, s2/{A-2) /-^(c)Nl/(A+2)-. 

which all depend only on 7, m. A, (5i and 82- Then, if < C+ < ^"2) it follows from (j5.87p - ()5.9ip that 

-/32/+' < QA^, ^) < -/3i/+\ iv, ^) G (0, C+) X (0, m), (5.92) 

|Q^(V?,V')| < PzV^^'^^\ {^,i^) e (0,C+) X (0,m), (5.93) 
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-2a2V?^+2 < g(^^ ^) < -^ai/+2^ (yp, V) G (0, C+) x (0, m) (5.94) 

and 

c* < A;HQ((^, V)) < c, (v?, ^) G (0, C+) X (0, m). (5.95) 

Set 

J+(3;) = A-^(Q($+(x),m)), 0<x<;+ 

and 

J((/>, V) = Q(C+'/', V'), (</>, V') G (0, 1) X (0, m). 

Due to (jSlMjl - lfOSj) . one can get that J+ G C([0,/+]) satisfies 

c* < J+(x) < c, < X < /+, (5.96) 

and J G C°'H[0, 1] x [0,m]) satisfies 

-2<T2(c^<52 + l/+)^^V+' < ^('/',V') < -^cTi(cJ+)^+V^+', ('/',V') G (0,1) X (0,m), (5.97) 

-/32(cY^<52 + l/+)^"^%^+i < J4</.,V) < -/3i(cJ+)^+2,^^+i, (,^,^) e (0,1) X (0,m), (5.98) 

M<P,ip)\ < /33(cy^<5i + H+)'^^'^V'^/'+\ ('A,V') G (0,1) X (0,m). (5.99) 
Now we choose 

^0 = _ ,,.) , min{ri, r2}, (5.100) 

c^/6^ + 1 

which depends only on 7, m, A, (5i and 62- Due to (j5.96p - (|5.99p . we get that if < Z+ < (5o, then 
< C+ < min{Ti,r2}, J+ G C([0,/+]) satisfies dSl]) and J G C°'H[0, 1] x [0,m]) satisfies (IS31D-(IS3S1) 
with (fODD . 

For < /+ < (5o, define 

^ = {(^+,^) G C([0,/+]) X CO'i([0,l] X [0,m]) : 

^+ satisfies (^, while ^ satisfies ([Qijl - ffsrsel) with (l5:90]l | 

with the norm 

||(,^+,^)||.r =max|||=2+||icx=((o,/_^)),||^||icx,((o,i)x(o,m))|, (^+,=S) G =y. 
Owing to the discussion above, we can define a mapping J from =5^ to itself as follows 

J((^+,^)) = (J+, J), (^+,^)G^. (5.101) 

Now the existence of sonic-supersonic fiows can be stated as 
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Theorem 5.1 Assume that /+ G C^([0,/4-]) satisfies (j3.41|) . // < /+ < 60, then, the problem 
(J3.47p - (l3.52p admits at least one weak solution Q satisfying 

(</.,V')G(0,C+)x(0,m), (5.102) 



wheredQ, ai{i = 1,2) and(3i{i = 1,2,3) are defined in (jS.lOOp . Proposition \5.2\ and (j5.90p . respectively, 
which all depend only on 7, TTi, A, 5i anc? (52- 



Proof. As mentioned above, the mapping J defined by (jS.lOip is from J^ to itself. It follows 
from (j5.85p . (|5.86p and the embedding theorem that J is compact. Therefore, Theorem 15.11 follows 
from (|5.92p - (|5.94p and the Schauder fixed point theorem provided that J is also continuous. 

It remains to show that J is continuous. Assume that {(=S_^ , =S'-"^)}^q C =5^ satisfies 



hm II^J^) - ^f ||loo((o,,^)) = 0, hm ||^(") - ^(°)||lo.((o,i)x(o,™,)) = 0. 



(5.103) 



(n) „(n) 



.(n) 



(n) 



Let $Y ' ^+ and C+' be defined by (IQijl - ffOH]) with ^+ = ^^"^ for ?z = 0, 1, 2, • • ■ and set 

Q(")((^,V) = ^(")(^,V'), (99,^)G(0,C?^)x(0,7n), n = 0,l,2,.... 



Then, (I5lil3]l leads to 



.(n) 



lim \\^ 

n— >oo 



C«) ^(0) 



$ 



+ II0,1;{0,/+) 



0, limCi"^ = cf, 

n— >oo 
j{0) 



Q*^ '' uniformly in (0, C+ ) x (0, m 



X'^'' -^ X|°^ uniformly in (0, C+^), Q^") 
It follows from the definition of J that 

J((^J'\^("))) = (JJ'\J(")), n = 0,l,2,--- , 
where 



(0)^ 



Ji"^(rc) = A;^(Q(")($i'^a;),m)), < x < /+, n = 0,l,2,. 



and 



J(")(,/),^)=Q(")(Ci"V,^), (</.,^)G(0,l)x(0,m), n = 0,l,2,. 
with Q^"^ being the unique weak solution to the following problem 



)(n 



Ql^^-{b{QnQ\;% = o, 






(n)N^(n)^ 



(0,^) = 0, 

(o,V') = o, 

((^,0)=0, 
(c^, m) 



,(n) 



C«), 



eV(x|"^(9.))cose+(x^"^(^)) 

6(QW((^,m)))V(Q(")(^,m))' 



¥'G(0,C 



in). 



(5.104) 
(5.105) 



((^,V')G(0,cf^)x(0,m), 


(5.106) 


V'G (0,m), 


(5.107) 


V' € (0,m), 


(5.108) 


^e(o,d"^ 


(5.109) 



(5.110) 
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for n = 0, 1, 2, • • • . So, to show show the continuity of J, it suffices to prove 

Hm II J(") - J(°)||Loc((o,i)x(o,m)) = 0, (5.ni) 

since this hniit, together with (j5.104p . also imphes 

hm ||jP-Jf ||ioo((o,,^))=0. 

We win prove ()5.1f II) by contradiction. Otherwise, there exists a positive number eq and a subsequence 
of {^^'^)}J^]^, denoted by itself for convenience, such that 

II J^'^) - J(0) ||L.o((o,i)x(o,m)) > eo, n = 1, 2, • • • . (5.112) 

Since S-^"^^ satisfies (|5.84p - (j5.86p for each n = 1, 2, • • • , there exist a subsequence of {=^(")}^;^, denoted 
by itself again for convenience, and a function S.* with ()5.84p - ()5.86p such that 

J(n) _^ ^* in L°°((0, 1) X (0, m)) as n ^ oo, (5.113) 

Jj"^ ^; and Jj"^ ^ =S; weakly * in L°°((0, 1) x (0, m)) as n ^ oo. (5.114) 

Set 

Q*(^,V') = ^*(^,V), (</',^) G (0,Cf ) X (0,m). 

Letting n -^ oo in (j5.106p ~ (l5.110D and using (j5.104p . (J5.105D . (15.1131) and (J5.114I) . one can get that 
Q* solves the following problem 

Q;^ - (fe(Q(°))Q;)^ = 0, ((^,V^) G (0,cf ) X (0,m), 

Q*(0,V) = 0, ^G(0,m), 

Q;(0,V')=0, ^^£(0,^), 

<5;(<^,o) = o, (/^€(o,cf), 

0*r ^ 9V(xf(v.))cose+(xf(v.)) (0), 

^V'^^'"^) = -Wo)(,,^)))^-(g(o)(,,^))' ^^(°'^- )• 

It follows from the uniqueness in Proposition 15.21 that 

Q*((^,V) = Q(°)(¥',V'), (<^,V') G (0,Cf ) X (0,m), 
which contradicts (|5.112l) and (|5.113p . Hence (|5.11ip holds. D 

5.4 Smooth sonic-supersonic flows 

In this subsection, we establish the existence of C^'^ sonic-supersonic flows under the assumption that 
/+ € (^^([0, /+]) satisfies ()3.4ip and ()3.7ip . The method is still a fixed point argument and the method 
of characteristics. To this end, one should establish a priori C^'^ estimates for the linearized problem 
(I5TD- (I5TID . 
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Lemma 5.4 Assume that {W,Z) is the unique weak solution to the problem. (j5.3ip - (|5.36p in Lemma 
\5.1[ Then, 



\{W^,Z^){ip,i;)\ < ^lo/, (¥^,V) e (0,C+) X (0,m) 
with /_iio depending only on 7, m, A, 61, 62, ^3, /-ti, /U2 and (32- 

Proof. Set 

(C/,y)(v;,V) = {W^,Z^){ip,4^), ((^,V) e (0,C+) X (0,m). 

Then, ([/, 1/) is the unique weak solution to the fohowing problem 

U^ + h^'\Q)U^ = \h~\Q)p{Q)Q^U, ((^,V) e (0,C+) X (0,m), 

V^-h'l\Q)V^ = h~\Q)p{Q)Q^V, ((^,V) G (0,C+) X (0,m), 

t/(o,V) = o, Ve(o,m), 

y(o,V') = o, Ve(o,m), 

t/((^,o) + y((^,o) = o, <^G(o,C+), 

U{ip,m) + V{ip,m) = h'{ip), V?G(0,C+). 
Decompose (L'^, F) as 

{U,V){^,^) = {U\V^){ip,i^) + {U\V^){ip,ij), (<^,^) € (0,C+) X (0,m), 
where {U^^V^) and (J7^,y^) solve the following problems 

Ul + 6i/2(g)[/^ = 0, {ip, V) € (0, C+) X (0, m), 

yj -6i/2(Q)V;i = 0, (<^,V) € (0,C+) X (0,m), 

yi(o,^) = o, ^^£(0,7^), 

^n</',o) + yn<^,o) = o, ¥'e(o,C+), 



and 



Ul + b'l\Q)Ul = h,-\Q)p{Q)Q^{U^ + U^), ((^,V) G (0,C+) X (0,m), 

y^2 _ 6V2(g)y^2 ^ h~\Q)p{Q)Q^{y^ + y2), (<^, V') G (0,C+) X (0,m), 

C/2(0,V)=0, V€(0,m), 

^'(0,V)=0, V€(0,m), 

c/2(<^,0) + 1/2(^,0) = 0, v'e(o,C+), 

respectively. It follows from the same proof as for Lemma 15.11 that 

\{U\V^){^M < ^/^lO^^ (v',^) e (0,C+) X (0,m), 
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where ;Uio depends only on 7, m, A, 5i, 62, 5^, fii, ^2 and /32- Then, Remark 15.31 gives 

\{u^v^){ip,2P)\ < ^m^\ {ip,ip) e (o,c+) X (o,m). 

Therefore, 

\{U,V){ip,i^)\ < fiio^\ {y^,tP) G (0,C+) X (0,m). 



D 



Now, we assume additionally that 

\Q^Av,i^)\ < f3Av\ \QM^,n < /^4¥''^/"+'/', {^,i^) € (0,C+) X (0,m) (5.115) 

with /34 > 1 to be determined later. Then, 

\iW^,Z^)i^,ij)\ < fiuf34ip\ (¥^,V) e (0,C+) X (0,m), 

where //n depends only on 7, tti. A, 5i, 62, /^i, ^2 and f32- 

Lemma 5.5 Assume that (J5.115P holds, < C+ < "^3, IK'it' — W*,z — Z*)\\ag < ^^j32 and 

\{{w-W*)^,{z-Z*)^){^,^P)\<N^\ {if,^) G (0,C+) X (0,m) 

with a constant N > 1, where {W*,Z*) is the unique weak solution to the problem ()5.3ip ~ ()5.36p in 
Lemma \5.1\ Then for any {f,ip) G (0, C+) x (0, ?n), 

m-W*),dZ-Zn,){^M < (^^^%^m^.r^^'-'^' + ^^N + ^,2y, (5.116) 

where {W,Z) is the unique weak solution to the problem (|5.19p - ()5.24p in Proposition 15. Jl and /Lii2 
depends only on 7, ?n, A, 5i, ^2, (^3, /^i; A*2 ctnd j32- 

Proof. Let (#', i2°) be the weak solution to the problem (j5.73p - ()5.78p in the proof of Proposition 
O Set 

(C/,y)((^,V) = (#'^,^^)(<^,V'), (¥',V') € (0,C+) X (0,m). 

Then, ([/, 1/) solves the following problem 

U^ + b'l\Q)U^ = h~\Q)p{Q)(2Q^ + W)[U + F^ +F^ + F^), (^,^) G (0,C+) x (0,m), 

V^^ - h^'\Q)V^ = \b-\Q)p{Q){2Q^ _ z)(F + G^ + G^ + G^) , {^, V') G (0, C+) x (0, m), 

?7(0,V)=0, i:e{Q,m), 

[/((/;, o) + y(<^,o) = o, v'e(o,c+), 

C/(¥P,m) + y((^,m)=0, VpG(0,C+), 

where 
F^ =(2Q<^ + W)-^W^{W + (VF* + Z*) + {z- Z*)), 
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F^ ={2Q^ + W)-'w(b{Q)p-\Q){b-HQ)p{Q)) - h-\Q)piQ)Q^ [W + {W* + Z*) + {z - Z*)) 



F^ ={2Q^ + W)-^W{z - Z*)^, 

6) (h-^( 6)7,(6)] 

'f 2 

+ i2Q^ + W)-^WiW* + Z*)^, 

G^=- (2Q^ - Z)-^Z^{{w - W*) + {W* + Z*) + 3^), 
G^ = - (2Q^ - Z)-^Z{w - W*)^, 

G' ={2Q^ - Z)-^z{h-\Q)p{Q)Q^ - h{Q)p-\Q){b-\Q)p{Q)) ^) {{w - W*) + {W* + Z*) + ^) 

-{2Q^-Zy^Z{W* + Z*)^. 

Direct calculations show that 

zpi 
\{F^G'){^,^P)\ < ^'^\ Nip\ {^,ij) € (0,C+) X (0,m) 

and 

\{F\G'){^,i^)\ < lii2^\ {^,ip) G (0,C+) X (0,m), 
where ^12 depends only on 7, m, A, (5i, ^2, ^s, fii, ^2 and /32. Then, Remark 15.31 yields 

n 

Remark 5.4 Under the assumptions of Lemma \5.5[ if, in addition, 

2^8/^2+^4 n A/4-1/2 , , 2/3l /CIIVN 

^^^mi^4r3 + /.12 < ^^^^^^^A^> (5-117) 

t/ien (15.116P implies 

\{{W -W%,{Z - Z*)^){^,^P)\ < N^\ {ip,ij) G (0,C+) X (0,m). 

Therefore, replacing ^ with 

J =[{w,z) G L°^((0,C+) X (0,m)) x L°°((0,C+) x (0,m)) : \\{w -W*,z- Z*)\y < fisp2, 

\{{w -W*)^,{z - Z*)^){v,i;)\ < N^\ ((/p,V)G(0,C+)x(0,m)} 

in the proof of Proposition \5.S\. we can show that 

Proposition 5.3 Let the assumptions of Proposition \5.S\ hold and {W, Z, Q) be the unique weak so- 
lution to the problem (j5.7p - (|5.14p . Assume further that < C+ !^ '^s? -^ ^ 1 and (J5.117P holds. Then 
{W, Z) G J and 

\Q^A^,n < (A^ + /Uio)v'\ \QM^^^)\ < W3(A^ + /^io)^'^/'+'/', {V,^P) G (0,C+) X (0,m), 

(5.118) 

where /iis depends only on 7, m, X, 61, 62, /^i, /U2 and (32- 
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Remark 5.5 Under the assumptions of Proposition [5731 if, in addition, 

A>■.„{(A^ + „„)(£^^)^^^,,(A' + ™)(iif±i)""*"^}. (5.119) 

then (|5.118|) leads to 

\Q^A^,n < I3J _ /^L— Y^%\ (if,^) G (0,C+) X (0,m), 



IQ^^IV^.V')! < /34( _ ^^ )'^^'^'^%^"/^+^/^ (VJ, V) G (0,C+) X (0,7n) 



We now choose 

_/ /32 \4/{A-2)/ c* X (5A+6)/(A-2) 

^^ ~ V/iii(l + /^i3)(2/i8/32 + M4)(^5/32 + 2f3i)) Uy^fTT^ 

A = (1 + ..,) (2«„ + ^i^<i55|±^ + i) (£4±i)"'"'^ 

and 

Pi 
Then, ([511711 and (fSHOD hold. Set 

•^0 = _ .... , minln, T2, ra} 
cV(52 + 1 

and 

^ = |(^+,^) G C([0,/+]) X C^'H[0,1] X [0,m]) : ^+ satisfies dSI]), while ^ satisfies 
dSSD-dSSSD with dOOjl . |^^0((/.,V)| < /34(cJ+)^+2^^ and 

Replacing 5o and =5^ with 5o and ^ in the proof of Theorem 15.11 one can prove that 

Theorem 5.2 Assume that /+ G C^{[0,1+]) satisfies ^JT\i and (J3TT]) . If < 1+ < 5o, then the 
problem (|3.47|) - (|3.52|) admits at least a solution Q G C"^'^([0,C+] x [0,m]), which satisfies (|5.102p and 

((^,V)G(0,C+)x(0,?n), (5.120) 

where 6o and //14 depend only on 7, m, A, (5i, 82 and 83. 
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5.5 Uniqueness of sonic-supersonic flows 

In this subsection, we study the uniqueness of sonic-supersonic flows. Owing to the strong singularity 
of the equation at the sonic curve, one can only prove the uniqueness of sonic-supersonic flows with 
precise estimates as in the existence (Theorems 15.11 and 15. 2p . In Theorem 15.21 we get a C^'^ sonic- 
supersonic flow under the assumption that /+ € C^([0,/+]) satisfles (j3.4ip and (j3.7ip . This smooth 
sonic-supersonic flow will be shown to be unique if /+ G C^([0, Z+J) additionally. Furthermore, it 
is also unique in the space of weak solutions given in Theorem 15.11 It is noted that the additional 
condition /+ G C ([0,^+]) arises from the nonlocal and implicit boundary condition (j3.5ip . 

Theorem 5.3 Assume that /+ G C"^([0,Z+]) satisfies ([OT]) and ([3TT]) . Then the problem ^JT}- 
()3.52p admits at most one weak solution Q G C'^'"'^([0, C+] x [0,m]) satisfying 

-M2^^+^ < Qi^, ^) < -Mi(/p^+2, ((^, ^) G (0, C+) X (0, m), (5.121) 

-M2ip^+^<Q^{ip,^)<-M^ip^+\ \Q4ip,iP)\<M2ip''^/^+\ ((^,V')G(0,C+)x(0,m), (5.122) 

where Mi and M2 (Mi < M2) are positive constants. 

Proof. Note that the problem (I3.47p - (|3.52p is equivalent to the problem (j3.16p - (|3.2ip . which has 
standard initial and boundary value conditions, and (j3.16p is strictly hyperbolic away from the sonic 
curve. So, it suffices to prove the uniqueness of weak solutions to the problem (j3.47p - (|3.52p for small 
C+ > 0. Thanks to Theorem [521 for sufficiently small C+ G (0,1), the problem (f3:i7p - (l332l) admits 
a solution Q G C^'^([0, C_|_] x [0, ?7i]) satisfying (j5.102p and (|5.120p . Therefore, one need only to prove 
that Q is the unique weak solution to the problem (I3.47p - (l3.52p . Assume that Q G C'^'^([0, C+] x [0, m]) 
with (|5.12ip and (|5.122p is a weak solution to the problem (f3:?7l) - ([332]) . Note that 

div(<5ip, -b{Q)Q^) = div{Q^, -Q^) = 

in (0, C+) X (0, m) in the sense of distribution. It follows from the theory of L°° divergence-measure 
vector fields ([l]), (j5.12ip and (I5.122P that Q^p and Q^ have L°° trace on [0, C+] x {ip} and {(/?} x [0, m] 
for each '0 G [0,m] and Lp G [0, C+]. 

Denote the functions given by (J3.45P corresponding to Q and Q by X_|_ and X_|_, respectively. For 
convenience, we use /i > and M > in the proof to denote generic positive constants depending 
only on 7, m. A, /+, 5i, 82, 6^, ||/| ||L-=((o,i+)), ^1 and M2. Set 

Q{^, i;) = Q{ip, -0) - Qiip, ij), {if, V) G [0, C+] X [0, m] 

and 

X{ip) = l+(v9) - X+(99), (v., V) G [o,C+]. 

Then Q G C70'i([0,C+] x [0,m]) and X G C^{[0,C+]) satisfy 

biQi^, V)) - b{Qi^, V')) = hiif, i;)Q{ip, ^p), (<^, ^) G (0, C+) x (0, m). 



and 



Q^{ip,m) = ei{ip)X{ip) +e2{ip)Q{(p,m), ip£{0,C+) (5.123) 

X\ip)=di{ip)X{ip)+d2iip)Q{ip,m), (^G(0,C+), (5.124) 
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where 

h{^, ij) = - f K"{tQ{ip, V') + (1 - t)Q{^, ij))dt, {if, ij) e (0, C+) X (0, m), 
Jo 

eiM = -— ^7 ,,,L,., ;- / Ei{tX+{^) + {l-t)X+{^))dt, ^G(0,C+), 

b{Q{Lp,m)))A^ {Q{ip,m)) Jo 

/■I 

e2{ip) = -Q'+{X+{ip))cose+iX+i^)) / E2(tQ(v',m) + (l-i)Q((/.,m))dt, (^G(0,C+), 

JO 

^i('^) = 7^T7^ ^/ Di{tX+{^) + {l-t)X+{^))dt, ^G(0,C+), 

^+ (<9(¥?,m)) Jo 

d2{ip) = cose+{X+{ip)) D2{tQ{ip,m) + {l-t)Q{ip,m))dt, ipe{0,C+) 



with 
and 



^i(s) = (e+(s)cose+(s))', Di(s) = (cose+(s))', 0<s</h 



Making use of the conditions (j3.4ip . (j3.7ip and /+ G C^([0,/+]) on /+ and the asymptotic behavior 
(J5.1(J2|) and (j5.12Up on Q and (j5.12ip and (|5.122p on Q near (p = 0, one can show by direct calculations 
that for {ip,ip) G (0,C+) x (0, m), 

and for 99 G (0, C+), 

|ei((^)| < iVVV2, |e2((^)| < ivV/2-1^ 1^/(^)1 < Af^A/2+1^ 

|e'2((^)| < AV/2-2^ |^^(^)| < aVA+1^ |^^(^)| < M(^-V2-i^ 
-M¥.3A/2+i < Q^i^^^rn) < -^(^3A/2+i^ |Q^^((^,m)| < M(^3A/2^ 

Here, the bound of e'^ depends on ||/+ ||Loo((o,ii )) and it is noted that this is the only reason for the 
assumption that /+ G C^([0,/+]). 

First we estimate X. It follows from (I5.124P and the regularity estimates of di and d2 that 



which implies 



and 



\XicP)\ < M / ip-^/^-'\Q{ip,m)\dip, (t> G (0,C+), (5.125) 

Jo 



X2(0)<M(^2/ ^-^-^Q-i[^^m)d^^ 0G(O,C+) (5.126) 

Jo 



|X'((/))|2 <M0^^+M ^~^-'^Q'^{^,m)dip + M<j)-^~^Q'^{^,m), 0g(O,C+). (5.127) 

Jo 
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The Holder inequality gives 

- \Q{cl>,i,)\di^ + |Q^,(,/),^)|d^ 

"I Jq Jq 

, I r<t> rm rm. 

^fn Jo Jo Jo 

<Mcj)^+^ / ip~^-^Ql{ip,ij)dipdij + M Ql{cP,iP)diP, 0G(O,C+). (5.128) 
Jo Jo Jo 

Then, one gets from (I5.126l) - (l5.128p that 

pep pm pep pm 

X\<l>)<McP' / (^-^-2Q2(<^^^)rf^rf^ + M,^A/2+2 / / <^-3A/2-3g2(^^^)^^^ 0G(O,C+) 

Jo Jo Jo Jo 

(5.129) 



and 



Jo Jo Jo Jo 



+ Mr^-^ Ql{<t),i^)di:, 0E(O,C+). (5.130) 

Jo 

Next, we improve the estimate (15.128P for ||Q(<^, •)llL°°((o,m))- It follows from the definition of 
weak solutions that 

d f^ [^ f ~ dQ - dQ \ 

— Q{(l),ilj)dilj = \^b{Q{ip,m)) — {ip,m)-b{Q{ip,rn)) — {ip,m)jdip 

{e3{ip)X{ip) + e4{ip)Q{ip, m))dip, (/) E (0, C+), 





where 



^3(<^) = -i^^TTX} ^ / Eiitx+i^) + (1 - t)x+{^))dt, ^ e (0, c+), 

A^\Q{ip,m)) Jo 

/■I 

e4(y^) = -eV(X+((^))cose+(X+((^)) / D2{tQ{ip,m) + {l-t)Q{ip,m))dt, (/9G(0,C+) 

Jo 

with 

\e3{^)\ < M^^-\ |e4((/9)|<M(^V2-i^ <^e(0,C+). 

Therefore, 

r4> 
Q{(t>,il^)dil^ <4> \e3{(f)X{(p) +e4{(f)Q{(p,m)\d(p 
Jo 

V2 _,w,^w,/ f^ ^o, . . \l/2 







<M,/)^+V2M x\^)d^)' +M0V2+1/2M Q^(^^^rn)d^y ^ 0e(O,C+), 



which, together with (15.1291) and (|5.128p . implies 

9 rd) rm. 

'0 Jo 
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(^(j) pin 

Jo Jo 



Thus 



snpQ\cP,-)<(- Q(</.,^)(i^ + / |gv'(</.,V')|dV' 



(0,m) 



(f) rm 



JO 

m 



<f> pm 

Jo 



Jo 
Moreover, it follows from (|5.125p and (j5.13ip that 



(5.131) 



X2(<^) <M(/>V2+2 / ^-3A/2-3q2(^^^)^^ 

Jo 

A+5 / / ,„-A-2^2 ^,„ „/,^^,„^„;, I ^^^A/2+2 / / ,„-3A/2-3^2 



<M4) 



^0 Jo Jo 



(5.132) 



which improves the estimate (15.129p . 

Now we are ready to prove the theorem by a weighted energy estimate. Let <j) G (0, C+] be given. 
Multiplying the equation of Q by (p~''^~^Q^ and —^~''^~^Qtp, respectively, and then integrating over 
(0,0) X (0,r7i), one can get that 

/•(p rm 







-A-1, 



f ^ ^Q^^Q^dipdip - ip {b{Q)Q^)^Q^d(pdii 

JO Jo Jo 

-t rm \ j^ \ r4> f-m j-4) 

~'~'^''^ +^/ / ^-""-^Qldipdi; - ip~^~^biQ)Q^Q^d^ 

'P=4> ^ Jo Jo Jo 







V~^~^Qld^ 



1 



^-^-%Q)Qld'il: 



1 />(/> i-m 



(^=0 z 



Jo 



{^-^''b{Q))^Q'^,dipdtP 



and 



m ^ p(f> pm 

if^^^^Q^^Q^difdip + / Lp~^~^{h{Q)Q.^)^Q^dipdii. 

JO Jo Jo 



(5.133) 
(5.134) 



It follows from the definition of weak solutions and a standard limit process that 



=\ I ^-"--^Qldi^ 



ip=<j, 



\ I 1 r4> rm. 
+ -^ I I if-^-^Qldipd^ - I ip-^-^b{Q)Q^Q^dip 



Jo 







1 r™ 

+ 2/ ^~^~'biQ)Qld^ 



1 pij) rm 



ip=4> 2 



JO 



(^-'^-^6(Q))^g^(i^dV- 



ip=m 

(5.135) 



and 

0: 



-A-l 



V QipQu^dip 



(A + 1)/ / ^' 
v>=<p Jo Jo 



^ '^Q^Qipdifdip + / / y' '^ ^Q^^Q^dipdip 
Jo Jo 



Jo 



4> rm 



i}=rn. Jo Jo 

A direct calculation and a standard limit process show that 



(f b{Q)Q^^Q^,dfd'ilj. 



^0 



f~^''b{Q)Q^Q^dip 



^=m Jq 



ip-^-'b{Q)Q^Q^di; 



ip=4, 



(f> j-m r<j) /""^ 

/ {ip~^~^b{Q))^Q,^Q,^dipd'il) + / / ip~^~^b{Q)Q^^,Q^dipd'tp. 
Jo Jo Jo 



Integrating by parts leads to 



if ^ ^hQQ^^Q^dipdijj 



'0 Jo 
-J ^-^-'hQQldi; 



-1 i-cp j-m 



ip^4' ^ 



JO 



{^~^~'hQ)^Qi^,d^di;. 



Summing up from (j5.133p to ()5.138|) yields 

= /i+/2 + Ji + J2, 
where 



A + 1 f^ r 



\ 4_ 1 ffp f-m 



JO 



-A-2//S2 



Q'+Q'-2Q^Q^]dipdi, 



Jo 



ip-^-'Q'd^d^p, 



(5.136) 



(5.137) 



(5.138) 



-1 /-(p rm 



-A-1 



A-1, 



Jo 

1 ptf /""^ 



(v5-'^-^6(Q))^Q^ + {v-^-'biQ))^Q'^ - 2{if-^-%Q))^Q^QAd^d^ 



Jo 

4> rm 



{^'^''hQ)^Q'^,dvd^l} 



'0 JO 
2/0 Jo 



if ( b{Q)Q^^Q^, - b{Q)Q^^Q^, + hQQ^^Q^ ) dtpdip 



{^~^~^b{Q))^Qldipd^ 

1 /-(f) rm 



1 r<f> rm 



JO 



<^ h{Q^ + Q^)Q^Q-^dipdip 



JO 



{(f h)^{Q^ + Q^)QQ^d(pdip 



(f) rm 



JO 



-A-1 



^ hQ^^QQ^dipdijj, 



1 r^n 

2 



Ji=^ I V-^-HQl + Ql-2Q^Q^]di; 







ip=<f> 



+ 



"P 



b{Q)Q'^ - 2b{Q)Q^Q^ + 6(Q)Q^ + hQQ'^ ) d^ 
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tf=cj> 



^ 



-^-'Qld^P 



1 



ip-^-^b{Q)Qi^di; 



ip=(j) 



^ 



-A-l 



KQi^ + Qip)QQ-,pdip 



ip=(j> 



and 



J2 = 



v? 



-A- 



^ [b{Q)Q<pQ^ - KQ)Qi>Qv + KQ)QvQ^ - b{Q)Q>pQ^) d(p 



ip=m 



f 



-A-1 



biQ)QipQip - b{Q)Q^Q^jdip 



•^=171 







9? hQjpQQ^dif 



= - \^~^-^hQ^Q^ 



1 



^-^-^b{Q)Q^Q^d^ 



+ - I {ip-^-^hQ^)^Q^d^ 



ip=m 



ip=m Jq 



ip-^-^b{Q)Q^Q^,d^ 



ip=m 



Therefore, 



^-""-^Qldipd^ 



/o Jo 

-1 /•lb rm 

^ I I -A-1, 



1 /-(p rm 

2 



{^~^~H{Q))^Qld^d^ 



JO 

1 f'" 
+ ,, ^~^~'b{Q)Qldi; 
V='t> ^ Jo 

1 r<f) rm 



V=<1 



l^-^-^hQ^Q^ 



{(p,ilj)={(f>,m) 



if ^ ^h{Q^ + Q^,)Q^Q^,d^dil) + 

^ '0 ^0 



2 -/o Jo 

+ / / ^-^-^hQ^^QQ^d^di) - - / '^^^1 

JO JO ^ JO 






ip=<j) 



1 



((^~^~^/iQ^)^Q2^(/j 



+ / ^-^-^b{Q)Q^Q^d^ 

4'=m Jo 



4'=m 



It follows from the regularity estimates of Q, Q and h, the asymptotic behavior of Q^ and Q<^^ on 
the wall near 93 = and the Holder inequality that 



-3A/2-3g2|^^^^^ 



/o Jo ^ 

<M / ip-^-^\Q^Q^\dipdil: + M I I if 
Jo Jo Jo Jo 

rm 







"P 



■0 pm f<p rm 



-A-4 



\QQ^\d^pdip + M I I LP 
Jo Jo 






+ M I ip-^~^ 
Jo 

<j> cm 

3A/4-5/2^2 



\QQ^\dil) 



ip=m 



<M I I ip-'-'^/'^-^/'^Q^^dipdiP + M 



Jo 

+ M 



+ M / ip~^~^Q^dip + / ip~^~^b{Q)Q^Q^dLp 

ip=<f> Jq i!=m Jq 

(p pm p(j) pm 

/ ^~^>^l^-y^Qld^d^ + M I I (p-5V4-7/2g2^^^^ 

Jo 



Jo 



10 



if-^-'-^Q^di) 



V=<P 



+ M / if-^-^Qldij 



v=<t> 



+ M / if-^-^Q^d^ + / ip~^-\Q)Q^Q^dip 

Jo ^=^ Jo 



ip=m 



(5.139) 



To estimate the last term on the right side of (|5.139p . one substitutes (j5.123p into this term, then 
integrates by parts and uses the Holder inequality and the regularity estimates of Q, ei and 62 on the 
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wall near ip = io get 



^=m 



LP b{Q{(f,m)){ei{(p)X{(f) + e2{(p)Q{<f,m))Q^{ip,m)d(f 



1 



A-l, 



0"-'-^6(g(0,m))ei((/.)X(0)Q(</.,m) + -0-^-^6(Q(0,m))e2(0)Q'(</>,m 



{ip~^~^b{Q{ip,m))ei{ip)yX{ip)Q{ip,m)dip - / ip~^~^b{Q{ip,m))ei{ip)X'{ip)Q{ip,m)dip 
Jo 



1 



ip b{Q{ip,m))e2{(p)) Q {(p,m)d(p 



<M(l)-'^\X{(t))Q{(t),m)\ +M(I)' 



-^~^Q\(l),m) + M / ip~^^'^\X{ip)Q{^,m)\d^ 
Jo 



+ M ip~^\X'{ip)Q{ip,m)\d^ + M v3~^~*Q^((^,m)d(/? 
JO Jo 

r<t> c<t> 

<M(l)^~^X'^{(t)) + M<p"^~'^Q'^{(l), m) + M ip~^/^~'^X^{ip)dip + M ip~^^/'^~^Q'^{ip, m)dip 

Jo Jo 



+ M (p^\X'{^)\'d(p + M ^-^-''Q\ip,m)d<p, 
Jo Jo 



which, together with (|5.130l) - (l5.132p . yields 
ip~^^^b{Q)Q^Q^dip 



(f) /•m 



JO 



<M(0^+2 + 0V2+5) I I ^-A-2g2 ^^^^ ^ M(0^/2-l ^ ^. ^ ^ ^ ^ 



<f> pm 
JO 



+ M0^/2-l / ^-3A/2-2g2^^ 



ip=4, 



(5.140) 



Substituting (j5.140p into (|5.139p and using (j5.13ip . we obtain 



<M(^^/^-i/2 / / ^p~>^~-^Q^^dipdip + M 
Jo Jo 







-1/2 



<f=<P 



(f) pm 



Jo JO 



"P 



-3A/2-3^2 



Q^dipdip 



+ Af</.^/2-l / (^-3A/2-2g2^^ 

io 



'P=<A 



Owing to A > 2, one can conclude that for sufficiently small (j) G (0, C+], 

Q^{ip,7p) = 0, Q^{ip,ip)=0, {ip,7p) e{0,(l))x {0,m). 

Therefore, 

Q{ip, i/j) = Q(99, V'), (v?, ^) G (0, </>) X (0, m). 
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